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1. Introduction 

The unification of tlie tlieory of critical phenomena and Euclidean quantum field theory 
has been immensely yielding for both areas. The theory of critical phenomena improved in 
a very important way our understanding of renormalization. The concepts of universality 
classes, critical points, marginal operators etc. have made many aspects of the renormal- 
ization procedure and the renormalization group equations much more intuitive for the 
field theorists, while the full machinery of Green functions and Feynman integrals have 
been a most important technical ingredient, which has influenced the theory of critical 
phenomena. 

There is a special area of field theory which is most intimately linked to statistical 
mechanics, but where it nevertheless was felt that one could not use the methods from 
statistical mechanics: the quantization of geometrical objects, i.e. the first quantization 
of the free relativistic particle, the relativistic string and the quantization of gravity. In 
all cases the classical action is defined entirely in terms of simple geometrical expressions 
which, however, when written in terms of an explicit parametrization, become ugly and 
difficult to treat by analytical methods: the length of a curve will involve the square root 
of the coordinates, just to take an example. Having chosen a parametrization one gets 
the additional problem that the final result should not depend on this parametrization. 
Rather absurdly, this is often viewed as a virtue: the theory has a large invariance, namely 
invariance under diffeomorphism. In reality it reflects our inability to quantize the phys- 
ical, geometric degrees of freedom. Much is lost compared to the beauty present if we 
simply view the system as a statistical ensemble of geometrical objects, the partition func- 
tion being defined as the integral over all such objects, the weight being the exponential 
of the classical action which itself is defined entirely in terms of the geometry. 

When the link is made between statistical systems and Euclidean quantum field theory, 
the discretization of space is usually a key ingredient. By discretizing ordinary space and 
restricting the volume to be finite we approximate the field theoretical problem by a finite 
dimensional problem, which can often be viewed as a generalized lattice spin problem. In 
the infinite volume limit we can look for phase transitions of the spin system, and if these 
are characterized by a divergent correlation length, it is possible to forget the underlying 
lattice structure: we can take the lattice spacing to zero compared with the correlation 
length, and in this limit we might recover a continuum field theory. Usually masses and 
coupling constants are defined not at the critical point but by the approach to the critical 
point. In this process space is just playing a spectator role. We do not demand any local 
invariance maintained and this is why we can discretized space without any problems. 
Lattice operators which break explicitly Euclidean invariance will be suppressed in the 
scaling limit. Local invariance usually mixes high an low frequencies and is much more 



difficult to discretize in a simple way. This is often used as an argument against attempts 
to discretize geometrical theories and it is presumably correct that it is rather fruitless 
to attempt a discretization of a given parametrization of the geometry. If we, on the 
contrary, choose to discretize the geometry directly there will be no problems at all. 
The main theme of these lectures will be that a discretization of geometry is natural, 
it fits perfectly with a statistical mechanics interpretation of the theories and the whole 
machinery of critical phenomena and scaling limits can be applied in a very powerful way 
to geometrical theories. In addition we will get a pleasant surprise: some of the theories 
are easier to solve directly at the discretized level than in the continuum, and this fact 
allows us to study the scaling limit in considerable detail. 

What is meant by "geometrical theories"? These are theories which describe the 
propagation of ((i-l)-dimensional manifolds by summing over an appropriate class of d- 
dimensional manifolds of which they are boundaries. The action might depend only on 
the intrinsic geometry of the (i-dimensional manifolds. If the manifolds are embedded in 
R^ the action might in addition depend on the extrinsic geometry. 

The simplest of such theories describes the propagation of point particles, i.e. we 
consider the theory of paths P{x, y) between two space points x and y which belong to 
R^ . The simplest action for the paths will be 



S[Pix,y)] = m dl + X dl\k\ + ---. (1.1) 

Jp{x,y) Jp{x,y) 

The ffist term only refers to an "intrinsic" property of the path, its length, which can be 
defined without reference to the target space where it is embedded. The second term, a 
curvature term, refers explicitly to the embedding. In principle we can add higher powers 
of the curvature and also torsion terms to the action. 

If we move up one dimension we get a theory which describes the propagation of one 
dimensional objects (strings), i.e. we consider the theory of surfaces S{li) spanned between 
boundary strings (or loops) /», i = 1, ..,n. The simplest action for such surfaces will be: 

S\Mik)]=iil dA(M)+xf dA(M)H^- (1.2) 

JMih) JMih) 

where M{li) denotes a manifold with boundaries /j and dA{M) the area element of the 
manifold. Here a number of interpretations are possible. We can choose to view the 
area as the induced area of a surface embedded in, say, R^ . This is analogous to the 
interpretation given in ( p, . 1| ) . Then it is possible to add further terms depending on the 
extrinsic geometry as indicated in eq. ( |1.2| ) where H denotes the extrinsic curvature 
term. In two dimensions we have in addition the possibility to formulate the theory 
without reference to a target space (we have the same possibility in one dimension, but 
the theory will be trivial unless matter fields are added, as will be discussed later). By 
viewing the variables as intrinsic we consider two-dimensional quantum gravity and the 
action can be written as 

S[M{U)] = f d'^y^ U - MR + X^R^ ■■). (1.3) 

In this formula Qab denotes an internal metric. The ffist term is still the area term now 
written in terms of the internal metric. R denotes the intrinsic or Gaussian curvature of 



the manifold. 

When we finally move to higher dimensions we usually have no interest in referring 
to some embedding space and we will discuss the propagation of {d — l)-dimensional 
manifolds 6, via (i-dimensional manifolds which have the 6j's as their boundary entirely 
in terms of intrinsic variables. The action will be 

S[M{k)] = f d'^^v^ (/i - Aii? + A2i?2 ■ • •) + b.t., (1.4) 

where b.t. denotes boundary terms. In a formula like (|1.4| ) the natural variables to 
consider are equivalence classes [gab] of metrics. In a continuum formulation it is quite 
difficult to work directly with such variables. One almost inevitably ends up with gab 
themselves. Somewhat surprisingly, the regularized (discretized) quantum theory offers 
the possibihty to work directly with equivalence classes. 

The natural way to define the quantum theory corresponding to the above classical 
actions is to use the path integral. (First) quantization tells us that the propagator 
G{bi, ..., bn) between d — 1 dimensional boimdaries will be 

G{k) = f dM{bi) e-^[*^(^»)], (1.5) 



where the summation is over d-dimensional geometries. From this point of view eq. ( |1.5| ) 
defines the theory of fiuctuating geometries. We have already defined the action. One 
has to contribute a meaning to the integration over geometries. A key ingredient in doing 
so will be to approximate in a natural way the smooth structures by piecewise linear 
structures. In this way ( p. . 5|) will be an ordinary statistical system. It is possible to 



discuss the critical properties of this system and they agree with the results obtained by 



continuum methods whenever they are known. In addition ( |1.5|) will provide us with a 
nonperturbative definition of the theory in cases were continuum methods seem powerless. 

In these lectures the intrinsic properties of our geometric objects will be described by 
the metric. However, as will be clear, this description can be replaced by any other, 
using more appropriate variables, if needed. The only requirement seems to be that these 
variables have a natural description on piecewise linear structures. 

In the following I will try systematically to develop the quantum theory of geometric 
objects starting from the simplest one-dimensional objects in (|1 . 1|) and ending with the 
four-dimensional objects in eq. (|1.5|) relevant to quantum gravity. 



2. Bosonic propagators and random paths 

2.1. Quantization 

The classical action of the free relativistic particle in R^ moving from x to y is, as already 
mentioned, expressed by: 

S[Pix,y)] = mo f d/, (2.1) 

JP(x,y) 



The classical equations of motion are derived by choosing a parametrization of eq. (pTl 



x(0 : [0, 1] ^ i?^, x{0) = x, x{l)=y. (2.2) 

With this parametrization one gets: 

5 = mo/'deV^, (2.3) 

Jo 

^^ -i(x7|x|) = 0, (2.4) 



where x = dx/ d^. The obvious solution to the classical equation of motion is a; = const., 
i.e. the straight line between x and y, but any reparametrization of the straight line: 

x'iO ^ x'ifiOl /(O) = 0, /(I) = 1, / > 0, (2.5) 

is a solution too. This is a reflection of the reparametrization invariance of the geometrical 
action ( PTT| ) 

First quantization of the system is implemented via the path integral. We get the 
propagator G{x, y) by summing over all paths connecting x and y and weighted by e~'^[^l 
The quantum aspect comes precisely from the fact that not only the path which solves 
the classical equations of motion contributes to the propagation: 

G(x, y; mo) = / VP{x, y) q-^^p^^^^v)^ . (2.6) 



Each path should be counted only once in eq. ( p.6| ) . Reparametrizations like ( p.5|) should 
not be counted as different paths. 

In order to contribute a meaning to JT>P we have to introduce a cut-off. Note that 
while we usually think about smooth paths, the action is in fact defined on a larger class of 
paths, the ones which are only piecewise smooth. Let us introduce a cut-off by considering 
only piecewise linear paths where each step on the path is of length a. This implies that 
we refrain from discussing structures smaller than a. Note that the cut-off introduced in 
this way by definition is reparametrization invariant, since it refers directly to a length 
in R^ . The possible length of the paths will now be a multiple of a : I = na and the 
action of such a path will be S* = mo/. For each piecewise linear path between x and y 
consisting of n pieces of length a we have to integrate over the possible positions of the 
n — 1 interior points compatible with the length assignment a. If we denote the vectors 
of the n linear parts by aci, ii being a unit vector in R^, we have: 

/oo ^ n 
VP{x,y) ^ E / n de. ^(^E^. - (y-x)), (2.7) 

n=l "^ i=0 



Gaix,y;mo) = ^ e— o''" / J] de, S{aJ2e, - (y - x)). (2. 

n=l "^ i=0 



We can calculate the propagator at the discretized level by getting rid of the 5- function 
by a Fourier transformation: 



ra=0 



i=l 



Ga{p; mo) = J dx e-^P(--^) G{x, y- mo) = E e"™"-" J J] de, 
We have 

de e-'P-*^ : 



-lap- Si 



(2.9) 



Zn ' 



J p-i (pg) 



2 



/(pa), 



.(W2))- 
where 

/(O) = Vol(5^-i), /(pa) ^ /(0)(1 - c^(pa)2 + ■■■). 

The final expression for Ga(p; mo) is 

n=0 



oo 

G'a(p;mo) = E(e-'"°'^/(H)" = T 



- e-'^oa /(pa) ■ 



(2.10) 



(2.11) 



This is an exact expression for the regularized propagator. We have to take a — i> to get 
the continuum hmit and it is seen that we get the free relativistic propagator if we at the 
same time renormalize the bare mass niQ such that 



-rriQa 



f{pa) 



^2 2 2 



I.e. mo 



log/(0) , 2 2 

h c mpj^a 



From eqs. (|2.10| ) and ( |2 . 1 1| ) it follows that: 
1 1 1 



Ga{p;rno) 



c^ a^ p^ + rripjj 



c^ a^ 



Gc{p;mph), 



(2.12) 



(2.13) 



where Gc{p', rripi^) denotes the continuum propagator. The divergent factor which relates 
Ga and Gc is a kind of wave function renormalization, but it has a physical meaning 
since the power of a which appears reflects directly the short distance behavior of the 
propagator, as we shall discuss in detail later. 

It is worth rephrasing the above results in terms of dimensionless quantities, and in this 
way make the statistical mechanics aspect more visible. Introduce fi = rrioa and q = pa 
and consider coordinates in R^ as dimensionless. The steps in the random walk will be 
of length one and (|2.11|) reads 



G',(g) = E «"""/(# 



n=0 



e-M/(g) 



(2.14) 



/x acts like a chemical potential for inserting additional sections in the piecewise linear 
walk. We have a critical value fic = log/(0) of the chemical potential fi. For fi > fi^ the 
sum is convergent for all q and the average number of steps in the random walk is finite. 
For fi < fic the sum is divergent for some range of q and the average length of the paths 



not defined. For fi ^ fic very long paths will dominate the sum and this is where we can 
take a continuum limit. We can write 

/(g)~e^=(l-cV), G,{q)^ ^—^, (2.15) 

and at this point we can introduce the physical length scale a, and the physical momentum 
Pph and the physical mass mph, both intended to be kept fixed for a — ^ 0, by 

/i - /ic = mJha^ q = cpp^a. (2.16) 

This defines a as a function of fj, by 

a{fi) = y/JT^^c- (2-17) 

Later we will later discuss relations like (|2.16D in great detail. 



It should be emphasized that the critical behavior we have found this way is universal. 
Any "reasonable" class of random paths should result in the same scaling limit. The 
piecewise linear paths are convenient because the results, even at the discretized level, 
are Euclidean invariant. If we choose to regularize the summation over all paths by 
considering the sub-class of paths which can be formed by links on a hyper-cubic lattice 
in R^ , the regularized propagator will be: 

Ga = Ee-"^°'^"E'^(«Ee. - {y~x)). 

Again the regulatized propagator can be computed by (lattice) Fourier transformation and 
in the scaling limit one gets the same results as for the piecewise linear random walks. 

2.2. One-dimensional gravity 

Let us now turn to a somewhat different quantization of the free propagator. The action 
(2.1) has a beautiful geometrical interpretation and from the discretized point of view 
there was no problems associated with the quantization, as explained above. However, 
the square root, which appears in eq. ( p.3| ) after a choice of parametrization, makes it 
very difficult to use the action in formal continuum manipulations. For this reason it 
might be preferable to use a different action which at the classical level is equivalent with 
the one given in eq. (|2.1|): 



LX «/ (J 






(2,18) 



In eq. (|2.18|) x^{^) denotes a path in R^ such that 

a;^(0) = x'^, x^'{l)=yf', 

and gab is an internal metric on the one- dimensional manifold given by the parametrization 
^. The indices a, b can only take the value 1, but we have written the action in a general 



covariant way. Like ( |2.1| ) the action ( |2.18| ) is invariant under the reparametrization ^ 



/(^) defined in eq. ( |2.5| ) provided g satisfies 

9^^\fi0) = j-,9i0- (2.19) 

This transformation rule ensures the invariance of internal distance: 

ds' = gab dC de = gii^ dr df; a, 6 = 1. (2.20) 

The classical equations of motion, obtained by considering gab and x^ as independent 



variables, agree with the classical equations obtained from eq. ( pTl 




0. 



However, it is by no means obvious that the quantum systems defined from eqs. (|2.1| ) 
and ( ^.18| ) agree. The quantum equivalent of eq. ( |2.6[ ) is: 



where x(0) = x and x(l) = y. The integration variable is equivalence classes of metrics, 
i.e. metrics which are related by reparametrization. This is indicated by the symbolic 
division by the "volume" of the diffeomorphism group. 

In order to define the functional integral in eq. (|2.21| ) we first introduce a reparametriza- 
tion invariant cut-off "a", i.e. we consider only paths x(^) which have no structure below 
the length scale ds = a, ds given by eq. ( |2.20| ). This can be achieved by restricting 
ourselves to piecewise linear paths where the internal length of the individual pieces is a. 
For a given smooth metric gabiO ^^"^ ^ given smooth path we can, if we want, approxi- 
mate the action by a corresponding action of a piecewise linear path: first we discretize 
the manifold [0, 1] parametrized by ^ according the above prescription. Here it is impor- 
tant to note that the only reparametrization invariant quantity which characterizes the 
manifold is the total length: 



/ 



I' de v^ (2.22) 



The length / is clearly invariant, and for any metric gab satisfying ( p.22| ) we can transform 
it to the constant metric g^^ = I'^Sab by choosing the function / in ( |2.19[ ) as 



/(O = y // dr V^- (2.23) 

A given discretized path will always have the length na. For a given / we simply take 
n = [I /a], and given gab we can calculate the points C,i on the manifold parametrized by ^ 



where the internal length from the boundary is ia: 



la 



J ' de V^. 



(2.24) 



For a given metric we have discretized the manifold. The continuum action can now 
be approximated on the discretized manifold in a natural way by using that d^{C,i) ~ 

a/ J g{ii) and ^/gg"'^ = ^/ y/g in one dimension: 



de 



1 / dx^ 

V9\ de 



+ f^Vg 



n-1 
i=l 



(a;^(^^)_a;A«(^,^_^))^ 



+ /^ 



(2.25) 



Note that the discretized action contains no explicit reference to the metric gab- The ^j's 
appear as one would have expected it from a constant metric. But this is consistent with 
the observation that any metric is equivalent to a constant metric and the fact that the 
action reparametrization invariant. 

We can now combine the information given above: the integration of equivalence classes 
of metrics reduces in one dimension to an integration over the length / and in the dis- 
cretized approach this integration is replaced by a summation over n = I /a. The weight of 



(discretized) configurations will be determined by eq. (|2.25|) . Let us at this stage change 
to dimensionless quantities, as advocated above. By redefining x and /i we can dispose of 
a' and choose a = 1, and the discretized version of eq. ( p.21j ) will be 



n D 



G,{x,y) = E^-'^Vn n dxf e- EL. «-<-.) 

n=l j=l u=l 



(2.26) 



A few remarks about the formula: we have replaced the index ^j with i since it is just a 
dummy label of an integration variable. We have also defined xo = x and x„+i = y. Note 
finally that the index g in the measure "Dgcf) in formula ( p.21|) is explicit present in the 
measure in eq. ( p.26|) : it is n. 

Since the integrals in eq. ( |2.26|) are Gaussian they can be performed explicitly: 



n D 

n n dxr e 

1=1 fj,=i 



and we can write 



j:::'i<-<.rf 



TT 



nD/2 



-{x~y)'^ /in 



n 



D/2 



(2.27) 



G^,{x,y) = J2 



^{fiO-n)n-{x-yf /in'^ 



n=l 



n 



D/2 



fio = D log a/vt. 



(2.28) 



Again we see that there is a critical point /io above which G^{x, y) is convergent and below 
which it is divergent. 

If we take the Fourier transformation of G^{x, y) we immediately get ( |2.15| ) and we 
conclude that the theories defined by eqs. ( PTT| ) and ( |2.18D equivalent both at the classical 
and the quantum level. 

The reader is invited to compare the above calculation of G{x, y) defined by eq. ( |2.21| ) 
with a calculation performed entirely in the continuum, where one first has to introduce 
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parainetrization and ghosts and in the end a cut-off in order to calculate ill-defined deter- 
minants. The advantage of being able to work directly with equivalence classes of metrics 
should be obvious. 

2.3. Scaling relations 

The free relativistic quantum particle is described as the scaling limit of the random walk. 
In the following we will consider models which cannot be solved as completely as was the 
case above. In general we will not be interested in the complete solution at the discretized 
level, but only in the scaling limit where we approach the critical point. Let us therefore 
discuss scaling relations using the random walk as an example. The relations to be derived 
will be valid in a much broader context and we will use them many times in the rest of 
these lectures. 

Consider a model for random walks from point to a; in R^ . Step number n will be 
characterized by an initial position a;„_i, and an (unnormalized) probability distribution 
P{v) for a step to Xn = Xn-i + v. Let us assume that P is only a function of |f |. 
These assumptions can (and will be) considerable relaxed in the following, but they are 
convenient for a first discussion. If we use the notation xq = and x„+i = x for a random 
walk of n steps from to x, our generalized propagator can be written as: 

/n n+1 „ 

n dx, n ^(1^^ - 3;-il). / dt; P{v) = e^^ . (2.29) 

Note that / dx G^(x) can be performed since the integration over probability distributions 
reduces to a product of single integrals: 



x(/i) ^] dx G,ix) = ^_^_(^_^^) . (2.30) 

x(/i) is called the susceptibility lil. It is now clear that the critical point is fic- For fi > fic 
( |2.29D is convergent for all x. For fi < fic ( |2.29| ) is divergent for all x. One can solve 
the model given by eq. ( ^.291 ) in the scaling limit since it follows from the central limit 



theorem that the convolution of P many times reduces to the normal distribution, i.e. 
we get precisely the Gaussian model considered in the last subsection. Let us, however. 



discuss some general properties of eq. (|2.29| ) which will be of use for the more general 
models to be considered in the following. 

Theorem: G^{x) falls off exponentially for ^ > fi^- 

p{x) = G fj_{x) / xif^) is the probability density for a random walk from to x, i.e. Ap{x) = 
p{x)Ax is the probability for a random walk from to x and we have the the inequality 
Ap(x + y) > Ap{x)Ap{y), simply because the random walks from to x + y which pass 
through X are a subset of all random walks from to x + y. This means that — log Ap{x) 



^ This notation is borrowed from spin systems, where the spin susceptibiHty is the second derivative of 
the free energy with respect to the magnetic field, but also has the interpretation as the integral of the 
spin-spin correlation function over space. 
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is a sub-additive function of \x\: 

-logAp(x) < -logAp(ax) -logAj9((l - a)x), < a < 1. (2.31) 

From the sub-additivity alone it follows that the following limit exists 

lim ~^QS^/-(^) = ^(^)^ ^(^) > 0. (2.32) 

\x\^oo \x\ 

Since G^(x) by the definition ( f^.29| ) is a decreasing function of /i it follows that m{fi) 
must decrease as /i decreases towards ^c- Since / dxG^{x) exists for fj,> fic'we conclude 
that m(yu) > 0. It is now clear that m~^{fi) will serve as the correlation length for the 
random walk and we can only obtain a correlation length which is large compared to the 
individual step length provided m{fi) -^ for ^ -^ fie- We will assume this is the case and 
introduce the following critical exponents: 

m(yu) ^ {ji — ficT for fi — > fie- (2.33) 

Gu{x) ~ I in ,^ for 1< |x| < ^-. (2.34) 



dx G^(x) ~ — for \i -^ \Xc. (2.35) 

(/^ - l^cV 

The mass exponent v^ the anomalous scaling dimension ri and the susceptibility exponent 
7 are not independent. They satisfy Fischer's scaling relation: 

7 = z/(2-r/). (2.36) 

The proof of this relation is simple. From the behavior assumed for G^{x) in eqs. ( p.33| )- 
( pj.35| ) we can cut off the integration over x in the susceptibility at l/m{fj,), i.e. eq. ( |2.35| ) 
can be estimated as follows: 



x(/i) ~ / dx 



\x\<l/m{t,) X^-2+r, m(/i)2-'7' 

Eq. (|2.36| ) follows from the definition of 7 and u. 

For our simple random walk case it follows from eq. ( p.30| ) that 7 = 1 and it follows 
from the results in the last subsection that u = 1/2. This implies that rj = and this is 
why it is called the anomalous scaling dimension: If rj is different from zero it is anomalous 
with respect to the free particle. 

The exponent u has a direct geometric meaning: 

dH = l/i^: (2.37) 

where dn denotes the Hausdorff dimension of the random walk. It can be defined in the 
following way: Let (L)^ denote the average length of a path from to x in the ensemble 
of paths defined by eq. (|2.29| ). This is illustrated in fig. 1. Let us define 

12 








X 



Fig. 1. A typical piecewise linear path between and x. 
(L)^ ~ x'^" for |x| -^ oo, m(/i) |x| = const. 
Note that for any power k we can write for large n: 



(2.38) 



Let us use this ioi k = 1: 



CkU. 



(L). 



Xi-l 



\)m=lPi^^-^i-l) 



E„ e-/^" / nr=7 dx, nr=i ^(^^^ - x^-l 



(2.39) 



i.e. we can write: 



(L). ~ H. 



g^G^{X) 

G^{x) 



'm'{fi)\x\. 



(2.40) 



For ft > fic fixed and |a;| ^ oo we just get (L) ~ |x| and the reason is clear: (n)^ is 
finite for fixed \x\ and just goes to infinity proportional to \x\. This implies that only 
paths which are close to the straight line from to x contribute for /i > ^c ■ However, 
the coefficient m'{fi) diverges as fi ^ fic and we are interested in the limit where m{fi)\x\ 
is constant. This is the limit where we can introduce a physical "lattice" length a{fi), a 
physical length Xph and a physical mass r/iph by 



m{^) = m^]^a{jji), Xph = xa{fi), i.e. a{fi) ~ (/i — ficY 



(2.41) 



such that Xph and m-ph are kept fixed for /i —>■ fic, where a{fi) —>■ 0. In this limit we have 
(// — ficY ~ \x\~^ and 



W. 



m'{fi) 



1 



fi- He 



~ X 



\lv 



(2.42) 



Since the random walk described by ( [^.29| ) has v = 1/2 we arrive at the well known result 
that du = 2, but ( ^^.37] ) will be valid if z/ 7^ 1/2. We will meet such situations for the 
smooth random walks considered in the next subsection. 
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It should also be emphasized that the "physical" length of (L)^ will diverge in the 
scaling limit if we define Lph = La{fj,). It has to be this way for dimensional reasons: 

The propagator itself will be singular in the scaling limit. As already mentioned it can 
be viewed as a kind of wave function renormalization. From the assumed short distance 
behavior ( p.34| ) we get from ( p.4lD 



G^(x)~a^-2+''G,(a;ph,mph) for a -^ 0, (2.44) 

or, for the Fourier transformed where q = PphO.: 

G';.(g)~a''-'a(Pph,mph). (2.45) 

2.4- Smooth random walks 

In the last subsection we saw that ordinary random walks have Hausdorff dimension 
dn = 2. This result is quite universal as follows from the general expression ( |2.29| ) for a 
random walk. In the scaling limit the different random walk representations all agreed 



with the direct discretization of the action ( |2.1| ). Let us consider the first non-trivial 
generalization of purely geometrical nature: 

S[P{x,y)] = fi f dl + xf dl\K{l)\, (2.46) 

Jp(x,y) JP{x,y) 

where k{l) denotes the curvature of the path P. 

Recall how the curvature of a curve in R^ is defined. Let x^(^) be a parametrization 
of the curve. Use as ^ the length / of the curve. Let t'^(/) denote the normalized tangent 
of the curve: 

t^" = —- = i:^, tH^" = 1, tH^" = 0. (2.47) 

It follows that t^ is orthogonal to t^ and we write 

i^' = Kn'', (2.48) 

where k is the curvature and n^ the principal normal. It lies in the osculating plane of 
the curve. The geometry is shown in fig. 2. It is seen that 

|At| =2sin— ^ A(/., (2.49) 

A0 is called the angle of contingency and we have 

K{1) = ^, n{l) = ^, (2.50) 
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t(/ + dl) 




Fig. 2. The geometry related to the calculation of curvature. r{l) — 1/k{1). 




Fig. 3. The angle 0{ei, 62) between successive tangent vectors in the piecewise linear random walk. 



where r{l) is the radius of curvature. It is the radius of the osculating circle which is 
defined by the quadratic approximation to the curve and has its center on the principal 
normal at the distance r{l) from x^{l). 



To quantize the theory defined by the action ( p.46|) we have to perform the path 
integrala. Again it is useful to regularize the sum over all paths by restricting the sum 
to be over all piecewise linear paths, the length of the individual paths being a (which 
we choose equal 1) [2]. If the path consists of n linear pieces, the i's piece will be char- 
acterized by the a unit vector Cj. Although there is no universal definition of curvature 
for such paths, eq. (|2.49|) is a useful guide, where we simply take A0j to be difference 
^(ei+i,ei) in angles between e^ and Cj+i. This is shown in fig. 3. The discretized version 
of the action (|2.46| ) reads: 



n-l 



^[P„,/i,A]=/i^|e,| + A^2 



sm 



6'(ei+i,ej) 



j=i 



1=1 



^ There is an extensive literature on canonical quantization, saddle-point calculations and large d expan- 
sions of the theory. Ref. [ 1] is a very incomplete list. 
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n-1 



= /in + A^/(^(e,+i,e,)), (2.51) 

where we have aheady generahzed the curvature term to a function f{6) satisfying 

/(O) = 0, f{e) > for < ^ < TT. (2.52) 

The path integral for the propagator from to x can now be written as: 

G{x; i^,X)= f VP{x) e-^[^(")''^'^l (2.53) 



/n n—1 



In order to solve this model let us introduce the probability Kx{e2,ei) for a step 62 
provided the former step was ei: 



7^^(62,61) = -Le-me{e2,&,))^ Nx= f dei g-^^^^^^^'^i)) . (2.54) 

Nx J 

With this notation we can write: 

„ n n—1 

G{x; fx,X)-T. («"'' ^^) y n de^ H ^A(e.+i, e,) 5(^ e, - x). (2.55) 



i=l i=l 



We can calculate the susceptibility since the additional integration over x just removes of 
the (5-function in eq. ( p.55[ ) and the integrations of Cj become trivial (/ deiKx{e2, ei) = 1): 

Xi^i, X) = Jdx Gix; /i, A) ~ ^_J_^^^ . (2.56) 

For a fixed A it follows from eq. ( |2.56| ) that the critical point /ic(A) is determined by: 

e'^=W = iV,, i.e. x(/i,A) ^-. (2.57) 

/i-/ic(A) 

In the {fi, A)-coupling constant plane we have a critical line yUc(A), as shown in fig. 4. The 
model ( |2.55| ) is defined to the right of the critical line and the scaling limit for a fixed A 
leads to the susceptibility ( |2.57| ), i.e. the susceptibility exponent 7 = 1. From the general 
arguments presented in the last subsection it follows that G{x; /i, A) falls of exponentially 
for large |a;| for any (/i, A) to the right of the critical line: 

hm Il^^i^^^ = m(/.,A). (2.58) 

\x\^oo \x\ 

However, in this model we have in addition an independent correlation between tangents: 

/n 
Yl dci [in ■ K{en, Cn-i) " " " K{e2, ei)ei]. (2.59) 



j=i 
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A 



/^c(A) 




^^ 



Fig. 4. The phase diagram in the (/i, A)-plane. The theory is defined to the right of the critical fine /ic(A). 



For a fixed A the tangent-tangent correlation function falls of exponentially with the 
number of steps: 



(e„-ei) ~e-'^'"'(^), m(A) > 0. 



(2.60) 



This implies that we for a fixed A just have an ordinary random walk in the scaling limit. 
After no ~ l/mt{X) steps the initial orientation of the a tangent will be lost. If we group 
together tiq steps we will have an ordinary random walk where there is no correlation 
between successive groups of steps. Can we ever get any non-trivial behavior? The proof 
of ( |2.60| ) will show the way to a non-trivial random walk behavior. 

In order to prove eq. ( |2.60|) let us view Kx{e2, Ci) as a kernel for an operator Kx acting 
on L2(5^-i): 



(i^A0)(e2) 



dei Kx{e2,ei)(j){ei), 



eL\S 



D-l\ 



(2.61) 



The kernel is symmetric and Kx{e2,ei) > and Kx is compact. The Perron-Frobenius 
theorem tells us that the largest eigenvalue is non-degenerate and that the corresponding 
eigenfunction is the only one which can be chosen positive. Since the constant function 
clearly is an eigenfunction of Kx with eigenvalue 1 we conclude that 1 is the largest 
eigenvalue. It is easy to show that —1 is not an eigenvalue. Finally 0a(e) = a ■ e is an 
eigenfunction since 



(i^A0a)(e) = / dei Kx{e, ei) a ■ ii = a(A) a ■ e , 



(2.62) 



The equality follows since the integral is linear in a and invariant under simultaneous 
rotation of a and e, i.e. proportional to a ■ e. From Perron-Frobenius it follows that 
a{X) < 1, and direct calculation shows: 



< a(A) ^1 for A ^ oo. 



(2.63) 
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We can directly apply ( |2.62| ) in eq. ( |2.59| ) to get 



(e„ ■ ei) = a"(A) = e-™'^, m,(A) = log^VT- (2-64) 

The only possibility to get a non-trivial scaling limit is to take A ^ oo simultaneously 
with /x -^ /^c(A) since the tangent-tangent correlator can only approach macroscopic 
distances for A -^ oo. At the same time it is necessary that ^ -^ /^c(A) since this is the 
only possibility for the two point function to be non-trivial according to ( p. 5^ ). Such a 
scaling limit exists. In order to construct the limit it is useful to introduce the propagator 
G{ef,ei]x) which depends explicitly on the first and last step in the random walk: it is 
defined by ( p.55|) except that there is no integration over ei = Cj and e„ = e/. G{ef, ii] x) 



can be viewed as a kernel for an operator G on L'^{S^ ^) precisely as K\{e2,ei) is the 
kernel of K\. Let |1) denote the constant function: 1 = 1(e). We have by definition: 

G(x;/x,A) = (1|G|1) = f defdeiG{ef,ei;x). (2.65) 

By Fourier transformation we can get rid of the ^-function in ( [^.55[ ). If we introduce the 
notation Afi = /i — yUc(A) and denote the Fourier transformation of G(e/, e,; x), G{x) and 
G{x; /i, A) by G{ef, if, q), G{q) and G{q; /i, A), we can write: 

G(q) = y e-^'^" fe-'«-^ Kx] " = ^ ^ (2.66) 

In this formula g ■ e is viewed as a multiplication operator, i.e.: 

(/(g-e)0)(e) = /(g -6)0(6). (2.67) 

We have in analogy with (|2.65 ): 



G(g;/i,A) = (l|G'(g)|l) = J de^ de, G(e^, e,; g). (2.68) 

While these expressions look somewhat formal, they allow a rather transparent discus- 
sion of the scaling limit. Let us consider a fixed A and take A/i and q to zero. Consider 
the matrix element: 

(1|1 - e-^"-''^-^ kx\l) ^ (l|A/i + iq-e+hq- ef\l) = A/i + cgl (2.69) 

Let us now introduce the ordinary scaling for the random walk: 

A/i ~ mp^a^ g~ppha, a(/i. A) ~ ^/i - /ic(A). (2.70) 

For fixed A there is a finite mass gap from 1 to the next lowest eigenvalue. This implies 
that the only matrix element which contributes to (p.68|) to leading order in a is (|2.69| ): 



G(g;/i,A) = (l| ^ ^11) 
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1 



a2 



-^ TT + 0(^/a) 

Pph + "^ph 



(2.71) 



i.e. just the ordinary propagator. However, if A ^ cxo the mass gap goes to zero and more 
matrix elements will contribute to {l\G{q)\l) and the expansion in (2.71) is no longer 
valid. By Taylor expanding 0(e) around a fixed vector cq and using that e~^-^^^^^°''^^^ will 
be peaked around Cq one can show: 

i^^^e-^W^', c(A) ~ A-Mor A ^ oo. (2.72) 

where L^ denotes the Beltrami-Laplace operator on S^^^. Let us now introduce a lattice 
length scale by 

c(A) = Apha(A), m(/x, A) = mpha(A) = A/z. (2.73) 

The last equation fixes /i as a function of A and defines the approach to /ic(A) for A ^ oo. 
From (2.71) we get 

G(q-fi,X) = (l\ i ^11) 

^ (I L .. , ^ , . rll)- (2-74) 



a(A) AphL2 + mph + ipph ■ e 

This is our final expression. The matrix element on the rhs is expressed in terms of 
continuum variables, and the scaling factor in front tells us that r] = 1 if we compare with 
eq. ( p.45| ). Likewise eq. ( |2.41| ) and m(/i. A) = A/x shows that u = 1, i.e. (Ih = 1- We 



have a new class of smooth random walks (du = ^)- We have already shown that 7 = 1 
and the exponents 7, v and rj is our first example of a set of non-trivial exponents. They 
satisfy Fischer's scaling relation. 

The ordinary random walk has a stochastic interpretation as a Brownian motion of a 
particle, i.e. each step is performed according to some probability distribution P, but 
is independent of the former steps. With the extrinsic curvature term added the step 
also depends on the direction of the former step. For a finite coupling constant this does 
not change the universality class of the random walk, but as we take A ^ 00 we enter 
a new class of random processes characterized by different critical exponents. The in- 
terpretation of these is that the velocity, rather than the position of the particle itself, is 
changed stochastically according to some probability distribution P. The path in an ordi- 
nary random walk will be continuous, but, with probability one, nowhere differentiable. 
If the velocity is stochastic the typical path will be differentiable and the first derivative 
continuous but nowhere differentiable. There is a number of stochastic processes which 
have this feature. The well known Ornstein-Uhlenbeck process [ 3] is one of them. It can 
be shown that the propagation of particles in such processes are described by propagators 
of the type (2.74), but we have to refer to the original articles for details [ 2]. 

The progagator (2.74) is also related to the propagation of a spinning particle with 
infinitely many components. This will be discussed at the end of the next subsection. 
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2.5. Fermionic random walks 

For the ordinary bosonic propagator we found the following representation: 



The fermionic propagator in D dimensions is obtained by the change 
1 1 



(2.76) 



where the D 7 matrices satisfy: 

{7^7n = 25'^^ (7^)^=7^ (2.77) 

The lowest dimensional representation of this so-called Clifford algebra is by matrices of 
dimension n = [D/2]. In two dimensions one can take the Pauli matrices ai and o"2. In 
three dimensions one can use all three Pauli matrices. The 7-matrices can be used to 
construct a representation of SO{D): The matrices 

s''' = l^'' = l[r,ll (2.78) 

satisfy the commutation relations for the generators of SO{D) and rotations will be gen- 
erated by 

Kicu) = e-^^-^V^ _ ^2.79) 

These rotations act on spinors, i.e. vectors in C", n = [D/2]. If 

u;>^^ = 6n^'^{e„ e,), n^-^ = ^'f 7 f ^% (2.80) 

2sm6' 

where 6 is the angle between ei and 62, the rotation (|2.79|) is a rotation with angle 6 in 
the plane spanned by ei and 62- An important aspect of the spinor representation of the 
rotation group is that a 27r-rotation gives -1. This well known fact will be of outmost 
importance for the fermionic random walk. 

Let x^{^) be a curve in R^ and let t^^{C,) be the normalized tangent vector, defined by 
( |2.47| )-( p.48| ). t^ and t^ span the osculating plane of the curve and we can write: 

u"^ = -{tH'' - tH") = nn^r (2.81) 

Recall that the curvature n is related to the angle of contingency by ( p.50| ). This implies 
that we can write 

^^"({1= den^"" (2.82) 

where n^^ is the antisymmetric tensor which defines the osculating plane and d9 is the 
angle between tangent vectors t{l) and t{l + d/). If we consider a discretized random walk 
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and as usual denote the two successive unit vectors as Cj and Cj+i and the angle between 
them as ^(cj+i, e,) the discrete analogy of ( |2.82| ) is 

u;'^"(ei+i, e,) = ^(e^+i, e,)n'^"(ei+i, e,) (2.83) 

Let if) G C" be a spinor. Let us imagine it "propagates" along the given path in such 
a way that it is always rotated according to the orientation of the curve, i.e. for each 
discrete step we perform a rotation given by the rotation matrix 

K{ei+i, Ci) = e--'"'-^''(^»+i'^0 . (2.84) 

The total rotation during a travel along the path P„ will be given by 

K{Pn) = K{en, e„_i) ■ ■ ■ K(e2, ei). (2.85) 

A formal continuum version of K{Pn) for a smooth path is 

K(P)=Pe-i/p^'^^^-^^« (2.86) 

where V denotes the path ordered integral. This factor appears in the famous Strominger- 
Polyakov representation of the fermionic propagator[ 4]: 

G{x,y)= fvP{x,y) e^^^^^-'V)^^ K{P{x,y)). (2.87) 



The definition is identical to the one for the bosonic particle except for the matrix K, 
which rotates a spinor "along the curve". The factor K{P) is rather ill defined and it has 
been difficult to use this expression. However, from the above definitions it is clear how 
to write down a well defined discretized version of ( ^.87[ ) [ 5]: 

/n n 

n de, n ^(e.+i, eO 5{Y. e, - x), (2.88) 

» i=l 4=0 



where Cq = Cj and e/ = e„+i. This expression is very similar to ( p.55|) for the bosonic 



particle and we can directly use the machinery developed in the last subsection. The 
matrix i^(e2,ei) can be viewed as the kernel for an operator K which acts on wave 
functions belonging to 7i = L{S^~^) x C"^, i.e. spinors on S^~^. If we normalize K we 
have: 

{k^m = J dei K{e,e,) ^(ei), J dei K{e,e,) = 1. (2.89) 

In the same way we can view G^{ef,ei,x) as the kernel of an operator G{x) on H. We 
get rid of the (5-function in eq. ( p.88| ) by Fourier transformation and express the Fourier 
transformed G{q) as 

Guiq) = ^ ^ (2.90) 

in the same way as for the scalar particle (see eq. ( p.66|) ). 
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Let us now define the following scaling limit: 

fi = mphO, q = pphtt, (2.91) 

i.e. the propagator becomes: 

1 



G,iq) = - 
a 



l-K 



— + mph + ipph ■ e + 0{a) 



(2.92) 



As usual p ■ e should be viewed as an multiphcation operator. As for the scalar particle it 
can be shown that 1 is the largest eigenvalue and that there is a gap to the next lowest 
eigenvalue (there is no extrinsic curvature which allows us to tune the gap to zero). It 
follows that only eigenvectors corresponding to the eigenvalue 1 will propagate for long 
paths and contribute to matrix elements of G^{q) for a —>■ 0. For the scalar particle we 
only got a non-trivial result from the scaling ansatz ( ^.91| ) in the limit where the extrinsic 
curvature term reduced the mass gap to zero. Here it is different because the eigenvalue 1 
is degenerate. It is 2?T,-times degenerate and the operator G^{q) acting on the Hilbert space 
Ti reduces in the scaling limit to 2n x 2n matrix acting on the eigenspace V corresponding 
to the eigenvalue 1. 

Theorem: The constant vectors and the column vectors of the matrix '-^^^e^ span the 2n 
dimensional eigenspace V oi K corresponding to eigenvalue 1. 

That constant vectors on C" are eigenfunctions of K of eigenvalue 1 is clear. By definition 
of the rotation we have: 

K{e2, ei)7 ■ eiK-\e2, hci) =7-62. (2.93) 

This shows that the columns of the matrix 7 ■ e are eigenfunctions of eigenvalue 1 since: 

K7-e = /de,K(e,e0 7-e, = 7-e/de,K(e,eO=7-e. 

We leave it as an interesting exercise to show that there are no other eigenvectors corre- 
sponding to the eigenvalue 1 (see [5]). 

The projection operators P^ = 4^(1 ± 7 ■ e) commutes with K and split H into "±- 
chirality spaces": 

Theorem: On V± the multiplication operators e^ can be replaced by ±7^-matrices. 

The proof is simply by calculation: Let ip{e) and 0(e) be two vectors in 7i. The scalar 
product on TC (and V) is 

(0|^)=y'de0*(e)„^(e)„ 
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A^ = M = 2Tt 

Fig. 5. The cancellation between two paths. For the first we get a phase factor e*"^^/^ = 1, while the 
second gives a phase factor ef-^^l"^ = — 1. 

where a is a spinor index. Let us now apply this to n vectors f^^-'(e) = P^p.^^)- These are 
orthonormal vectors on V^_|_. It follows that 

We can finally state that the scaling limit (|2.92| ) is given by 



^.(g)-- p%^^^'^-" 1 ■ (2.94) 

We recognize two copies of the Dirac operator corresponding to Y±. This doubling is 
needed by the Nielsen-Ninomiya theorem. 

On the way we have determined the critical exponents for the fermionic random walk. 
From eq. ( p.92| ) and G^ici = 0) = / &xG{x) we get 7 = 1. From m(/i) = /i we get v = \ 
and from Gfj_{q) ~ 1/a we get rj = 1. The fermionic random walk has the same critical 
exponents as the smooth random walk. This implies that dn = ^ and effectively we have 
smooth paths. However, we have no extrinsic curvature term to produce the smoothness. 
Rather, it comes about because of cancellations between a large number of bosonic paths. 
To show this let us consider two dimensions (but the mechanism is the same in higher 
dimensions). From (^.82|) it follows that we can write: 



uj^^s^"" = — (T3 (2.95) 

di 

and the eigenvalues of cr-^ is precisely the split of H in TC±. On one of these spaces K{P) 
( |2.86|) becomes a phase factor and ( p.87|) allows a scalar interpretation: 



G{x) = fvP{x) e-™/<^'+5/de^ (2.96) 

where 6{l) is the angle of the tangent relative to a fixed direction in the plane. It is 
now essential that we have spinors since the —1 resulting from a 2it rotation leads to 
a cancellation between intersecting and non-intersecting paths as shown in fig. 5. Only 
smooth paths survive this cancellation. 

Even if the mechanism for producing the smooth paths seems very different for the 
bosonic and the fermionic particles there is a connection. As is rather clear from the 
derivation of the bosonic propagator (2.74) the tangent vectors are effectively performing 
a random walk on S^~^. This is just what one would expect from a classical spin and 
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this is made clear from the expression for the operator K\ in the scahng hmit: 



K. 



Q-\\iL a 



(2.97) 



Since K\ is the transfer matrix for proper time evolution along the path and since the 
Hausdorff dimension is one, L^ can be viewed as the continuum Hamiltonian, and we know 
that it is precisely the Hamiltonian for a classical spinning particle. This interpretation 
is substantiated by canonical quantization of the particle with extrinsic curvature. The 
Lagrangian contains higher derivatives and this implies that x/\x\ will be a new additional 
canonical coordinate which will serve as a the spin. We have to refer to [ 6] for further 
discussion. Here it is worth emphasizing that this result is quite natural in the discretized 
approach. In addition it is seen that the wave function of such a classical spinning particle 
is an infinite component object as is clear from eq. (2.74). It is closely related to the 
infinite component spinor introduced by Majorana in an attempt to avoid the negative 
eigenstates of the Dirac operator. From this point of view it is possible to consider the 
Dirac propagator in a more a more general context. If the matrix A'(ei,e2) had been 
generalized to contain an extrinsic curvature term, one would have obtained a scaling 
limit similar to (2.74), the only difference being that 



Aph-^ 



Aph(-^ 



(2.98) 



where s^ is the spin operator while J^ refers to the total angular momentum operator. 
For Aph — > oo it will project to the finite dimensional subspace which we obtained above 
for the ordinary Dirac operator but in a limiting process where the Hausdorff dimension 
always is one. 

It is natural to ask the following question: Although we have found a nice discretized 
version of ( p.87| ) and have shown that it leads to the Dirac propagator in the scaling 
limit, ( p. 871 ) itself is not a path integral of an action. Is there a path integral which leads 
to ( p.87| )? The answer is yes! The supersymmetric generalization of the bosonic action 
( |2.18| ). The supersymmetric partner of the field x^ is a Grassmann variable ip^ while the 

supersymmetric partner of e(^) = J g{i) is a "gravitino" field x(0 ^^"^ ^^ have 



Sb ■ 

Op ■ 



de 



-X + e 
e 



e e 2 \dO 



(2.99) 



It is possible to show that the path integral over the Grassmann variable x ^i^nd ip results 
in the factor K{P) in ( p.86|) . Details can be found in [4]. 



3. Random surfaces and strings 

3.1. Definition of the model 

The theory of random paths described the relativistic particle. We expect the theory of 
random surfaces to describe the relativistic string. The strings sweep out a surface while 
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they propagate. The path integral is a sum over all such surfaces with a weight given by 
the classical action. It is our goal to define this sum and analyze it in detail. 

As for the relativistic particle we have two actions which are equivalent at the classical 
level. The first action is geometrical, only determined by th£ area spanned between 
the starting position and the end position of the closed stringa. Let M(/j) denote the 
manifold with boundaries k, i = 1, ...,n, F{li) a corresponding surface in R^ and a;'^(^) 
the coordinates of F in R^. The action is 



S[Fik)] = f^ dAiF) (3.i: 



2/o,,\2 /o,, o,,\2 






An alternative description is obtained by introducing an internal metric Qab, a,b = 1,2 
on M{li) and use the following action: 



S[g,x] = ^f d^^v^ 
a' Jm(l) 



a JM{k) 



9 g^-g^+^ 



(3.2) 



The classical equations for the actions (3.1) and ( |3.2| ) agree, but it is not at all obvious 
that the quantum theories are identical. In the case of strings we have many "natural" 
objects, in contrast to the situation for the free particle where one only has the two-point 
function. A surface can join an arbitrary number of strings and one is led to consider the 
n-loop amplitude G(Zi, ..., /„) between the n loops or strings. The formal path integral 
expression for the n-loop is written as: 

G{h,...,ln) = JvF{h) e-^[^('')], (3.3) 

where the integration is over physical distinct surfaces F in R^ , or 

where the integration is over all equivalence classes of metrics on M{li) and all embeddings 

x'^iO. 

In order to define (|3.3| ) and ( |3.4| ) we have to introduce a reparametization invariant 
cut-off. We follow the procedure outlined for the relativistic particle. In the case (|3.3| ) 
it amounts to use as a building block a smallest triangle in R^ and glue together these 
triangles in all possible ways to surfaces with the given boundary conditions. Alterna- 
tively one could consider the hyper-cubic lattice surfaces where the surfaces are made of 
plaquettes [ 10, 11]. Many of the results we will obtain in the following are valid (and 
easier to prove) for these models than for the model defined by eq. ( p.4|) . However, eq. 
( |3.4|) relates closer to quantum gravity, since it is just two-dimensional gravity coupled to 
D free scalar fields x'^, and for this reason it is convenient to use here eq. (3^) rather than 



eq. (|3.3| ) as we are going to consider quantum gravity in some detail [ 7, 8, 9] (again 



For simplicity we consider here only the theory of closed strings. 
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there is a large number of articles describing this approach [ 13, 12, 14, 15], just to 
mention some of the articles which concentrated on computer simulations of the model). 
In the case of the random walk the reparametrization invariant cut-off introduced was 
related to the shortest distance. In the two dimensional case it is natural to combine 
shortest distance and smallest area in a single cut-off. The fundamental building block 
will in this way be an equilateral triangle with edge length a. At this point we encounter 
a new problem compared with the one-dimensional situation: The gluing is in no way 
unique. In the process of gluing together triangles to form a two-dimensional manifold 
the order of a given vertex (i.e. the number of triangles to which the vertex belongs) is 



almost arbitrary. In the case where we use eq. (|3.3|) and glue together triangles directly in 
R^ the answer is clear: We should include all distinct different ways of gluing compatible 
with the boundary conditions since we will get different surfaces in R^ . But also for the 
model ( p.4|) where the triangles are defined with respect to the internal metric the freedom 
of gluing will go hand in hand with the fact that a closed surface, apart from the total area 
(the equivalence of the total length of the path), also has a new local invariant: R{C) • the 
Gaussian curvature. R{^) cannot be changed by a reparametrization of the surface. In 
the following it will be argued that the sum over triangulations in a precise way captures 
this new degree of freedom. 

Let us make a short digression and discuss curvature. The curvature (or Riemann) 
tensor R^bc can be defined in terms of the metric. It describes the deviation from fiat 
space. This is manifest in the formula for parallel transport of a vector S'^ around an 
infinitesimal closed curve: 

^Sa = \R''abcSdjx'dx\ (3.5) 

The once contracted tensor Rab is called the Ricci tensor and the scalar obtained by 
contracting the Ricci tensor is called the scalar curvature R: 

p pc p pa 

-fl-afe — J^ ach^ -'^ ~ -'^ a- 

In two dimensions there is only one independent component: 

Rabcd = -{dacgbd — daddbajR, Rab = -QabR, R = '^K. (3.6) 

where K denotes the Gaussian curvature on the surface. K has the simple geometrical 
interpretation as being the product of the principal curvatures associated the normal 
planes intersecting the surface which is assumed to be embedded in R^ in order for this 
interpretation to make sense. But R itself is of course independent of this embedding. 
Eq. ( |3.5| ) simplifies: 



ASa = -RS'' dAab, d9=-RdA = K dA, (3.7) 

where dAat = \§{xadxb — Xbdxa) is the area tensor and d9 the infinitesimal angle by 



^ By different ways of gluing we have in mind that the resuhing (abstract) triangulations are different 
in the way defined below. 
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Fig. 6. A geodesic triangle on a smooth surface. 

which S"" has been rotated during the parallel transport. There is a nice integrated version 
of this relation, known as (one of the versions of) the Gauss-Bonnet theorem: Let T be a 
geodesic triangle on the surface, i.e. a triangle where the sides are geodesic curves, with 
angles /3i,/52 and Ps. The sum of the angles is no longer vr but the deviation from the 
Euclidean value is given by the integral of the Gaussian curvature over the interior of the 
triangle: 



K dA = pi + p2 + p3 - T^ = £t, 



{3.t 



where Et is called the excess angle of the triangle. The situation is illustrated in fig. 6. 
Parallel transport of a vector along the boundary of the triangle will rotate the vector by 
the angle Et- 

Due to Regge [ 16] we can define curvature and parallel transport in a natural way on 
piecewise linear surfaces. The curvature cannot be located in the interior of the triangles 
since we view the interior as flat. Since the curvature is defined as an intrinsic geometric 
quantity it is clearly bending invarianTO. Since we can bend the surface around an edge 
without changing anything we cannot use the edges either. In this way we are lead to 
locate the curvature of the piecewise linear surfaces at the vertices. To each vertex v we 
associate a deficit angle Ey by: 



Ey = 27r-^aj,(t), 



(3.9) 



tBv 



where the summation is over the f-angles of the triangles to which v belongs. This is 
illustrated in fig. 7 where a geodesic triangle with the vertex v of the piecewise linear 
manifold is shown. It follows from fig. 7 that Et = Ey on the piecewise linear surface, i.e. 
we can write: 



[xdA^Ey, I K dAr-^Y^^v 



(3.10) 



^ Historically the Gaussian curvature was defined in terms of the principal curvatures ki and K2 , which 
explicitly depended on the embedding. It came as a surprise that K = kiK2 is independent of the 



embedding (bending invariant). 
Gauss. 



This is a "Theorema egregium" 



■'most excellent theorem" , wrote 
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ai + q;2 + «3 = 27r — e^ 



Fig. 7. A geodesic triangle on a piecewise linear surface. The vertex v is an interior point of the triangle, 
(fig. 7a). Fig. 7b shows the excess angle after the piecewise linear neighborhood has been cut open along 
a link and unfolded in the plane. 

The Ihs of the equations are intended to be vahd for smooth surfaces, the rhs for piecewise 
hnear surfaces. In the last equation we have generahzed the formula to any region where 
the boundary is a piecewise geodesic curve and the summation is intended to be over 
interior vertices in the triangulation. In particular the formula will be valid for closed 
surfaces. 

Let us split e„ in a curvature term and an area term by assigning the following area to 
the vertex v. 



A. 



t3v 



K„ 



6y 



(3.ii: 



i.e. the area of each triangle is distributed equally among its three vertices. With these 
definitions we have finally: 



I dA ~ ^A, JKdA ^Y.KvA, 



(3.12) 



Regge originally intended to use these formulas by constructing a sensible sequence of 
piecewise linear approximations to a given smooth surface such that 



t1 J 



(3.13) 



Y^f^K^A^JfiOKiOdAiO. 



(3.14) 



On the same grid one can discretize the covariant action of the free scalar fields in 
a natural way. For a given triangulation we denote the vertices by indices i,j,k,.... 
The triangulation is characterized by its coincidence matrix, which specifies the neighbor 
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vertices to a given vertex (see below for definition), and by the length kj of the links 
between vertices i and j. A natural coordinate system is introduced on the piecewise 
linear surface by assigning coordinates yi to the vertices i such that kj = \yi — yj\. The 
l/j's live in some ambient space /?". The interior of the the triangle (ijk) is parametrized 
by barycentric coordinates and for a field defined at the vertices by 0j we use in the 
same way the linear extension to the interior of the triangle: 

y = ey^ + eyj + {l-e-e)yk (3.15) 

Hy) = ev. + ev, + {i-e- e)<pk (3.16) 

The metric will be defined by 



QyO QyC 



9ab = ^;^^;^, a, 6 =1,2 (3.17) 



and it is straightforward to calculate Qah and g da 



/ 72 l(n 1/2 _ /2 N\ 

Hk 2\''ik ^ ''jk Hj) ('^^R\ 

Uab — I 1 /,2 I 72 _ /2 \ /2 I ' 1,0.10; 

\ 2 \ ifc """ jfc ijJ jk / 

g^'da<pd,<P=- fe(0^ - 0,)(0. - 0,) + 

g L 

^ifc(0i -0i)(0fc-0i) + 4(0i -0i)(0fc-0i) , (3.19) 



\2 1 fi2 ;2 I 72 r2 , r2 r2 1/74 ,74 ,74 



where 

5' = (2A.yfc)^ = - f Z^^jfc + ^ij^kj + ^iJij - 2^^tj + ^tk + ^jfc) j • (3.20) 

After this digression let us return to the problem of regularizing the integration of 
Riemannian structures of two dimensional manifolds. The reparametrization invariant 
regularization suggested above consists of constructing all piecewise linear manifolds ob- 
tainable by gluing together equilateral triangles of side- length a and assigning to these the 
metric structure offered to us by Regge calculus. In the case of equilateral triangles the 
formulas above simplify a lot: Let n^ denote the order of vertex f in a given triangulation, 
i.e. the number of triangles which contain vertex v. From eq. ( p.ll|) we get: 



1 v3 2 A 27r 

A^ = -ny—a, R^Ay = —{6-ny). (3.21) 

Let us for simplicity of notation absorb a/3/4 in a^ and put the resulting a = 1 in the 
following. The formulas for curvature and for the Gaussian action then read: 

A, = -n,, R^ = 2n -, (3.22) 

3 fly 
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J2A, = Nt, 5]A,i?, = 47rx, x = Nt-Nl + Nv. (3.23) 

V V 

Here Nt, Nl and Ny denote the number of triangles, links and vertices in the triangu- 
lation, and x denotes the Euler characteristic of the rnanifold. (|2.33| ) is the discretized 



version of the Gauss-Bonnet theorem mentioned aboveH. Finally the action of the free 
Gaussian field given by eq. (3.19) in the case of equilateral triangles becomes 

d'e V99''%<pd,(P -^ 5:(0, - <Pjf (3.24) 

m 

From ( |3.22|) one observes that different triangulations, i.e. triangulations which cannot 
be mapped onto each other by a simple relabeling of the vertices, lead to different local 
curvature assignments and consequently inequivalent metric structures. In two dimensions 
a closed manifold is characterized entirely by its Euler number. Given a manifold we 
want to integrate over equivalence classes of metrics. Since all the different triangulations 
we can construct by gluing the equilateral triangles together correspond to inequivalent 
metrics it is clear that one should sum over all such triangulations. By this prescription 
one approximates a continuous integration over metrics by the summation over a grid 
of points in the space of inequivalent metrics. The conjecture is that this grid becomes 
uniformly dense when the number of triangles Nt of the triangulations goes to infinity. 
We shall later on verify this conjecture. 

Although we have used Regge's prescription for assigning curvature, the philosophy 
outlined is very different from the one which motivated Regge. In the classical Regge 
calculus the objective was to approximate a given smooth surface by a piecewise linear 
manifold. A fixed triangulation was chosen and the link length treated as the dynamical 
variable which should be adjusted to get the best approximation to the given manifoldQ. 
In particular, different link assignment will not necessarily result in a different metric 
assignment, as is clear by considering triangulations of the plane. Clearly there is a lot 
of room for moving the vertices around (and thereby changing the link length) without 
changing the metric at all. Integration over link length is not an integration over equiv- 
alence classes of metrics, but involves a highly non-trivial Jacobian. Here we are not 
interested in approximating specific manifolds, but in using different triangulations to 
label different equivalence classes of metrics. 

The regularized definition of the multi-loop Green functions (|3.4| ) reads: 

G,{h,...,ln)= E e"""""/ n dx. e-S(«.)(---^)\ (3.25) 

TeT(h,...M i£T/{h,...,l„} 

In this formula T denotes an abstract triangulation, defined by its vertices i and a table 



^ It is easy to prove eq. ( 2.33 ) from Euler 's formula which states that x — Np ~ N^ + Ny for any polygon 
net covering a surface of topology characterized by x- If oii^ use Euler's formula on a triangulation 
{2Nt = 3Nl) and combine it with ( ^3.22 ) one arrive at the formula in ( ^.25 ). 



^ The given smooth manifold in the Regge approach was the one given by solution to Einsteins equations 
with suitable boundary conditions and it would be an extremum of the Einstein-Hilbcrt action. An 



important feature in this context is the convergence (3.14) (in higher dimensions where the Einstein- 
Hilbert action is not a topological invariance), which ensured a good approximation to the action if the 
triangulations were chosen well in accordance with the geometry of the problem. 
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which tell us the neighbor vertices to a given i. This information is encoded in the 
coincidence matrix Ct which is an Ny x A'^v' matrix where the {ij) entry is —1 ii i and j 
are neighbors, if they are not neighbors and n^,, the order of the vertex, on the diagonal. 
The links are defined to be pairs {ij) of neighbor vertices and the triangles are defined as 
triples of neighbors {ijk) such that each link (ij), which is not a boundary link, belongs to 
precisely two triangles (ijki) and (ij7c2)- The summation is over (izjferent triangulations. 
Two triangulations are considered as identical^ if there is a map between the vertices 
compatible with the assignment of links and triangles. U has a two-fold meaning as a 
fixed polygon loop in target space R^ and an abstract boundary in the triangulation 
T and T(/i, ...,Z„) denotes a suitable class of triangulations with the given boundaries. 
Usually we have in mind all triangulations of a given topology x- However, occasionally it 
is convenient to enlarge the class of simplexes considered, such that they strictly speaking 
do not form a combinatorial manifold. Local "irregularities" of this kind should not be 
important, since they are related to short distance effects which should not play any 
role in the continuum. More serious is the restriction on topology. Formula ( |3.25| ) is 
very tantalizing in the sense that it has no reference to topology. Is it possible that eq. 
( p.25|) provides a non-perturbative definition of the summation over topologies? It has 
always been an annoying aspect of the continuum formula ( |3.4| ) that we only know how 
to interpret it for a given manifold, i.e. a given topology in the two-dimensional case. A 
summation over different topologies has to be performed by hand. It turns out that eq. 
(\3.2^ can be used to study the summation over topologies, but not directly as it stands. A 
special limit, the double scaling limit has to be taken. This will be discussed later. At the 
moment we will always restrict the class of triangulations T to mean triangulations with 
a fixed topology, usually the simplest, the spherical topology. In the following we will 
also use the notation spherical topology for surfaces with boundaries where we recover 
the sphere after closing the boundary. 

It is often useful to consider a number of special cases of ( |3.25|) . If there is no loops at 
all we talk about the partition function for closed surfaces: 

Z{ii) = ^ ^e-^^^ / n dx, e-^«)("'-"^)'. (3.26) 



c 

Ter "^T -^ ier/{io} 

In this formula is included an additional symmetry factor St for the triangulation. It 
is similar to the additional factor which appears in vacuum Feynman diagrams, and it 
reflects the additional symmetry which can be present for surfaces without a marked 
boundary: A permutation of the vertices might leave unchanged the links and triangles 
and in this way not change the surface. 5*^ is equal to the order of the automorphism 
group of the graph T. A vertex zq is excluded from the integration in order to kill the 
mode associated with translational invariance. The Gaussian integration is independent 
of this choice and alternatively one could have chosen to fix the center of mass. 

Another limiting case arises if we contract the loops to points, i.e. marked vertices. 
Strictly speaking this cannot be done in a continuous way on the triangulations. We 



^ Often the word equivalent is used for such two triangulations. However, we are going to reserve this 
notation to triangulations which have a common subdivision. 
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Sn 



Fig. 8. To each surface SN,i with one marked point Xi corresponds N —1 surfaces SN,i,j with two marked 
points since we can put yi at the A'^i other vertices. This is true up to symmetry factors. 

denote the n-point function G(xi, . . . , x„). It will be given by 



G,{xu...,xn)= E ^— / n dx, e-^W)^^-^^^ (3.27) 

TeT{iu...,i„) "^T •' ier/{n,...,i„} 



where the symmetry factor can be different from zero for the 1- and 2-point function. Note 
that the 1-point function (which by translational invariance is independent of the target 
space point Xi) is (essentially) equal to (minus) the derivative of the partition function 
with respect to /i, while the integral of the 2-point function is (essentially) equal to the 
double derivative of the partition function: 

G^(xi) ~ -Z'ifi), J dx G^{x,y) ~ Z"{fi). (3.28) 

The equations follow from the observation that differentiation of Z{fi) multiplies each tri- 
angulation with a factor Nt coming from e~^^'^ . The 1-point function is the summation 
over marked triangulations, but there are Ny of these for each triangulation without a 
marked vertex, up to symmetry factors which play no role for the generic large triangu- 
lation. Relations like ( p.28| ) will be valid in the limit where triangulations with large Nt 
dominate. This is illustrated in fig. 8. As we shall see shortly the model will have a critical 
point He- For /x above the critical point all G^ will be analytic functions of /i, but at the 
critical point they will contain non-analytic parts. These are the universal parts which 
have our interest. They are determined by the large N^ part of the triangulations since 
finite A^t's only produce analytic contributions. The word "essentially" above and the 
symbol "~" in ( p. 2^ ) refers to this non- analytic part determined by the large Nts. The 
line of argument can be extended to the n-point function. Let the generalized susceptibil- 
ity x^"'^ (/x) be defined as the integral over n — 1 of the arguments, x^"''' (/^) is "essentially" 
equal to the derivative of x^'^~^'^ {lA since the derivative brings down a factor Nt in the def- 
inition (|3.27|) while x^'^'' (/^) contains an additional marked point which produces a factor 
Ny (~ Nt/2) in the counting of surfaces. It is summarized in the formulas: 



^W 



ifx) = I dx,--- dx„_i G^(xi, . . . , a:„) ~ (-l)"^Z(/i). (3.29) 
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3.2. Physical observable 



Eqs. ( |3.25|) - (|3.28| ) defined the regularized loop-functions. The following theorem ensures 



the existence of a critical point /ic like for the random walk: 

Theorem: For surfaces of spherical topology exist a critical point fic such that G^{{li}) 
is convergent for fi > fie and divergent for ft < fie- /^c is independent of the boundary- 
loops {k}. 

Conjecture: The theorem is true for surfaces of any fixed topology x and /Xc is the same 
for all X- 

The theorem will not be proven here, but it is not difficult to show that Z{fi) is well 
defined for // sufficiently large. Two steps are needed. First we have to bound the 
Gaussian integral for a given triangulation. In fact it is not difficult to show that there 
exists constants q and c„ such that 



e" 



c,Nt ^ f Yl dxi e'^c^)^'"'"''^' < e'-^^ . (3.30) 

ieT/{io} 



Secondly we have to use that the number of triangulations of a fixed topology is expo- 
nentially bounded. Let us denote the number of triangulations with topology x? which 
can be constructed from N triangles as A/'( A^, x) '■ Constants di and du exist such that 

^d.Nr < X^N^^ ^) < e'^"^^ . (3.31) 

Combining eqs. ( p.31|) and ( |3.3CI|) we conclude that for any x'- 



ci + di < fic< Cu + du- (3.32) 

For the rest of the theorem we refer to the original article [ 7]. The conjecture that fi^ is 
independent of x is almost certainly true. It has been proven for sums over triangulations 
coupled to matter with central charge c < 1. 

It is worth to notice that the exponential bound on the number of triangulations plays 
an important role for the existence of the critical point. If we try to define the summation 
over all topologies directly from eq. ( |3.25| ) we will fail due to the entropy of triangulations. 



The number of triangulations of A^ triangles A/'(A^), with no restriction on topology grows 
faster than factorially and from eq. (|3.3CI|) it follows that the sum ( p.25| ) is ill defined for 
any choice of /i! Later we will discuss some attempts to make sense of the sum ( |3.25| ) 
after all. 

After the existence of a critical point is established it is of interest to study the critical 
behavior of the n-loop and n-point functions when we approach the critical point. 

The most important quantity in this context is the mass gap since it determines the 
possible scaling. The mass gap can be defined by the exponential decay of the two-loop 
function as the distance d between the loops goes to infinity: 

m(/i) = - lim ^ALIA^ (3.33) 

d-*oo d 
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Fig. 9. The sum over surfaces from I to l^ is larger than the restricted sum there is a "bottle neck" at 
distance di. 

where la denotes the loop / displaced a distance d. 

Theorem: The two-loop function falls off exponentially with the distance between the 
two loops. 

The argument is the same as for the random walk: since G^(/, ^di+dj) is the unnormalized 
probability density for a propagation of a string from / to ldi+d2 it will, multiplied by 
an appropriate normalization factor which converts it into a probability, be larger than 
the corresponding product of G^{l,ldi) and G^ild^^ld+d-i)) since this product (correctly 
normalized) imposes the constrain that the surfaces from / to ldi+d2 should pass through 
Id^- This is illustrated in fig. 9. I.e. the (correctly normalized) two-loop function is sub- 
additive and since the normalization factor is independent of d the limit ( |3.33| ) exists. It 
follows from the definition that m{^) is an increasing function of /x. Unfortunately there 
exists no proof that m{ii) goes to zero for /x -^ /ic- We will assume this is the case. We can 
then define the same critical exponents and fractal dimensions as for the random walk, 
but let us first introduce a new critical exponent p corresponding to the string tension. 
We define the string tension cr(yu) by the exponential decay of the 1-loop function G^{Ia) 
for a large planar loop I a which encloses an area A ~ /^: 

G^(/a) -^^e-'^^'^)^. (3.34) 

The proof that the (^^(/a) falls off exponentially is again based on sub-additivity: Let A 
denote both the planar region enclosed by Ia and its area. The number of surfaces which 
has Ia as boundary is larger than the sub-ensemble where we divide A in two pieces with 
area Ai and A2 by introducing a new boundary in the middle of A which we force the 
surfaces to respect. We conclude that G^{Ia) > G^{Iai)G^{Ia2) and this implies ( p.34| ) 

Why is cr(/i) called the string tension? We can view G^j_{Ia) as the partition function 
Z{A) for an ensemble of surfaces which are allowed to fluctuate, but where the boundary, 
i.e. the frame, is fixed. The Gibbs free energy of the system is F{A) = — log Z{A). The 
string tension is defined by change in free energy per unit area if we change the area A 

AF = aAA. (3.35) 

From p. 34]) it will precisely be our a{fi) for large A. 



Let us now introduce the scaling parameters: 
(1): The critical mass exponent u is defined by the assumed scaling of m{fi) to zero. This 
allows us to introduce the physical mass and a length scale which goes to zero at the 
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critical point. 

■m{fi) ~ (yU — ficY, m{fi) = mpha(/i). (3.36) 

(2): The short distance behavior of the 2-point function is characterized by the anomalous 
scaling exponent 1]: 

G^(x,y)~|x-y|2-'^-^ l<|x-i/|<l/m(/i), (3.37) 

(3) The susceptibility x(/i) and the susceptibility exponent 7^ is defined by: 

x(/i) = dx G^{x,y) ~ — . (3.38) 

(4) The string tension cxlfi) is defined by the exponential decay of the 1-loop function and 
we assume 

a(/i) ~ (/i - /ij". (3.39) 

From dimensional arguments we expect p = 2v. However, we will prove that o"(yu) does 
not scale to zero for the simplest model. 

(5) The extrinsic Hausdorff dimension du of the ensemble of surfaces is defined by 

{Area)^r^{NT),r^r^", (3.40) 

where the average area is over an ensemble of surfaces with two marked points a distance 
r apart in target space R^ . 

Let the topology of the surface be spherical except for possible boundaries. If the mass 
scales to zero it follows, by arguments identical to the ones presented for the random walk, 
that 

7, = z/(2-77), dH = l/u, (3.41) 

i.e. Fischer's scaling relation and the relation between the mass exponent and the Haus- 
dorff dimension. For most systems one has that < r^ < 2. In ordinary statistical systems 
?7 = 2 at the infinite temperature limit, while r/ = is the Gaussian approximation. For 
such systems it is clear that z/ > implies 7 > 0. In the following we will often encounter 
systems where 7^ < 0. It is important to realize that such systems exist. As a simple 
model one can take the closed random walk with two marked points Xi and Uj kept fixed 
in R^ . It has sovne analogy with the 2-point function for random surfaces^. Since it is 
clearly the product of two ordinary random walks between x and y the mass, which is 
determined by the exponential decay at large distances, will be twice that of the ordinary 
random walk, i.e. it has the same critical exponent v. The short distance behavior will 
be 

G,(x)~^^, i.e. r^ = d-2. (3.42) 



® It does not describe the propagation of the desired geometrical object, the particle, between x and y, 
while the two-point function for the closed string indeed describes the propagation of the string between 
X and y. 
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For d > A it follows that 77 > 2. It is easy to show directly that 7s = 2 — d/2, i.e. Fischer's 
scaling relation is valid even if 7 < 0. 

There exists a theorem which indicates that it might be difficult to obtain a positive 7^ 
for random surfaces. 



Theorem: 7s < 1/2 for the random surfaces model ( p. 251 ) with spherical topology. 



Conjecture: 7s > implies 7s = 1/2 for the random surface model ( |3.25| ) with spherical 
topology. 

The theorem and the conjecture are known to be true for the hyper-cubic random surface 
model [11]. Rather than giving the rigorous arguments let me present the underlying 
geometrical reason for the theorem and the conjecture. But first a remark about the 
technical point which should be dealt with if the arguments should be made exact. In 
order to apply the cutting and sewing arguments below one would have to introduce 
somewhat more complicated objects than the n-point functions we have been considering 
until now: the correct objects to consider are n-loop functions where each boundary loop 
consists of, say, three links. Fix the center of mass of the boundary vertices of loop i to 
be Xi, the corresponding point in the n-point function, but integrate over the positions 
of the boundary vertices compatible with these constraints. The corresponding n-loop 
function will be a function of Xj, precisely as the ?7,-point function and we expect that 
they coincide in the scaling limit where the "bare" distances Xj all scale to infinity while 
the "physical" distances Xph = x a{fi) stay fixed as yU — > fic, i-e. a(/i) —>■ 0. The reason 
is that any contribution where the distances between the boundary vertices are larger 
than 1 will be exponentially suppressed by the Gaussian action and in the scaling limit 
distances of order 1 mean physical distances of order a{fi), i.e. of the order of the "lattice 
spacing" . Had we considered surfaces made of plaquettes living on the hyper-cubic lattice 
or surfaces in R^ build directly from the gluing of equilateral triangles these problems 
would be absent and the arguments to be presented would be exact in the sense that one 
would not have to integrate over positions of the boundary vertices. In the following we 
will ignore these complications for the sake of argument. 

Consider the generalized susceptibility x''"^/^) ^s defined by eq. ( p.29|) . For n > 2 we 
have the obvious geometrical inequality (see fig. 10) 

X^-H/^) > (x^'H/i))^ (3.43) 



Assume 7s > 0. Eq. (|3.38|) implies that the 2-point function diverges as /i — ;► fie, i.e. the 
non-analytic part dominates and for fi ^ fic ^^ is legal to use eqs. (|3.29|) and (|3.38 
O: 



m 



>7 ^-^-^ for fi^f,,. (3.44) 



(/i - /ic)^^+"-2 - (/i - /Xc)"T^ 

We conclude that 

7. < 1 ^, n>2. (3.45) 

n — 1 

It is tempting to apply the formula for n = 2, in which case we get 7s < 0. However, 
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Sn ^=2 \ / S; •" i=2 




Fig. 10. The sum J2s / ' ' ' °'^^'^ ^^^ surfaces with n punctures is larger than the restricted sum J2s' / ' ' ' 
over over surfaces S'„ with n punctures where the surfaces S'^ are characterized by a "joint" where the 
n "bubbles" connected to the Xi's get together. Since we integrate over the position of the joint this 
effectively implies a factorization in integrated 2-point functions from Xi to the joint. 



( |3.43|) is not valid for n = 2 since there is not a unique decomposition in "joints", as 



illustrated in fig. 11 which gives the correct decomposition for n = 2. 

Let us consider the following two random surface models: In one model we allow the 
gluing of triangles such that a minimal loop-length on the surface can be two. This 
can only happen if the surface is pinched in a bottle neck consisting of these two loops. 
Such bottle necks are precisely what we have in mind in fig. 11. We call this class of 
triangulations 7^. The other class of triangulations differs from 7^ only by not allowing 
such two-loopsQ. In this class, which is denoted 7^ we can still have bottle necks, the only 
difference is that the length of the bottle neck loop will be three. In the scaling limit we 
clearly expect no difference between the random surface models constructed from 7^ and 
73 since the bottle neck loops anyway will be of the order of the cut-off. Denote the n-loop 
and n-point functions of the two models by G^ and Gp,. Up to the technical complications 
mentioned above (but rigorous for the other classes of random surface models mentioned) 
we have the following identity for the 1-point function, (or more precisely, for the one-loop 
function where the boundary consists of two links in the way described above). 



^ = ^I-^\og{l + G,) (3.46) 

The interpretation of the following: In the class of triangulations 7^ each link serves as a 
potential source of a bottle neck from which a new baby universe can grow. For a given 
triangulation T G 72 we can cut away the maximal size baby universes and close the 
corresponding two-link boundary. This will leave us with a triangulation which belongs 
to 72. In this way we get all triangulations of 72 with one boundary by summing over 
T3 and for each link either do nothing or add a whole one-loop universe, i.e. G^ itself. 
Finally we note that Ni = 3Nf/2. 

Eqs. (3.46) and (3.46) define the relation between the two models. Let us use that 



^ In terms of the dual (p'^ graphs the difference is that self-energy diagrams are excluded. In both cases 
(/)^ tadpole diagrams are excluded. They correspond to one-loop diagrams. 
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Fig. 1 1 . The decomposition of the 2- loop function for the class of triangulations in parts which cannot be 
cut in two part along a 2-loop bottle neck such that the two boundary loops are separated. It is possible 
to perform the sum over the "irreducible" components as shown in the lower part of the figure. 

7s < 1, i.e. that G^ = const. + (/i — ficY~'^'' is finite at the critical point. Differentiating 
eqs. (3.46) and (3.46) after /i and /i we get 

di,^ 3m_ (347) 

2 l+Gp 

i.e. the algebraic version of the fig. 11. The factor 3/2(1 + G^) multiplying each bottle 
neck is a combinatorial factor associated with the outgrow of the baby universes at the 
bottle neck. 

By universality the two models based on 7^ and T^ have the same critical exponents. 
Let us assume that 7^ > 0. It implies that xif^) ~^ C)0 for /i — > /ic- The same is true 
for x(/^) for /i -^ fie, the critical point of the T^ model. We conclude from (3.48) that 
/i(/Xc) > fie since xifiif^c)) = 3(1 + G'^^)/2 < 00 and we can Taylor expand the rhs of 
(3.48) around fi{fic)'- 

x(/i) ~ _ \ . ~ J: (3.49) 

fi - fl{flc> V/U - fie 

The last equation follows from (3.47) which shows that the transformation from fi to fi is 
non-analytic in fi^. At the critical point we have dfi/djl = 0, and since x(/i)/(l + Gp) is 
monotonic decreasing we have (Pfi/dfi'^ > 0: 



fi — fie = const. {ji — ii{jic) , fi — fJ'ifJ'c) = const. v^yU — fie- (3.50) 

The above line of arguments is rigorous for the hyper-cubic random surface model [11]. 
For the Gaussian model considered it is plausible but not completely proven due to the 
technical assumptions mentioned. In a later section we shall see that in more elaborate 
theories with more than one coupling constant it is possible to find a loop-hole in the 
argument, and indeed a different critical behavior. 
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The arguments presented above can easily be generalized. The important relation (3.48) 
is valid not only for xil^) but for the Fourier transformed Gfj_{p) of the two-point function 
Gfj^{x, y). Recall that 

x(/i) = y dx G^X^, y) = G^{p = 0). 

The extension of (3.48) to p ^ has the simple graphical interpretation that a momentum 
p is flowing through the bubbles in fig. 11. We have 

G..W = -^f|g. (3.51) 

2 1+Gp 

Since /i did not approach fie for fi ^ fi^we can in this region expand Gfi{p) as 

Gi^{p')=xm)-cp' + --- (3.52) 



where c is constant. The rhs of ( |3.51| ) can be expanded around p = and after the use of 



(3.48) and little algebra this leads to 

GJp) , , -, ^ o • (3.53) 

This relation shows that if'js > 0, i.e. xif^) diverges for /i — >■ ^^i we have: 

"^(/^) ~ x(Ai)~2, i.e. i^=^ls and r] = 0, (3.54) 

where the last relation follows from Fischer's scaling relation 7^ = z/(2 — r]). In fact ( |3.54| ) 
shows that for strings embedded in R^ the mass m{fi) scales to zero if and only if x(yu) 
is divergent at the critical point. 

A as final application of technique leading to (3.48) let us consider the string tension. 
We consider a large loop I a where the boundary length \Ia\ ~ \A and the sum of all 
random surfaces with this loop asboundary. Again we can cut away two-loops and we 
get directly the analogue of (3.46)t°l 

G^iU) = GpiU), i-e. or(^) = a{fi). (3.55) 

This relation tells us that the string tension does not scale to zero if 'js > since ft does 
not go to fie for fi ^ fie ^^ that case. Let us instead write: 

a(/i) = ctq + c(yu - yUc)^. (3.56) 

This still defines p as a critical exponent and differentiating ( p.55| ) after /i and using (3.47) 

■1 X(/^), i-e. p = l-7s- (3.57) 

an 



^^ It is assumed that one cannot have two-loops directly at the boundary. If this assumption is dropped 
there will be a perimeter term (1 + G^)''-*' on the rhs of (3.55). This term plays no role for the string 
tension argument. 
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Fig. 12. A typical triangulation of a square loop. When mapped into R^ we integrate over the interior 
vertices, while the boundary is kept fixed. 

If we combine this with 7^ = 1/2 we get p = 2i/ as one would expect from dimensional 
analysis if the string tension was scaling. It seems still to be satisfied with the definition 
( ^3.561 ). We will return to this definition when we consider strings with extrinsic curvature. 

3.3. Non-scaling of the string tension 

We have defined the string tension cr(/i) as the exponential decay of the one- loop Green 
function for large loops (see (|3.34|) ). Above we presented arguments in favor of a non- 
scaling string tension. The arguments did to constitute a proof since we had to make 
certain technical assumptions. It is therefore important that there exists a simple rigorous 
proof of the non-scaling of the string tension [17]. 

Theorem: a(/i) > for all /i > Pc- 

Let us for simplicity assume that we have a large square loop I12 of area L^ : 

G^{Il2) ~ e-"^^^^' for L ^ 00. (3.58) 

It is possible to bound the Green function G^{Il2) in the following way: The points at 
the boundary are kept fixed and are not integrated over, as it is also assumed in the 
general notation for the loop Green functions. It is natural to imagine that the density 
of boundary points is proportional to the length of the perimeter but it is not essential 
for the following. Let T be one of the triangulations (see fig. 12) in the sum (|3.25|) . The 
action 



S[X:T] = ^(Xi 






ihj) 



can be bounded because of the following decomposition : 

S[X, T] = SminiT, Il2) + S[x', T'] 



(3.59) 



(3.60) 



40 




Fig. 13. The surface obtained by contraction the boundary loop of a graph like the one in fig. 12 to one 
point. 

For the given (abstract) triangulation T we let SminiT, 1^2) denote the minimum of S[x, T] 
as a function of the coordinates x^{i) of the vertices i G T/dT. T' denotes the triangu- 
lation where all boundary points are identified. For the surface embedded in R^ we can 
view it as a contraction of the boundary loop Il2 to a single point 0^, of order \dT\ oc L. 



This is illustrated in fig. 13. The decomposition ( |3.6CI| ) follows from the quadratic nature 
( p39D of 5[a;,T]. 
The loop Green function G^{Il2) for the large square loop 1^2 can now be written as 



G,{1 



mI^l^ 






-SminiT,lr2)-IJ.NT 



J n dx(z 



-S[x,T'] 



(3.61) 



ieT/dT 



Next we note that the sum of squares of the length of any two sides of a triangle is > 
2 times its area. It follows that 



SminiT, 1 1,2) > 2L 



(3.62) 



and from (|3.61|) we can write 



G,{1 



A'^L^, 



<e 



-2L^ 



G^i^i 



(3.63) 



In eq. (|3.63|) G'^(Ol) denotes the loop Green function where the loop 1^2 is contracted to 
one point of order \Il'^\ oc L. Since /ic is independent of boundaries this Green function 
has the same critical point as ordinary Green functions like G^{Il2) and it can be bounded 
by 



G^(Ol) < e^(^)^ 



(3.64) 



where c(/i) is finite for /x > /ic. This is a consequence of 7^ < 1/2, which implies that the 
one-point function is finite at the critical point. 
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From the definition of tlie string tension it finally follows that 

a(/i) > 2. (3.65) 



People performing strong coupling expansions will recognize estimates like ( |3.63| ) as 
typical strong coupling estimates. What usually happens is that the function G^{Ol), 
which is based on a strong coupling approximation, becomes dominant before one reaches 
lie- However, in this case we can control it all the way down to /ic since 7^ < 1/2. 

What are the consequences of this non-scaling of the string tension? We have to assume 
that m{n) scales in order to take the continuum limit. As discussed above this uniquely 
fixes how the lattice spacing scales to zero as a function of /i: a(/i) ~ m{fi). It follows 
that the physical string tension scales to infinity since we have (for dimensional reasons) : 



apha2(/i) = cr{fi) > 2. (3.66) 

Since the physical string tension scales to infinity fluctuations including any surfaces 
having an area different from the minimal area for the given boundary 1^2 will be strongly 
suppressed. When we approach the critical point fic we will be left with a class of surfaces 
consisting of a minimal surface, depending on the Green loop function in question, and 
singular, spiky, branched polymers growing out everywhere on this surface. Such polymers 
are essentially one-dimensional objects with no or very little area. 

For large dimensions D there is little doubt that this picture is correct. It might still 
be that it is too coarse an approximation to consider the surfaces strictly as polymers for 
lower dimensions. According to the theory of such polymers the generic values of 7^ and 
u for polymers are 7^ = 1/2 and z/ = 1/4 as we will now explain. 

3.4- Branched polymers 

Let us return to eq. (3.48) as shown in fig. 11. If 7^ = 1/2 we concluded that the modified 
theory based on the class T^ of triangulations is not critical for fi ^ ficiTi the model based 
on T2. This implies that the individual bubbles in fig. 11 are not critical, i.e. they are of 
lattice size. The only way the number of triangles can grow to infinity is by the successive 
gluing of bubbles. All dynamics lie in this gluing and it seems that we get a perfect model 
of this dynamics by consider a model of branched polymers: Each individual bubble is 
represented as a link with an associated chemical potential which we denote fi as in the 
original model and with a weight factor /„ associated with the joining of n bubbles at 
a vertex Vn- In the original surface theory we have the possibility of gluing n bubbles 
to n links which share two vertices. This is the motivation for introducing the factor /„ 
associated with a branched polymer vertex f„ of order n. If we consider only spherical 
surfaces (as we will do in the following) the branched polymer graphs will be tree graphs. 
There can be no loops. The partition function for the branched polymers can be written 
as 

Z{^^) = T.^^-''''''X{fni.) (3.67) 



BP SbP 



V 
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Fig. 14. The equation for rooted branched polymers. 



where the summation is over all branched polymers, i.e all tree graphs, and the product 
is over all vertices in the tree graphs. It is convenient to consider rooted trees, i.e. with 
one marked link, since the symmetry factor Sbp in this case drops out. It corresponds to 
the one-loop function G^ considered for the full surface theory, and we will also denote 
the corresponding BP function G^. It satisfies the self-consistent equation [ 7]: 



G, 



e->^{l + f,G, + fsGl---) 



which is shown in fig. 14. We can solve for /x as a function of G 



M- 



1 + f2G, + fsGl 
G„. 



F{G,). 



(3.68) 



(3.69) 



This relation is shown in fig. 15 and for the weights /„ positiveliil we conclude that the 
lowest value /Xc of yU for which eq. ( p. 691) has a solution is the minimum of F{G^): 



/i - /Xc ~ c{G^^ 



C^fj.) ) 



I.e. 



G„ — G 



Mc 



Cy/H - He- 



(3.70) 



Since G^ is the one-point function it is expected to have the critical behavior G^ ~ 
(/x— /ic)^"'^'' and from ( p.70| ) we conclude that the generic value of'js for branched polymers 
is 1/2. 

One can consider branched polymer models directly in i?^, for instance with Gaussian 
interactions between the vertices. As long as the branched polymers are tree graphs it is 
possible to perform the Gaussian integration explicitly, as for the random walk, and for the 
calculation of 7^ one immediately gets the above considered model. Finally v = l/A and 
?7 = 0. This has the following interpretation: Consider a branched polymer propagating 
from X to y. For a given branched polymer there is a unique path of minimal length 
along the links going from x to y. This path can be viewed as a random walk path and 
the summation over branched polymers is a summation over all random walks from x to 
y where each vertex can be the source of an outgrowth of a rooted branched polymer. 
This observation allows us to solve the problem as a random walk problem, only is the 
chemical potential e~^ renormalized by the factor (1 + G^), i.e. we get a random walk 



^^ There is a loop-hole in this argument if we allow infinite branching and the weights /„ in addition 
satisfy certain conv ergence relations which move the critical point ^c out to the radius of convergence of 
the rhs of eq. ( 3.69 ). I refer to the original articles for discussion [ 7]. 
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Fig. 15. The graphic solution of the branched polymer model, 
with fi given by: 



i.e. expanding around /Xc 

1 



fi- fie 



1 



dG„ 



{jJ, - flc) ~ V/^ - /ic- 



1 + G^ d/i 
The exponential decay of the random walk is given by 



(3.71) 



(3.72) 



(3.73) 



The exponent u = 1/4 is due to the non-analyticity of the coupling constant transforma- 
tion (|3.71|) - (|3.72|) , a phenomena we have encountered a number of times by now. Since we 
still have u = l/dn we conclude that the Hausdorff dimension of the ensemble of branched 
polymers is 4, the "product" of two random walks. 

3.5. Extrinsic curvature terms (I) 

It is natural to ask if there are more elaborate random surface models where both the 
mass and the string tension scales. One can view the situation much as in the random 
walk case. The generic random walk had u = 1/2, but by tuning an extrinsic curvature 
term to infinity it was possible to enter a new universality class of smooth random walks. 
We saw that this universality class was related to spinning particles in the sense that 
one got the same critical exponents as for the Dirac particle and that the random walk 
with extrinsic curvature term could be viewed as a classical spinning particle, the tangent 
vector of the path playing the role of the classical spin. For the Dirac particle one could 
further relate it to an underlying world line supersymmetry. 

Could a similar scenario be present for random surfaces? The theory of random surfaces 
is not so well understood yet and we cannot present analytic arguments in the same detail 
as was the case for the random walk. However, in many ways the situation is identical to 
that of the random walk [ 18, 19]. First, it is possible to add extrinsic curvature terms to 
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the action in a natural way. In the context of statistical mechanics of membranes these 
terms are well known (see the lectures of Peliti and Nelson )li3. The terms suppress the 
branched polymer outgrows and their presence could result in a phase transition where the 
new phase is characterized by smoother surfaces. In the case of the random walks we had 
to take the bare coupling constant of the extrinsic curvature term to infinity in order to 
reach a new phase. For the random surfaces there are strong indications that the transition 
takes place for a finite value of the coupling constant. In addition numerical simulations 
indicate that both the mass and the string tension scales to zero at the transition point. 
Therefore this point is of interest if we want to discuss continuum limits of the random 
surfaces theory. One could further ask if there is any hint of extrinsic curvature terms 
coming from a fermionic surface theory. The answer is yes. An old result of Wiegmann 
[ 21] shows that the integration over fermionic variables will produce such terms. To be 
more precise the results are the following: If we consider a fermionic string theory, i.e. a 
string theory where we have local worldsheet supersymmetry, i.e. our bosonic variables 
x'^{^) have supersymmtric partners ip'^{^), one can explicitly integrate out the fermions 
and arrive at the effective bosonic string action: 



Seff = Skos + rj d^^Vh [{D^nfY + (e^9,e 



MA^/^^2 



+ 



— Tr (- [ d^^ dAAA + — [ d^^A A A A a) . (3.74) 

8 V2 J Sir JD ) 

The terms in eq. (3.74) refers explicitly to the extrinsic geometry of the surface defined 
by the bosonic variables x'^(^). r depends on the fermionic representation, k is an integer, 
e^, a = 1, 2 are unit tangent vectors of the surface, nf are D — 2 normal vectors and Da 
the covariant derivative with respect to the generic normal bundle associated with the 
surface: 

A^^. = nfaX, A'' = A%Mij, Mijeso{D~2). (3.75) 

The Lie algebra elements A*^ will be generators of parallel transport in the (-D-2)- 
dimensional vector space orthogonal to the tangent space. The last term in eq. (3.74) is 
a WZW-like term. / d^^ is the integration over a three-dimensional manifold which has 
the two-dimensional manifold as its boundary. This term, which in Euclidean space is a 
pure phase factor, is analogous to the term Pexpz J uj'^'^s^'^ /2 for the Dirac particle. It is 
a very interesting question whether it serves to cancel the contributions between various 
rough surfaces in the path integral, leaving only smoother ones, as was the case for the 
fermionic particle. The term (Dafif)'^ definitely acts in favor of smoother surfaces. 

To see this recall the following facts from the classical theory of surfaces embedded in 
D dimensions. Let hab denote the induced metric and F^^'^ the corresponding connection: 

hab = daX^dbX^, Tab;c = - {dahc + dbkac " d^hab) = dadbX^dcX^. 

where the indices a, b,c... are lowered and raised with hab and the inverse h"'''. The second 



^^ Also from the point of view of the continuum string theory there is a vast hterature and many different 
motivations for including extrinsic curvature terms. In [ 20] a very incomplete list of references is provided. 
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fundamental form is given by 

Ki-ab = -dau'^dbx'' = nfdadbx''. 

The basic equations which are satisfied if x^{^) represents a surface embedded in R^ are 
the Gauss- Weingarten equations: 



dad,x^ = T^,'d,x^ + K,,abnt 
dan': = -n\darfj - Ki-abh^'d^x'^ 

Using the definition ( |3.75D the last equation reads: 

Let Da denote the ordinary covariant derivative with respect to the connection T^^^'^ and 
□ = DaD"'- Using the covariant Constance of hah'- D°'hab = one can check that: 



UX^ = -l=da^/hh-'d,X^ = h'^\dad, - TJd,)x'' = Kl 






The final result is that 

[Ux^f = {h^'K,.,ab)' = {h^'DanfD.nff. (3.76) 

The principal curvatures ni and H2 of the embedded surface is determined by the second 
fundamental form: If we define the mean curvature H = {ni + fi;2)/2 and the Gauss 
curvature (as already mentioned) K = R/2 = K.1H21 we have 

{h'''Ki,abf = ^H\ {h'''K,,abf - h''^h'''Kr,abKi,cd = R/2 (3.77) 

To summarize the situation we seemingly get a number of terms which can act to 
produce smoother bosonic surfaces if we integrate out the fermionic degrees of freedom 
for a string with world sheet supersymmetry. Here we will consider only the simplest of 
these terms: the extrinsic curvature term. 

We have two different versions of the action available: one which refers exclusively to 
the geometry of the surface F embedded in R^ and one which is a hybrid between terms 
referring to extrinsic and intrinsic geometry: 

^[F;/x,A] = / dA(/i + Ai^2), (3.78) 

JF 

S[g, X- fi,X] = J d^^v^ [/i + g'^'dax'^dkx^' + Xg'^'DanfD.nf] . (3.79) 

In the last equation the induced metric enters via the equation of Weingarten. Let us 
for simplicity discuss the situation in D = 3 (the results are easily generalized to any 
dimensions) and let us choose the action given by eq. ( |3.79|) since it is easier to use in 
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/^c(A) 




/^c(A) 



l^c{X = 0) 

Fig. 16. The phase diagram for the model (3.81) 

numerical simulations. In D = 3 there is only one normal n^ to the surface and DaUi 
reduces to daU. If we want to discretize the system there is no unique way to include such 
higher derivative terms but as in the case of the random walk it can be done in a natural 
way. Let T be an abstract triangulation and let z be a vertex and [ijk) a triangle. Xi is 
the coordinate in R^ and {xiXjXk) will define a triangle A in R^ and a normal tia- For 
a given triangulation T, i.e. a given choice of equivalence class of internal metrics g, the 
discretized version of eq. ( p.79| ) is 



S[T, X] fi, A] 



(ij) 



(A>,Aj) 



(3.80) 



and the for the path integral we can write: 



(3.81) 



T(^r 



idT/iio} 



This is a theory with two coupling constants, as for the random walk with extrinsic 
curvature. Qualitatively the phase diagram looks identical to the diagram for the random 
walk except that in the two-dimensional system we have the possibility for a phase transi- 
tion for a finite value of X as shown in fig. 16. Will there be a transition for a finite value 
Ac of A? It is not known rigorously, but extensive numerical simulations seem to support 
the idea that we have such a transition and that surfaces for A > Ac are flat [ 22]. Could 
this be the transition we asked for where the string tension scales to zero together with 
the mass. Again numerical simulations seem to support this idea [ 23, 24]! 

Let us briefly discuss the possible scaling at the critical point. For A = we have seen 
that the string tension a{^, A = 0) does not scale to zero. Let /ic(A) denote the critical 
line in the (/i, A)-plane (see fig. 16) and let A/i = /i — /ic(A) parametrize the approach to 
/ic(A). Let us assume the string tension has the form: 



aifx, A) = (To(A) + c(A)A/i2'^(^) for Afi -^ 0. 



(3.82) 



47 



Since cro(O) > we need to have cro(A) — ;> for A — i> Ac in order to have a scahng of the 
string tension. Let us finally assume that 

ao(A) ~ (Ac - A)° for A ^ A^. (3.83) 

As discussed above a number of times the continuum limit is dictated by the exponential 
decay of our correlators. Here we consider two: the one- loop function G^^x{Ia) which falls 
off like e~'^^^''^^^ and the two-point function G^^\{x,y) which falls off like e""^*^^'^^'^"^'. 
There is no reason not to expect the mass ??i(/i, A) to scale to zero for A/x -^ for all 
A < Ac and we can write at Ac: 

a(/i, Xc)A = c(Ac) A/i^KAc)^ _ ^^^^^^ 

7n{fi, Xc)x = c{\c)^^^ X = mphXph. (3.84) 

The scaling limit is one where (■)ph are kept fixed while A and x goes to infinity as A/i — >■ 0. 
We can introduce two different exponents v{X) and t'(A). But as we have already seen 
for the string without extrinsic curvature they have to agree if the string tension scales: 
i/(Ac) = t'(Ac) and they will be related to the Hausdorff dimension: z/ = l/du. 

From ( p.83|) and (3.84) it is thus possible to define a consistent scaling limit for A ^ Ac 
and A/i — > such that m?{^, A) — > and (T(/i, A) — *> while the ratio stays constant. The 
only requirement is that 

AA" ~ A/i'^ (3.85) 



It should be emphasized that it is indeed possible to measure the exponents v{X), t'(A) 
and a by Monte Carlo simulations of the statistical system (3.81) [ 24]. In such simulations 
it is convenient to transform from the "grand canonical" ensemble (3.81) to a "canonical" 
ensemble where the number A^ of triangles is fixed. This is a Legendre transformation 
in /i and A^. There is no space to discuss the details, but it should be mentioned that a 
very nice feature of this transformation is that the string tension o"(yU, A) just becomes the 
expectation value of the simplest local observable, the Gaussian action itself: 

^^ ' N->oo A 

where the average is taken in the canonical ensemble and the formula is true apart from 
subleading corrections in 1/A. 

The numerical simulations lead to a value of u close to 0.25 i.e. dn ^ 4: (more precisely 
3.5 < dfj < 4.5). It is also possible to measure the string susceptibility at the critical 
point [ 25]: 7s (Ac) ~ 1/4. The same values were obtained some years ago in a related 
hyper-cubic model [ 26] and suggests universality between the two models. We will return 
to the interpretation of 7s(Ac) ~ 1/4 in the context of conformal field theories with central 
charge c > 1 coupled to two-dimensional quantum gravity. 
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3.6. Supersymmetric random surfaces 

Part of the motivation to consider the surfaces with extrinsic curvature came from the 
study of surfaces with worldsheet supersymmetry. Is it possible to implement the local 
supersymmetry at a discretized level? The answer is no. Local supersymmetry cannot 
be put on a lattice since the generators relate to translations and rotations which are 
explicitly broken by the presence of the lattice. However, if we consider the Green-Schwarz 
formulation of the superstring the situation is somewhat different. In this formulation 
there is no local worldsheet supersymmetry, but space-time supersymmetry. This is no 
problem in the discretized approach which is Euclidean invariant. In fact one can directly 
write down an action which is supersymmetric [ 27, 28]. 

Recall the continuum formulation of the Green-Schwarz superstring: The simplest su- 
persymmetric action is obtained by the replacement: 

dax^ ^W = dax^" - iOYdaO, (3.86) 

in the bosonic string action: 

/ d'^^g'^'dax'^dax'' ^ I d^^^g'^'daU'^daU'^. (3.87) 

This action possesses an obvious supersymmetry if ^(^) like x'^ is a worldsheet scalar, 
but an anticommuting space-time spinor. The global supersymmetry is generated by the 
infinitesimal transformations: 

x^" -^x^'-^{eYe-e*YO), e^e + e. (3.88) 

At the discretized level we have assigned a bosonic variable Xi to each vertex in a given 
triangulation T. We now assign an additional fermionic variable 9i to the vertex and the 
discretized supersymmetric action is: 

St[x, 9] = (xf - 4 - n^jf (3.89) 

nf^ ^ '-{ea^e^ - e^YO^) (3.90) 

The path integral would now be written as a summation over all triangulations and for a 
given triangulation the integration is over all xf s and ^j's: 



^(/^) = E ^ / n dx, d^, d^, e-^-[-'^] . (3.91) 



The problem with the above action at the continuum level is that it is not understood 
if it leads to the theory we want. Even at the classical level the continuum system seems 
quite impenetrable. The Dirac bracket prescription leads to complicated expressions which 
seem impossible to disentangle. A naive counting of degrees of freedom shows that in 
spite of the global supersymmetry there is not a perfect match between the fermionic and 
bosonic degrees of freedom. In dimensions D = 3,4,6 and 10 it is possible to add an 
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additional term to the action which siniphfies the situation drastically: 

d^^ e^^daX^ie-i^dbO - dbOYO)- (3.92) 



This term has the correct symmetry properties in the above mentioned dimensions if the 
fermions are chosen as Majorana spinors in D = 3 and 6, Majorana or Weyl spinors in 
D = 4 and Majorana- Weyl spinors in D = 10. In all of these cases the resulting fermionic 
degrees of freedom will he2{D — 2). If we naively state that the bosonic degrees of freedom 
will he D — 2, since two degrees will be absorbed in reparametrizations of the surface, we 
still have an incorrect number of degrees of freedom. It can be shown that eqs. ( |3.87| ) 



and (|3.92|) together lead to an additional local fermionic symmetry on the worldsheet, the 
so-called K-symmetry, which effectively allows a decoupling oi D — 2 fermionic degrees of 
freedom. Unfortunately the local nature of the k, symmetry makes it difficult to enforce 
at the discretized level. A suggestion for a discretized version of ( |3.92| ) for a given metric, 
i.e. a given triangulation, has been [ 27]: 

t Y. ^^3k<^%- (3-93) 

A{jjfc)eT 

However, we have no exact /t-symmetry even with this term, so from this point of view it 
is not clear that we need to add the term. 

To summarize, it is nice that one can write down a regularized version of a superstring 
partition function, which even at the discretized level has space-time supersymmetry. 
However, it is not based on the action used in the continuum. Nevertheless one should 
be aware that the continuum action ( p.87| ) was discarded not because it was wrong, but 



because a simpler alternative arose after adding the term given by eq. ( p. 92 ) to the action 



( p.87| ). The action ( p.87| ) might still serve well in a non-perturbative framework, like the 
one presented here. 



4. Matrix models and two-dimensional quantum gravity 

4.1. Matrix models 

For the bosonic string without extrinsic curvature term the results of the discretization 
was somewhat disappointing in the sense that we did not find a string theory. However, 
it is worth to recall that it is possible to calculate 7^ using continuum methods [ 29]: 

c - 1 - Jic- l)(c-25) 
7. = -^^ (4.1) 

The formula is not valid for c > 1. The bosonic string has c = D and we are probing a 
most difficult region. 

For c < 1 some results are known from Liouville theory [ 30], e.g. ( |4.1| ). The concept of 
discretization works very well for general covariant theories in the same region. A number 
of aspects of the theories can be solved both by continuum methods and directly at the 
discretized level. Historically many of the results were obtain first using the discretized 
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approach. The exact solution allows one to study in detail the scaling limit. In this 
section I will outline how to solve some of the theories with c < coupled to quantum 
gravity, using very elementary tools. 

The continuum action of two-dimensional quantum gravity coupled to matter is given 
by: 



Seh[g; fj', G] + Sm[(p, 9;M = J d^^v^ 



1 



fX 



AttG 



R{0 + Cm{<P,9;X) 



(4.2) 



Here 5'eh denotes the Einstein- Hilbert action, Cm an invariant Lagrangian density, (j){^) a 
matter field and A a coupling constant, /i and G are the cosmological and the gravita- 
tional coupling constants, respectivly. Since the integration over the curvature term is a 
topological invariant for closed two-dimensional surfaces, only the cosmological term will 
play a role, except if we sum over topologies. If we consider a specific manifold, i.e. if 
the topology is fixed and characterized by its Euler characteristic x? and if A denotes the 
area of the manifold for an equivalent class of metrics, the Einstein-Hilbert action is 



Seh[g]fJ',G] = fiA 



X 
G' 



The quantum partition function can be written as 



(4.3) 



(4.4) 



We have already considered the discretization of ( [4.4| ) in the case of D Gaussian fields (p^. 
Let us consider pure two-dimensional quantum gravity, i.e. eq. (|4.4| ) without any matter 
fields 0. Recall the discretized translation of the gravity part: 



Vol(diff) ^ ^ 



(4.5) 



S^Y^_9]^^,G] ^ fiNr- 



X{T) 
G 



and the partition function 



can be written as: 



Z[/i,G]=Ee 



-lJiNT+x{T)/G 



(4.6) 



(4.7) 



As usual the summation is over a suitable class of (abstract) triangulations. As noted in 
the last section eq. ([4.7|) is not well defined since we can split it as follows 

(4.8) 






where 7^ denotes the abstract triangulations constructed from K triangles. The number 
of inequivalent triangulations which can be build from K triangle grows faster than K\. 
For a given K the Euler characteristic x{T) > —K/2 and it is clear that the sum 
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is divergent for all /i. In the continuum it is not known how to define the path integral 
except for a fixed topology. Eq. ( [4.8| ) shows that even in the discretized approach 



where we have introduced a cut-off we can still not define such a sum in a straightforward 
manner. We can split (|4.(j|) in a sum over triangulations of fixed topology: 



Z[^^M = ll^''"'ll^-''' E • (4-9) 

X K TeTKix) 

The number of triangulations of fixed topology has an exponential bound: 

^ =Ar{K,x)=e''^^K^''~^{l + 0{l/K)), (4.10) 

TeTKix) 

contrary to the factorial bound valid when all topologies are included. As a consequence 
we can make sense of eq. ( |4.9| ) if we restrict the topology: 

ZM=E^-'''mK,x). (4.11) 

K 

For a given x the critical cosmological term is the minimal value of /i for which the sum 
( [4.11|) is convergent. This value is precisely the value /ic in ( [4.10|) and it can be shown 



to be independent of x- For a fixed topology we can try to define a continuum limit by 
approaching /i^ from above. By combining eqs. ( |4.1CI| ) and ( [4.11| ) we see that sufficiently 



high derivatives of Zy,{fi) will diverge for n —^ /ic. If we define 



\-^ /x- 


Z^ifi) 




we find for fi 


—>■ fic and n > 2 - 


-7x: 


{K\ ~ 


1 




(/i-/i,)Tx-2+«- 





(-1)" — -- —-^ (4.12) 



(4.13) 

This is an indication that large K will dominate for /x — > /i,, and that it makes sense to 
introduce a scaling parameter a such that A = Ka? is viewed as the physical area of our 
two-dimensional world. 

While these results have been obtained by the mathematicians [31] by explicit counting 
the number of ways to glue together triangles to form closed combinatorial manifolds, it 
is convenient from the point of view of physics to make this counting "automatic" . This 
is done by representing the triangles by means of Hermitian matrices [ 32, 33]: Label the 
vertices of the i^^ triangle by abstract indices ai,f3i,''yi and attach an Hermitian matrix 
(paiPi to the oriented link from Oj to /?,. In this way we can attach the scalar quantity 

0a.ft0ft7>'^7.«. = Tr 0^ (4.14) 

to each of the K triangles. The Gaussian integral 



e-iT^'^^'J-f-TrA'' (4.15) 
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^2 



Fig. 17. The matrix representation of triangles which converts the gluing along links to a Wick contraction, 
where 



a</3 a</3 



(4.16) 



can be performed by doing all possible Wick contractions of ^-fields. This corresponds 
to performing all possible gluings of surfaces of K triangles. The reason being that each 
Wick contraction in the Gaussian integral glues together two links: 



5a,ft(/)a,/3. 



)-/ 



■\\4>M^ 



^a,P,(PajPj 



Sa,/3jS/3,ay 



(4.17) 



This is illustrated in fig. 17 
eq 



After all Wick contractions are performed on the Ihs of 



( [4.15|) the K triangles have been glued together in all possible ways. The surfaces 
created in this way will consist of disconnected parts, but we get the connected graphsEJ 
by taking the logarithm of all graphs. Furthermore we can calculate the contribution from 
a particular graph constituting a closed surface: in the process of successive gluing we 
pick up a factor N, N being the number of indices, whenever a vertex becomes an internal 
vertex in the process of gluing together links by Wick contractionsEJ. This means that we 
get a total factor A^^, where V is the number of vertices. If we make the substitution 



Tr 



A^ 



Tr 



(4.18) 



it is seen that the factors of A^ for K triangles combine to A^^ ^^"^ = N^, since the Euler 
characteristic for a triangulation of K triangles, L links and V vertices is 



^ = V -L + K = V- K/2. 



(4.19) 



^^ Among the connected graphs created in such an unrestricted gluing will be graphs with one-loops and 
two-loops. These are strictly speaking not combinatorial manifolds but they have a clear identification as 
"surfaces" with a specific x- As already mentioned we do not expect such short distance phenomena to 
play any role in the scaling limit and this is substantiated by the known fact that although the value of 
IJ.C in eq. (4.10) will depend on the particular class of graphs we consider (one-loops excluded, two- loops 
excluded etc), the exponent 7^ will not. This is in accordance with the general behavior in the theory of 
critical phenomena: the positions of the critical points are not universal, only the critical exponents. 
^^ It is important to stress that A'^ is a formal expansion parameter which should always be taken to 00 
at the end of a calculation since the analogy with surfaces is based on the fact that the indices of the 
different vertices are independent. 
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By multiplying each triangle Tr 0'^ by 1/yN the weight of a given triangulation only 
depends on its topology. In addition the sum over all triangulations exponentiates. Col- 
lecting this information we can write: 

where Z{ji, G) is defined by ( |4.9| ) and 



Z{g, N)= dcp exp — Tr 0^ + -J^Tr 0^ (4 21) 




provided we make the identification: 

^=logiV, f^ = -\ogg. (4.22) 

Below it is shown that eq. ( [4.21| ) allows a l/N"^ expansion. This expansion is therefore 
identical to the topological expansion ([4.9|) of the random surfaces and eq. (|4.21| ) is an 



attempt to perform a summation of this expansion. A glance at eq. (|4.13| ) gives an idea 



of the physics involved in this summation: Since 7^ < 1/2 the partition function Z(fi) is 
finite at the critical point fic and for Z{fi) itself we get: 

qX/G 

Z(u, G) = y^ — , , ^ + less singular terms. (4.23) 

The only way we can imagine at all a summation over x for fJ' —>■ f^c is, to conjecture that 
1s{x) ~ 2 = ex and that the bare gravitational coupling constant G is renormalized: 

11,1 . ^ 

— = --log^ -. (4.24) 

Gr G (/i - Her 

Although the Einstein- Hilbert action is topological in two dimensions it will play a non- 
trivial role if we attempt to perform the summation over topologies. The conjecture 'js{x)~ 
2 = ex turns out to be true, and c > 0. This implies that the continuum limit fi —>■ fi^ 
forces G — *■ if we want to perform the summmation (|4.23|). This limit is called the double 



scaling limit. [ 34]. If we introduce the renormalized cosmological constant A by 

li- lic = Aa^ (4.25) 

this a is precisely the lattice a which goes into the discretizations of triangles by Aph = 
Nj^a"^, where A is the continuum area of the surface. Equation ( [4.25| ) tells us that that 



cosmological term has an additive renormalization. In terms of A^ the double scaling limit 
reads: 



,2c 



Na^" = eGR for a -> 0, A^ -^ 00 (4.26) 

We will later verify the conjecture 

7,(X) -2 = ex, (in fact c = 5/4.) (4.27) 
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A few comments are necessary at this point. Clearly a formula like ( |4.21| ) makes no 
sense as it stands. The matrix integral is not convergent and we have only defined and 
identified it with a sum over triangulations in a perturbative sense, i.e. by expanding 
exp((yfTr 0^/(3viV) in a power series, and performing the resulting Gaussian integral, i.e. 
gluing the triangles together via the Wick contractions. It is of course a very interesting 
question if it is possible to make sense of the integral ( [4.21[ ) in a non-perturbative way. 



The original continuum functional integral which led to ( [4.9| ) is vague when it comes to 
the question of summing over different topologies. In the case of string theory, and in the 
case of random surfaces which are intended to be a regularized version of string theory, 
we have to sum over all topologies in the way indicated in ( [4.9|) by unitarity: A closed 
string can split in two, which can later join again, in this way changing the topology 
of the surface from that of a sphere (with two boundaries) to that of a torus (with two 
boundaries). In the case of gravity it is not clear that such a change is required if we 
disregard any connection between gravity and string theory. It is nevertheless tempting to 
assume that the summation should be performed. In the general case of random surfaces 
we have already mentioned that our regularized (discretized) approach has little to say 
about the sum over topologies: It is a non-perturbative regularization of the string path 
integral for a fixed topology, but perturbative in topology. In the special case of strings in 



(i = 0, i.e. pure two-dimensional gravity, a closed formula like ( 4-21 ) seems to offer some 
possibility for a non-perturbative definition of the sum over all topologies. An obvious first 
suggestion is to define the functional integral by analytic continuation. By the rotation 
— i> e*'^/^0 we can define a matrix integral which has the same perturbative expansion as 

(EH) 



Z{g^ N) = e-^'^^'/e f ^^ ^^^ l-e'^/^ Tr 0^ + z^=Tr J . (4.28) 



Contrary to ([4.21|) this integral is well defined. If we expand the interaction term in 



powers of g and perform the Gaussian integrals we get identical results to the "Wick- 
gluing" underlying the formal expression ( [4.21 ). The problem with an expression like 



( [4.28|) is that we cannot be sure it is real. In fact it is not: it can be shown that it 



contains a complex part which is non-perturbative in the coupling constant g: it has the 
form 

lmZ{g, N) ~ e"™"^''-/^ (4.29) 

i.e. it will never show up in the perturbative expansion in g. It seems from these consid- 
erations that we have not yet succeeded in a satisfatory definition of a non-perturbative 
summation over topology, but it shows the potential power of the discretized approach 
that one is able to discuss these questions at all . 

A second, and much more simple minded remark is that the problem with the definition 
of the matrix integral ( f4.21| ) is not due to the unboundedness of the Tr 0^ term. Although 



from a geometrical point of view it is natural to use triangles as building blocks, in the 
context of two-dimensional quantum gravity one could use squares, pentagons etc.. Had 
we chosen to glue together squares, we would have a term gTr (p'^/N in the action instead 
of the cubic term, but it would still appear with the wrong sign, i.e. the action would be 
unbounded from below, since we want all "surfaces" build from squares to appear with a 
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positive weight in our functional integral. In general the gluing of n-gons will be generated 
by the matrix integral 

Z(^„, N) = l d^exp (-^Tr 0^ + ^J^Tr r) (4.30) 



For the purpose of a general (perturbative) analysis of the matrix integral ( |4.21|) it is 
convenient to consider the generalization of ( [4.30|) to an arbitrary set of coupling constants 

^(^i,^2,...)=/d0e-^^^W (4.31) 



where 



oo 



Vi{9^}) = E -'^"- (4-32) 



In eq. ( [4.32| ) we have of convenience scaled —>■ yNcf). In this way the topological nature 



of the expansion is still preserved: All two-dimensional complexes of Euler characteristic 
X will have a factor N^ associated with them. The interpretation of ( |4.31|) is intended to 
be as before: we have in mind a Gaussian integral around which we expand, i.e. (72 > 
and Qn < with the sign convention used in ( [4.32| ). The convenience of considering an 
arbitrary potential is that the general coupling constants Qn act as sources for terms like 
Tr 0", and by differentiating Z with respect to Qn we can calculate expectation values of 
"observables" like Tr 0". (Tr (p'^/N) has the following obvious interpretation: It represent 
the summation over all "surfaces" which have a n-sided polygon as boundary. This follows 
from the gluing procedure realized by Wick contractions of the Gaussian integrals. In a 
similar way 

^ (Tr 0"Tr 0^) - ^ (Tr 0") (Tr 0^) (4.33) 

will represent the sum over all connected two-dimensional complexes which connect one 
boundary consisting of n links with another boundary consisting of m links. Since two- 
dimensional quantum gravity describes the amplitude between one-dimensional geome- 
tries such expectation values are precisely what we are looking for (in the end we will of 
course have to take some kind of scaling limit in order to make contact with continuum 
physics). Let us define the generating functional for connected loop correlators. The 
expectation value of an arbitrary observable is defined by 

(Q(0)) = |/d0e-^T^^W Q(0). (4.34) 

The generating function for s-loop correlators, which we will also, somewhat inaccurate, 
denote the s-loop correlator, is defined by 

00 (Tr 0'=! ■■■Tr 0'==\ 
W{z„...,z,)^N-^-' E ^ .,.^1....,.+;^°"" (4-35) 



^fcl+1 . . . ^fcs + 1 
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where conn refers to the connected part as defined by ( [4.33| ), or its generahzation to more 
correlators. One can rewrite ( ^.35| ) as 

W{zi, ...,zs)= N'-^ (tt ^- ■ ■ ■ Tr ^-\ ■ (4.36) 

In particular, we can drop the index conn for the 1-loop correlator: 



where we have introduced the free energy F by 

Z{{gi}) = e^'^«^>» (4.38) 

and the so-called loop insertion operator by 
d °° k d 



The name "loop insertion operator" is natural since it follows from the definition (|4.35| ) 
that 

Wizu ...,Zs)= ...") ...^f . ■ ■ ■ -^T^W{z,) (4.40) 

dV{Zs) (lV{Zs-i) dV{z2) 

and this equation shows that if the 1-loop operator is known for an arbitrary potential, 
all multi-loop correlators can be calculated. 

The 1-loop correlator is related to the density p(A) of eigenvalues defined by the matrix 
integral as follows : 

p(A) = ^^<5(A-A,)^ (4.41) 

where Aj, i = 1, . . . , A^ denote the N eigenvalues of the matrix 0. With this definition we 
have 

1 (Tr 0") = r dA p(A)A", Vn > 0. (4.42) 

iV J —oo 

and therefore 



oo 



W{z) = I dA -^ (4.43) 

J — oo z — A 

For N —>■ oo exist, as we shall see, consistent solutions where the support of p is confined 
to a finite interval [y,x] on the real axis. In this case W{z) will be an analytic function 
in the complex plane, except for a cut at the support of p and it follows from Schwartz's 
reflection principle that 

2mp(X) = lim W(X + ie) - W(X - ie) (4.44) 
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4-2. The loop equations 

Amazingly few assumptions enter in the derivation of a set of equations which allow us 
to solve, as an expansion in large A^, the above defined matrix model [ 35, 36, 37, 38]. 
Let us explore the invariance of the matrix integral ( |4.31| ) under field redefinitions of the 
type: 

(j)^(j) + e<^". (4.45) 

We consider e as an infinitesimal parameter and one can prove that to first order in e the 
measure d0 defined by eq. (^4.16|) will transform like 

n 

d^ -^ d(/.(l + £ ^ Tr /• Tr 0""^). (4.46) 

fc=0 

The action will change as 

Tr V(0) ^ Tr V{(j)) + eTr 0"V^'(0). (4.47) 

We can use these formulas to study the transformation of the measure under more general 
field redefinitions of the form 



oo -A: 1 

+ ^j:ir.=^ + '7r^- (4-48) 

k=0- 



pk+1 P — <P 



This kind of field redefinitions only make sense if p is chosen on the real axis outside the 
support of the eigenvalues of (p. As mentioned above we will verify that this scenario is 
realized for N —>■ oo. Under the transformation ( 4.48| ) the measure and the action change 
as 

d0 -^ d0 ( 1 + £Tr Tr -^—] (4.49) 

V p-(f> p-(f>J 



Tr VU) -^ Tr VU) + eTi \^—V'U)] . (4.50) 

\P-(t> J 



The integral ( [4.31| ) will be invariant under a redefinition of the integration variables by 



eq. ( [4.48| ) and the change of measure and action has to cancel to first order in e. By use 
of eqs. ( [4.49| ) and (^4.50|) this leads to the following equation: 



I'^K^^J-^^ (^HJ^-^-^'-O. (4.51) 

The first term in this equation is by definition 

N'^W{p)W{p) + W{p,p). (4.52) 
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Fig. 18. The integration contour C and the cut from y to x. When deforming the contour to infinity we 
get two contributions: one from the circle Coo and one from the circle Cz around the pole z. 



The second term in eq. ( [4.51| ) can be written as an integral over the the 1-loop correlator 
by means of the density of eigenvalues: 



l/Tr^' 



dAp(A) 



du V'(uj] 



p — \ Jc 2m p — uj 



W{ij). 



(4.53) 



The second equality is obtained by the rewriting 



dAp(A) 



dcu 1 V'{uj) 
c 2TTi \ — UJ p ~ UJ 



da; V'iuj] 



c 27ii p — UJ 



dAp(A) 



\ — UJ 



where the curve C should enclose A but not p. However, for the change of the order of 
integration it is essential that we can choose C such that it encloses all eigenvalues and 
not p. This is possible if the density p has compact support on the real axis. With this 
assumption we finally can write (|4.51|) in the standard form, known as the loop equation: 



duj V'{uj) 
c 2ni z — UJ 



W{u;) = W{zf + —W{z,z 



(4.54) 



where z is outside the interval [y, x] on the real axis. In addition, since we have seen 
that W{z) is analytic outside the support of p, eq. ( [4.54| ) will be valid in the rest of the 
complex plane by analytic continuation, again provided C does not enclose z (see fig. 18). 

We can solve the loop equation by standard contour integration if we ignore the last 
term on the rhs of eq. ( [4.54| ). This is in accordance with a perturbative expansion in 1/N, 
since the correlators are normalized to order 0(1) in N. The approximation where we 
ignore the last term on the rhs of ( ^.54| ) is called the large N approximation. Translating 
back to the language of surfaces the 1/iV^ expansion is the expansion in topology of the 
manifolds, and the large A^ approximation implies a restriction to spherical topology. 
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Let us for simplicity assume that the potential is quadratic: 

y(0) = ^02_ (455) 



From the definition ( ^IST]) W{z) falls off like 

Wq{z) = - + 0{z'^) (4.56) 

and C can be deformed to oo only picking up a simple pole {zWo{z)) at z, (see fig. 18). 
The contour at infinity contributes with 

^" ' -1 (4.57) 



/Coo 27ri 2; — u; 
and eq. ( |4.54| ) reads 



zWo{z) - 1 = WoizY, Wo{z) = l[z- ^J{z-2)iz + 2)j (4.58) 



Wo{z) has a square root cut at the real axis {[y,x] = [—2,2]) and according to ( |4.44| ) the 
eigenvalue density is related to Wo{X + ie) — Wq{\ — ie) and we have that 

p(A) = :l.y(2-A)(2 + A) (4.59) 

which is Wiener's famous semicircle law. 

For a general potential we can find a solution which has essentially the same structure 
as for the quadratic potential. By deforming the contour C out to infinity we pick up a 
pole term V {z)Wo{z) at z, and at Coo we get the following integral 

I ^}^W,(.). (4.60) 

JCoc 2,711 z — uo 

which is a polynomium Q{z) in z, as is seen by expanding the integrand in powers of 1/ui: 



^ z^-'^ ^ Wk 



H-r^ w^o(c.) = E:^, (4.61) 

i=i ^ fc=i 



Only the term with total power —1 will contribute, i.e. if V^ is a polynomium of power n 
then Q{z) will be of order n — 2: 



n 



Q{z■g,■Wk) = -Y.9^ E ^'^k. (4.62) 



j=2 j + fc=n-l, 

j>0,fc>0 



It follows that the loop equation does not determine uniquely the n — 2 "moments" 
Wi, . . . ,Wn-2, and the polynomium Q{z) is consequently an arbitrary polynomium of 
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order n — 2 which enters into the resulting second order equation for VFo(-2) 
WoizY = V'{z)Wo{z) - Q{W,, z, g,), 
Wo{z) = ^{V'{z) - ^{V'{z)y-AQ{z) (4.63) 

Not every choice of Q{z) is allowed since Wo{z) cannot have a cut away from the real axis 
by construction. Let us here assume that the eigenvalue density is as close in structure 
as possible to ( |4.59| ): 



^{V'{z)y-AQ{z) = M{z) ^{z-x){z-y) (4.64) 

where M{z) is a polynomium of degree n — 2. The eigenvalue density is 



p(A) ~ M{X)^{y-X){x-X), (4.65) 



and in this sense it is close to Wiener eigenvalue distribution ( [4.59|) . Nevertheless an 



interesting critical behavior appears when zeros of M(A) accumulate near x (or y) [ 39]. 
We can write: 



Wi 



>i^) = \ {V'{z) - M{z)^{z-x){z-y)^ (4.66) 



and the requirement that Wq{z) = 0{l/z) uniquely determines M{z) since (|4.66| ) allows 
us to write 

_, , / dujM(uj) r dujV'iuj) 1 ,^ „, 

M(z)=(p -^=9 -^ I 4.67 

Jc^ 2711 z-uj J Coo 2m z -uj J(^^^ _ x){uj - y) 

as the part of the integral which involves Wq{z) vanishes. By expanding the last integrand 
in powers of l/{oJ — x) we can write: 



M{z) = Y.Mk{z-x 

k=l 



,fc-l 



f do; V'(uj) 
M,[x,y;g,] = f —- -^ -r- (4.68) 

Wq{z) can now be expressed in terms of the "moments" M^. 

1 / oo \ 

W,{z) = - [V'{z) - ^{z -x){z-y)Y. Mfc[x, y- g,]{z - x)''-'\ . (4.69) 

The endpoints x and y of the cut is determined selfconsistently by the following boundary 
conditions 

M_i[x, y, gi] = 2, Mo[x, y; g^] = 0, (4.70) 
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These two equations follows by contracting the contour in the last integral in ( [4.67] ) to C. 



In this way we pick up a pole term at z which precisely cancels V'{z) in ( [4.66| ) and we 
get yet another expression for Wo{z): 



By expanding in 1/z and using that Wq{z) contains no constant term and Wi = 1 we get 

(EB- 

4-3. Complete solution to leading order in 1/N'^ 

In principle the complete solution at spherical level is given by (|4.69| )- (|4.7CI| ). These 
equations define Wo{z) and we can apply the loop inserting operator to obtain any multi- 
loop correlator. Quite surprisingly one can obtain an explicit formula if we change the 
matrix model slightly. Instead of Hermitian matrices we use general complex matrices 
and consider a potential: 

y(0t0) = Y: ^Tr (0t0)-. (4.72) 

n=l "- 



The one-loop correlator is given by: 



^o(^)"^E ^^V^;;/^ - (4-73) 

This model has a surface representation [ 40]. Tr (0V)" represents a 2n-gon where the 
boundary links have alternating black and white colors, corresponding the and (p^ . The 
effect of Gaussian integration with respect to the complex matrices is to glue together such 
"checker-board" polygons just as Hermitean matrices glued together ordinary polygons. 
The only additional rule is that only "white" and "black" links can be glued together 
since (0^) 7^ while (0^) = and ((jy^cj)^) = 0. Such short distance differences in gluing 
should be unimportant in the continuum limit. 

We can write down the loop equations for this model. Since the potential is symmetric 
with respect to ^ —(f) we have that y = —x and we can write: 

^^P4wi.)^WH.H'^, /^^p).2, (4.74) 



c Ani z^ — uj"^ N'^ Jc Ani ^uS^ — x 

where the last equation expresses that Wo{z) r^ 1/z ioi \z\ —^ 00. 
The solution can be written as 



W< 



,{z) = ]- {y\z) - M{z)Vz^-x^) 



.., ^ f dcU UJV'(u!) 

JCoo Am (cj2 _ z^)^yuJ^ - x^ 
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If we expand the denominator of M{z) in powers of z^/uj"^ we can write: 



Miz) = J2Mk[x,g,]i 



z^-x^f-\ 



fc=i 



Mk[x,gi\ = i 



duj ojV'iuj) 



(4.76) 



Tc 47ri (a;^ — x^)'^+^/2 ' 
and the position x of the cut is determined by the last equation in (|4.74| ): 

Mo[x,(?,] = 2. (4.77) 

While these formulas look rather complicated it is a pleasant surprise that things simplify 
a lot for the higher loop correlators. In order to apply the loop insertion operator we note 
that 



d_.^A;d d x^ d 

AV{z) ^~\ ^2fc+l d^ " dV{z) ~ Mi{z^ - x2)3/2 dx^' 



(4.78) 



where 



d 



dV{z) 



E 



k d 



dV'iuo) 2ujz 



^2k + l Qg^' 



dV{z) z"^ — uo 



(4.79) 



It is now straightforward, though tedious, to apply the loop insertion operator. We find 

[38, 41]: 



Wq{Zi,Z2) 



1 



A{zl - zi) 



2N2 



\ 



^1 ^ I 2 

Z2 — x-^ \ 



Z2 — X^ 
Zi — x^ 



2ZiZ2 



(4.80) 



1^0(2:1,22, 2:3) 



X 



I6M1 /(^2 _ ^2) (^2 _ ^2) (^2 _ ^2) 



(4.81) 



Wq{zi, . . . , Zs) 



I dV "^ 1 



n 



X 



Ml dx^J 2x2 Ml ^^^ 2{zl - x2)3/2 



(4.82) 



Notice the following [ 41, 42] (see also [ 43, 44, 45]): 

(1): The above formulas are valid for any potential V. All dependence on the coupling 

constants are hidden in Mi and x. 

(2): From ( |4.76| ) it follows that dMk/dx"^ = {k + l/2)Mk+i. For an arbitrary potential 

V and s > 1 the s-loop correlator Wo{zi, . . . ,Zs) is a simple algebraic function of Zk and 

only depends on x and Mi, . . . , Ms-2- 
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(3): The same statements are true for the Hermitean matrix model, except that we 
havetwo independent end points x and y and therefore two independent set of "moment" 
Mfc and Jk defined by: 

Mk[x,y,gi] = f -—- , =, Jk[x,y,gi]= Mk[y,x,gi]. (4.83) 

Jc 27n [u — xY^^l^^ijj — y 

(4): If we iterate the loop equation after 1/iV^ this simplicity continue to hold. For genus 
g surfaces the Hermitean s-loop correlator Wg{zi, ..., Zg) is a simple rational function of 
{zk — y) and {zk — x) and depends apart from x and y on at most 2(3(7 — 2 + s) moment 
Mki Jk, k < 3g ~ 2 + s. This is even true for s = if g > 0. As an example it can be 
shown that the free energy Fg for genus g = 1 is 

Fg=i = -^logMiJi(x - y)\ (4.84) 

This formula is valid for all potentials V{(f))\ 

The proofs of the above statements are all trivial (but sometimes tedious) and involve 
nothing beyond elementary linear algebra. 

^.^. The scaling limit 

In the case of the simplest potential (Tr 0^ for the Hermitean matrix model, Tr (0V)^ 
for the complex matrix model), we have one independent coupling constant g if we fix the 
coupling constant in front of the Gaussian term. Eq. ( [4.22| ) gives the relation between the 
bare cosmological coupling constant /x and the coupling constant g of the matrix models. 
We have seen that there is a critical /ic such that the continuum limit should be taken for 
II -^ He- Corresponding to /ic there will be a gc- If we now introduce the lattice spacing a 
the area of our universe for a given triangulation is given by N^a"^, N^ being the number 
of triangles in T. It is natural to introduce the renormalized cosmological constant A by: 

fi — fic = Aa'^, i.e. gc — g ^ Aa^. (4.85) 

In this way the renormalization follows the renormalization of one- dimensional gravity 
and is in accordance with the general additive renormalization of dimensionful coupling 
constants. In case we extend our model and consider the gluing of arbitrary polygons (but 
with positive weight) we expect nothing new, except that the critical point gc will now 
be a (n-l) dimensional hyper-surface if we have n coupling constants. The identification 
of this hyper-surface is easy in the present formalism. Let us for simplicity consider the 
complex matrix model. The density of eigenvalues p(A) is given by ( [4. 651 ) and ( |4.7(j| ): 



p(A) ~ M{X)Vx^ - A2, M(A) = J2 Mk{x^ - A')^-\ (4.86) 

k 

The only non-analytic behaviour is associated with the endpoints of the distribution 
|A| -^ X. If Mi[(yfj] is positive the behavior will be identical to that of the Gaussian model. 
Ml involves a relation between the positive coupling constant in front of the Gaussian 
term and the negative coupling constants of the polygons. It can be fine tuned to zero on 
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a n — 1 dimensional hyper-surface. All higher Mk will be negative if the coupling constants 
Qk are negative. We conclude that the critical hyper-surface is determined by: 

Mi[a;(^,),^i] = 0. (4.87) 

A glance on PFo(^i; ■■■iZs) corroborates this observation since it will be singular precisely 
when Ml = 0. Let us denote a critical point on the hyper-surface by g'[ and the corre- 
sponding endpoint of the eigenvalue distribution by Xc- In statistical mechanics the masses 
and running coupling constants are defined by the approach to the critical point. We have 
already seen examples of this in the case of random walks and non-critical strings. Let 
us therefore move slightly away from the critical hyper-surface by scaling the coupling 
constants according to 



g, = g'r{l-Aa')=gf + 6g,. (4.88) 

Corresponding to this change there will be a change x"^ -^ x^ — 6x^. It can be calculated 
directly from ( f4.77| ) using dMk/dx"^ = {k + l/2)Mk+i and Mi[xc, gfl = 0. Let us introduce 



the notation: 

M^ = M,[x,, gl], i.e M^ = 2, M^ = 0, M^ ^ 0. 
Expanding Mo[x,gi] leads to: 

2 = Mo[x, gi] = M^{1 - ka^) + -M^ {5x^f + ©(a^), (4.89) 



or 



(^^') = TTiiT^ Aa^ i.e. {5x^f^5g,. (4.90) 



3Mo^ 



In addition Mi[x{gi)., gi\ will now be different from zero: 

Mi[x,g,] = ^Jx' + 0{a'). (4.91) 

By a trivial rescaling of the cosmological constant we can there define it by 

a;2 = x^ - aVX. (4.92) 

Let us now calculate su. As discussed in detail in the last section 7^ will be related to 
the s-point function byll3: 

where ( [1.85| ) is used. We get the contribution to the s-point function in the scaling limit 
from ( [4.82| ) by first performing the contour integral with some finite powers of Zi. This 

^^ In the string case we had to integrate the s-point function over aU space, but here is no target space 
and the integration simply drops out. 
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produce a non-singular function of x. Using ( |4.90| ) and ( [4.91| ) we can take a — > and we 
get: 

l^o(^i,...,g~(AaY/'-^- (4.94) 

We conclude that 'y.str = ~l/2 

Let us now turn to objects which are not readily calculated in Liouville theory, correla- 
tion functions between genuine macroscopic loops [ 41]. These are object of fundamental 
interest in a quantum gravity theory. 

From the explicit formulas for the multi-loop correlators it follows that the complex 
variable z appears in the combination z"^ — x^ and it is natural to introduce a scaling of 

z^ = xl + an. (4.95) 

With this notation we can write in the scaling limit : 

M{z) = Ml + M2{z^ - x^) + . . . = aM:^{Tx - -VA) + 0{a^) (4.96) 

TT serves the same role in the scaling limit as z at the discretized level. Knowing 
Wo{zi, ...Zs) allow us to reconstruct the multi-loop correlators consisting of discretized 
boundary loops of length ni, ...,ns by multiple contour integration. In the scaling limit 
the physical length of these loops will be li = Uia, i.e. they will go to zero. If we want 
genuine macroscopic loops we have to scale rij to oo at the same time as a — *> such that 
li is constant. By substituting ( |4.92| ) and ( |4.95| ) in the expressions for Wo{zi, ...Zs) we get 



an expression Wo{'Ki, ...tTs] A) and we can reconstruct the corresponding multi-loop ampli- 
tude Wo{li, ....,/s; A) by an inverse Laplace transformation. The reason why the contour 
integration is changed into an inverse Laplace transformation is that the cut [—x, x] (or 
[y,x] in the Hermitean matrix model) by the substitution ( 4.92|) and ( |4.95 ) is changed 



into a cut ] — oo, -\/A]. The contour integration around the cut can now be deformed: 

/ f[dz,z^^^Wo{zu...,Zs)^ (4.97) 

/ n dTTi e'»"» WoiiTi, . . . , TT,; A) , c> VA, 

Jc—ioo • 1 



c+ioo ^ 



2(?iiH hris) 



i=l 

since we have 



z"" = xf{l + avr)'/" ~ xf e"', I = n/axi. (4.98) 

The highly divergent factor x^*^'^"' ^'■^)/"-^c jg ^ wave function renormalization of the macro- 
scopic boundaries. It is to be expected since it is possible in the continuum to add a term 

SaM = ^f dl = \fd^g^O (4.99) 

JdM J 
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to the action. This is just the induced one- dimension gravity on the boundary of the 
manifold. Since A has the dimension of mass we expect in the discretized version that the 
bare couphng constant will undergo an additive renormalization (like the cosmological 
constant itself), i.e. we will have to cancel a term like x^". 

From eq. ( ^.82| ) we immediately get the generating functional for macroscopic multi- 
loop amplitudes: 

1^0(^1, . . . , ^s) ~ a^-'^/'ti;o(7ri, . . . , TT,; A), (4.100) 

Wo(7ri,...,7r,;A) = — — r ^ TT ^-— , s > 3. (4.101) 

The expressions for wo{7!', A) and wq^tti, 7^2] A) are slightly more complicated since VFo(-2) 
contains a non- universal part (y'{z)/2). and we will not give them here. 
We can calculate the inverse Laplace transform (4.97) 



pc+ioo * 

m;o(/i, ...,/.; A) = / Y[ dni e^''» u;o(vri, . . . , vr,; A) 



4-5. Generalizations 

In the last subsection we considered the spherical limit of ordinary discretized gravity. 
All graphs appeared with positive weight and the scaling limit was independent of the 
precise class of graphs used: Any polynomial 

no 

V{<P^<P) = E -(</'V)", ^1 > 0, gn<0,n>2 (4.103) 

leads to the same scaling limit. 

Let us lift the constraint in eq. ( [4.103| ). In this way we clearly move away from pure 



gravity since some of the polygons are glued together with negative weight. However, as 
already mentioned, formulas like ( [4.80D -( P:.82D are valid for any potential where gi > 0. 
Without the constraint it is possible to get a new critical behavior in the scaling limit [ 
39] since one can fine tune 

Mk[xig^),Xi] =0, A; < m, Mra[xig^),g^] ^ 0. (4.104) 



The hyper-surface in coupling constant space satisfying ( [4.104[ ) is called the m multicrit 



ical hyper- surface and pure gravity corresponds to m = 2. On this surface the eigenvalue 
density p(A) ~ M{\)\/x'^ — A^ changes from the generic Wiener form to : 

p(A) ~ (x' - a2)™+^/^ (4.105) 

Let us move away from the critical surface in a way similar to (^]8^): 



g, = g'r{l - (Aa2)"/2) = gC + 5^^. (4.IO6) 
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We want to inaintain ( [4.92| ) and ( [4.95| ) which allow us to identify a with the link-length. 
This being the case, ( [4.104| ) forces us to replace ( [4.90| ) and ( |4.91| ) with 



{Sx 



2\m 



^9i 



Mi[xi{gi),gi\ ~ {6x 



2\m~l 



(4.107) 



and one is finally led to the power a™A™'^ used in ([4.106|) . We can repeat the arguments 
used for pure gravity and find corresponding to (|4.94|) : 



Wo{zi, 



{Aa 



2\m+i 



(Aa 



2\m/2 



(4.108) 



If we compare this formula for s = to the free energy F ~ (Aa^)^ '^^ we conclude: 
3 



Is 



-m + 



(4.109) 



Let us tentatively compare this to a (p, q) minimal conformal field theory coupled to 
gravity. The central charge is c = 1 — 6(p — qY /pq and from (|4.1|) we have: 



7. 



1 

•i 

P 



q > p and co — primes. 



(4.110) 



( 4.109| ) leads to the identification of the m*^ multicritical model with a (2,2m-l) minimal 



conformal theory coupled to quantum gravity. It might be confusing that the s-point 
function Wo{ii, ...,is) does not behave like ( [4.93| ). But the explanation is perfectly in 
accordance with the identification suggested, {p, g)-models have operators of negative 
scaling dimension, the most negative being Aq = {1 — {p — qY)/'^pq- After coupling to 
gravity they might produce a more singular behavior than the cosmological term. The 
change in potential (|4.106|) excite all operators, including possible negative dimensional 
ones. These will dominate over the cosmological term for m > 2 and the scaling dimension 
of the most negative one integrated over the whole surface produces the observed behavior 
(Aa^)™"/^. This implies that each time one inserts "an arbitrary puncture" on the surface, 
i.e an excitation which involves the negative dimensional operator, it will multiply the 
partition function by A~^, where A = A*"/^. From this point of view A acts like the 



cosmological term in the unitary theories and we can rewrite eq. ([4.40|) as: 



Iyo(^l,...,^.)~(Aa'")^+- (Aa™) 



from which one would be tempted to define a 

1 

m 



It is of course a matter of definition whether one uses ([4.109|) or 7^. 



The second point concerns the expansion beyond the spherical limit. Using the loop equa- 
tion it can be performed without any problem. Let us here just mention the result in the 
double scaling limit: Approaching the m*^ multicritical hyper-surface ([4.104|) according 
to 5gi ~ 0{a^) leads to the following behavior of the moments Mk. 



Mk 



Uka 



1 < k < m. 



Mu = 0(1] 



k > m. 



(4.111) 



68 



The free energy has a genus expansion [ 46, 42]: 

^=tj^F,. n=E (".■■■".)„ ti..:;;,-. (4.112) 

where J2j aj = 3g — 3 + s and (i = x — y is the length of the cut. The second equation 
is vahd only for g > 0. For g = there are additional non- universal parts. The numbers 
(«! ■ ■ ■ «s)„ are independent of the multicritical point and can be identified with so-called 
intersection indices in the moduli space Aig^s of Riemann surfaces of genus g with s 
punctures. 

This expansion points towards two apects of the theory which we have no space to 
discuss further 

(1): The double scaling limit where A^^a^™"^^ is kept fixed should be identified with the 
renormalized gravitational constant via ( [4.26 ). In fact we get c = {2m + l)/4. 



(2): The strong indication of topological nature of the theory due to the surprising and 
beautiful appearance of the intersection indices of moduli spaces of Riemann surfaces, 
(see [ 46] for a review, and [ 45] for the connection to the so-called Kontsevich model.) 

What we would rather like to stress here is the perfect analogy to standard statistical 
mechanics. The "masses" of the theory (/i^) are fixed by the specific approach to the 
critical surfaces. These are hyper-surfaces in an (in principle) infinite dimensional coupling 
constant space and they are all of finite co-dimension, precisely as we expect in the general 
analysis of critical phenomena. 



5. The mystery of c > 1 

5.1. The Ising model 

The coupling of any matter fields to two-dimensional quantum gravity at the discretized 
level is in principle simple. First consider the theory in ordinary two-dimensional space. 
Most two-dimensional field theories can be "latticized" , i.e. they can be formulated as 
statistical models on a regular two-dimensional lattice in such a way that the continuum 
limit is recovered as the scaling limit where the "bare" coupling constants are scaled to 
a critical point. At this critical point a correlation length is divergent and one can forget 
the underlying lattice. In order to couple the theory to gravity we formulate the model on 
random lattices and take the annealed average over a suitable ensemble of random lattices 
which can be identified with surfaces of a certain topology. Usually we have in mind 
triangulations, but as we saw in the last section, it is possible to glue together a large 
variety of polygons without changing the critical behavior of pure gravity. To the extend 
that it is possible to formulate the matter theory on such polygons, we expect the same 
universality after coupling to matter. After taking the annealed average there might still 
be a phase transition and at the transition point it might be possible to define a continuum 
limit of the theory. This continuum limit will then qualify as an explicit realization of 
the original theory coupled to quantum gravity. In general the critical properties of the 
theory defined by taking the annealed average over the class of random lattices will differ 
from the corresponding critical properties of the theory defined on a regular lattice. This 
change is interpreted as the influence of quantum gravity on the matter fields. However, 
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the weight attributed to the random surfaces is influenced by the presence of the matter 
fields and at the critical point the non-analytic part of the partition function can change. 
Critical exponents like 7^ change too, and this has the interpretation as a back reaction 
of the matter on gravity. 

In sec. ^ this program was realized, although the identification with the (2,2m — 1) 
minimal conformal field theories were made in a rather indirect way for m > 2. They all 
corresponded to c < for m > 2, i.e. to non-unitary field theories. Conceptually it is 
more interesting to consider unitary theories, i.e. theories with c > 0. The Ising model 
was the first model with c > where it was possible to take the annealed average and in 
this way calculate the critical exponents of the theory after coupling to gravity [ 47, 85]. 
Later the results were confirmed by the use of continuum methods. 

Following the strategy just outlined we define the model on a regular lattice. This is 
just the ordinary Ising modeM. Next we define the Ising model on a random triangulation 
as follows: 



Z{^^,(3) = E — e-'^^- Y: exp ^E(^^^. - 1) (5-1) 

TGT^T i^^i \2(^^.^ J 

where i,j refers to the triangles in T and J2{ij) is over pairs of neighboring triangles. 
We assume that the triangulation has spherical topology. The explicit solution of this 
model is made possible by mapping it on a two-matrix model. To each triangle i is 
associated a spin o"j variable which can take two values (±1) and the model ( p.l| ) has the 
following representation in terms of triangulations: we have to glue together two kind of 
triangles (with labels ±) in all possible ways compatible a given topology (here taken to 
be spherical). The weight of the gluing along links will be 1 if the triangles are identical 
and e~^ if they are different. In addition we have the usual weight given by the total 
number of triangles. This generalized gluing process can be realized by a two matrix 
model where the two kind of triangles are represented by two different matrix potentials 
(j)^ and ip^, (f) and ip being Hermitean matrices as in the one-matrix model. In order to 
perform the gluing we use a Gaussian part: 

Tr (02 + ^2 _ 2 e-/3 0^) . (5.2) 



2(l-e-2/3^ 

The inverse of this quadratic form will be given by 
1 e-^\ 

g-/3 I J ^aif3iSa2,l32 (5-3) 

where a, (3 refers to the matrix indices, while the 2x2 matrix refers the indices ±. It 
glues triangles with different spin with weight e~^ as desired. 

^^ On a regular lattice the spins are usually placed at the sites of the lattice. Alternatively one could place 
the spins in the centers of the siniplexes or hyper-cubes etc. depending on the structure of the lattice. 
By connecting the centers we can view the latter assignment as a spin system on the dual lattice. When 
it comes to critical properties of the system one would not expect any difference since we will usually be 
interested in long distance behavior where we want to forget about the underlying lattice. Historically 
the spins on dynamical triangulations were put at the centers of the triangles and we will follow this 
convention. 
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This two-matrix model can be solved explicitly, as was the case for the one-matrix 
model. We will not discuss the solution here, only mention the result. For large /3 (low 
temperature) the system is magnetized. The spins are aligned. All triangulations have the 
same magnetic energy for this spin configuration and this means that gravity can fluctuate 
as if there where no spins at all: The fractal geometry will be as in the pure gravity case 
and 7s = —1/2. For very low (3 (high temperature) the spins will fluctuate wildly and not 
care about any underlying lattice structure either. The magnetization will be zero, the 
geometry will still be independent of the spin system and 7s = —1/2. As /5 is increased 
from zero the spin system will start to interact with geometry. This is possible since the 
magnetic energy is proportional to the length of boundaries between spin clusters. On 
a regular lattice a large spin cluster of area A will have a length L > y/A. This is not 
necessarily so on a dynamical lattice. A glance at fig. 19 shows that we can have very large 
spin clusters separated by boundaries of only a few links. In this way it is clear that the 
matter system will have a tendency to deform the triangulations towards geometry with 
small "bottle necks" like the ones shown in fig. 19. This gives us a direct visualization of 
the back-reaction of matter on geometry. For the Ising model this back-reaction is not 
sufficiently strong to the change 7^ before we reach a /3 so high that it is favorable for 
the spins to align: we have a transition. At the transition 7^ changes to —1/3. Above 
the transition it jumps back to —1/2. Not only does the spin system affect the fractal 
structure of the geometry. The fluctuating geometry will change the critical properties 
of the spin system. Intuitively one would expect it to soften the transition and this is 
what happens. The second order transition for the Ising model on a regular lattice is 
changed to a third order transition and the specific heat exponent a is changed from to 
—1. The critical exponents f3 for the magnetization and the spin susceptibility exponent 
7 (not to be confused with 7^) are also changed. Needless to say, the exponents calculated 
this way agree with the KPZ exponents of a c = 1/2 conformal field theory coupled to 
two-dimensional quantum gravity. 

5.2. Multiple Ising spins 

It is now trivial to couple many Ising spins to quantum gravity. We simply put several 
independent species of Ising spins on each triangle. On a regular lattice this does not lead 
to any new physics since they are non-interacting. The only change is that the central 
charge will be c = n/2, n being the number of Ising copies. On dynamical triangulations 
they will interact in a non-trivial way via the back-reaction on geometry as described 
above. It is clear how one can formulate the models as multi-matrix models. However, we 
cannot solve these. If the number of copies n of Ising models is larger than 2 we enter into 
the region of strong interaction between geometry and matter in the sense that the KPZ 
formula ( |4.1| ) breaks down. In this region it is still not fully understood what happens. 
This is the same region where the non-critical string is defined. However, for n — i> 00 it is 
possible to analyze what happens using mean field theorytZ. 

As a starting point we use some very important inspiration from the extensive numerical 
simulations of multiple Ising models on dynamical triangulations [ 49, 50, 51, 52]. The 



^^ The discussion here fohows [ 53], since the interpretation in terms of random surfaces is very simple. 
However, the same mean field arguments have appeared in slightly different contexts in a number of 
papers [ 54, 55, 56, 57, 58]. 
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simulations indicate that: 

1) There is still a critical point (3c below which there is no magnetization and above which 
the system is magnetized. 

2) Above I3c the geometry seems to be that of pure gravity. Below (3c there is a region 
where the situation is not clear, and where the fractal structure of the surface is very 
pronounced. For sufficiently small (3 the geometry again is that of pure 2d gravity. 

3) For a large number of spins it seems as if the system increasingly fast will be totally 
magnetized when (3 > (3c- I.e. for (3 > (3c, but quite close to (3c, the system will essentially 
be identical to the one aX (3 = oo. 

4) For a sufficiently large number of Ising spins 7s (/3c) will be positive. 

In view of point 2) and 3) it seems reasonable to attempt a description of the model 
(for a large number of Ising copies) in a region around (3c and for large (3 in terms of 
an effective model which has only spin excitations with minimal boundaries separating 
± regions since these excitations have a minimum energy. Let us for convenience of the 
following arguments specify the class of dynamical triangulations T2 used in the annealed 
average to include all gluings of triangles such that the corresponding surfaces are spherical 
and the minimal length of a closed loop of links is two. Such minimal loops are shown in 
fig. 19. For this class of triangulations a minimal boundary between spin clusters will have 
length two and we can view the surfaces as glued together of ± decorated baby universes, 
i.e. parts of the surface connected to the rest by a small loop of length two. In other 
words we consider the summation over all triangulations in the class 7^ and on a given 
triangulation we consider self-consistent iterations of the minimal spin-energy excitation. 
This is obviously identical to the first term in a low temperature expansion on a regular 
lattice. Here it will interact non-trivially with the geometry. 

Let us in this mean field model consider the one-point (or more precisely one-loop) 
function G(/x, /3), where the boundary just consist of two (marked) links. The boundary 
eliminates the symmetry factor 1/ St in eq. (|5.1| ) and we have 



G(/.,/3) = Y. ^-'""^ E'exp (f E(^^^. - 1) I (5-4) 

Tex; Wi} V fe) 



where 7^' denotes the class of triangulations where the boundary is a loop consisting of 
two (marked) links and Z^io-j' denotes the summation over the restricted class of spin 
configurations. Similarly, one defines n-point functions that are essentially derivatives 
w.r.t. fi of the one-point function. The susceptibility x{^J'^ P) is defined as 

x(/^,/3) = -^ (5.5) 

and the string susceptibility exponent 7s (/3) is given by 

x(/i,/3) = 7 — — -^^ +less singular terms, (5.6) 

(/i - fic{(3)r^"> 

for /i — > ficiP), where fidP) denotes the critical cosmological constant as a function of /3. 
The corresponding quantities in pure 2d gravity will be denoted by Xo(a^), li^^ and jjq. 
Recall that the reason why the Ising model on the dynamical triangulations could be 
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Fig. 19. A graphical representation of eq. (p.7|). The complete one-loop function allows a recursive 
decomposition into baby universes with definite spin assignment. 

solved was the ability to map it to a random surface model (gluing two kind of triangles). 
For our mean field model we have this mapping explicitly given to us by the recursive 
decomposition in baby universes of alternating spin orientation as illustrated in fig. 19. 

Using this decomposition procedure and summing first over baby universes with differ- 
ent spin configurations we obtain 



G(/i,/3)= 5:e-'^^^(l + e-2/^G(/i,/5) 



(5.7) 



rer' 



where Lf denotes the number of links in T, i.e. Lf = 1 + \Nf. In eq. ( p.7| ) the factor 
g-2/3 represents the coupling of a baby universe across the phase boundary to the rest of 
the surface and the factor 1 in the parentheses originates from the empty baby universe. 
We can now write 



G(/x,/?)= Y.^' 



-pNj 



TeT' 



G'o(/i), /x = /x--log(l 



-''G{^, 



(5.1 



Note that the last equation can be written as 
/i = /i + - log (l + e-'^ Go 



(5.9) 



which expresses /i in terms of known functions of fi and (3 since pure gravity can be solved. 
From eqs. ( ^.8|) and (|5.9| ) we get 



x(/^,/5) = Xo(/^) 






(5.10) 



dfj, 
dfj. 



e'^+Goifi) 



e'^ -ilxoifi) - Goifi))' 



(5.11) 
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Fig. 20. The phase diagram. The partition function is defined (convergent) to the right of the critical 
hne pLc{P)- The infinite volume limit is obtained by keeping j3 fixed and approaching /i(/3) from the right. 
For (3 —f oo /i(/3) approaches /ig for pure gravity. j3c denotes the phase transition point. 

It is clear that the derivation and the equations are very similar to the ones used to prove 
7s = 1/2 for the non-critical strings. Indeed, we will find the same results here, but the 
additional coupling constant (3 allows a new critical behavior as will be clear. 

Since the string susceptibility exponent 7^°^ = —1/2 < in the case of pure gravity 
both Go(/io) and Xo(a*o) are finite. This implies that a /3c exists such that the denominator 
in (|5.11|) is different from zero for all /i > jJiciP) provided P > Pc ■ Pc will be the phase 
transition point. Fig. 20 shows the phase diagram in the (/z, /5)-plane. 

Let us first consider region I in fig. 20. The only chance for a singular behavior is 
that fl{fj,c{P) = A*o and in this case the leading singularity has to come from Xoifi)- We 
conclude: 



/i(/ie(/9)) = /io and 7,(/3)=7f for (3 > P,. 



(5.12) 



This is the phase where the model is magnetized, where the spin fiuctuations are small 
and where the geometry of the surfaces is not affected by the spins. 

The number P^ is characterized by being the largest P for which dfi/dfl equals zero for 
some fl, i.e. (by (|5.11|) ) the largest P for which the equation 



.2/3 



-Xo(/^)-Go(/i) (= E(3A^t/2 



l)e 



-p.Nj 



has a solution. If we define fidP) by 
= for P<P,, 



dp. 



Ac(/3) 



(5.13) 



(5.14) 



then ficiP) obviously solves eq. ( |5.13|) and we have 

PciPc) = /^o, PciP) > yUo for p < pc. (5.15) 

Let us now consider region II in fig. 20 where P < Pc- If we use eq. ( |5.10|) and (|5.11| ) 
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this implies that the Ihs will be singular for /i -^ fJ'dP)- In fact, since fidP) < A^o both 
Xoifi) and Go{jl) will be regular around fidP) and we can Taylor expand the Ihs of eq. 

X(/X, P) ~ —^7^ ~ / "" ,^, (5.16) 

To derive the last equation we have used ( p.ll| ) and (|5.14|) which tell us that 

/i = /^c(/3) + const.(/i - lic{(3)f + ■ ■ ■ . (5.17) 

We conclude that 7s(/?) = 1/2 for (3 < (3c- In this phase baby universes are dominant. 
Effectively we have branched polymers and the total magnetization of the system is zero. 
Let us finally consider the system at the critical point (3c- This point is characterized by 
the fact that ftdPc) coincides with fiQ- Although the singularity of x(yU, Pc) for fi — > fidPc) 
is still dominated by the zero of d^/dji we can no longer Taylor expand fi{fl, Pc) around 
ficiPc) (= /Uo) since the functions in ( |5.11| ) are singular in /iq. On the other hand we can 
use the known singular behavior of the pure gravity functions Go and xo at /xq to deduce 
from ( |5.11| ), remembering that 7^*^-' < 0, 

^ ~ (/x - ficm)-''''\ i.e. /i - ficiPc) ~ (A - MPc))-''°'^'- (5.18) 



From (|5.10| ) we finally get, using (|5.18| ) 



xifJ^,Pc) ~ (0) ~ ,^,, (0),, (0)^,, . (5.19) 

and we have derived the remarkable relation: 

-7(0) 1 

%m = J; , = TT- (5.20) 

The model can easily be extended to include the coupling to a magnetic field and 
one can explicitly verify that the system is magnetized for P > Pc and has zero total 
magnetization for P < Pc- Further the transition is a third order transition (like in the 
full Ising model on dynamical triangulations) . The model captures the essential features 
observed numerically for a large number of Ising models on dynamical triangulations. It 
also strongly suggests the existence of a region below Pc where the surfaces consist of 
numerous totally magnetized baby universes with different spin orientations such that 
the total magnetization is zero. The spin transition observed is just the final alignment 
of the spins of the baby universes. 

The following should be noticed: 

(1): For a finite number of Ising models the mean field approximation above clearly fails 
for /3 — >■ 0. In the infinite temperature limit 7^ has to return to the value —1/2. In 
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the mean field approximation 7^ stays equal 1/2. For high temperatures there will be 
important spin configurations which are different from the ones considered in the mean 
field approximation. If we for sufficiently high central charge c assume the existence of 
an interval below (3^ where [3^ = 1/2, it is an interesting and unsolved question how the 
transition to 7^ = — 1/2 takes place. Is there a single transition, a continuous change or 
a cascade of transitions, and how are these transitions characterized? 

(2): The transition at Pc had < 7^ < 1/2 and it was related to 7^°^ for pure gravity by 
eq. (15.201) . The formula refiects that the individual baby universes are totally magnetized, 
i.e. have the fractal structure of pure gravity. This can be generalized [ 59]: 



Theorem: given a multiple spin model with < 7^ < 1/2 it follows that 

7s = -^^ = -^, In = --. (5.21) 

-7„ + 1 n + 1 n 



The interpretation of eq. (|5.21| ) is that the individual baby universes have a fractal 



structure corresponding to a unitary conformal field theories with c < 1 coupled gravity. 
Such theories are characterized by a c^ = 1 — 6/n(n+ 1) and a 7^ = —1/n. Formula (|5.20| ) 
represented the simplest example, c = 0, but it can be shown, using renormalization 
group arguments, that (|5.21|) is the only possible critical behavior if 7^ > 0. It rules out 
the possibility that 7^ can change continuously while larger than zero. It also tells us 
the possible solution to the c > 1 disaster in the continuum approach: The interaction 
between matter and geometry becomes so strong that the surfaces disintegrate dynamically 
into baby universes which individually have c < 1. What is missing in our understanding 
is how to predict a specific c„ < from c > 1. The only thing we know is that for c —>■ 00 
it follows that c„ = 0. 

5.3. Random surfaces with extrinsic curvature (II) 

Let us return to the random surface theory with extrinsic curvature. Recall that numerical 
simulations indicated that 7^ ^ 1/4. Can we understand the theory in the general setting 
outlined in the beginning of this section? To investigate this possibility consider the model 

on a fixed regular triangulation T: 

Zr{\)=f n dx, e-^(-)^^^-^'^-^'-'^-«.-/""'-""^-^\ (5.22) 

ieT/{io} 



In ( [5.22| ) i denotes a vertex in T which is mapped into Xi G R^ (but it is easy to formulate 
the model in R^) and wa is a normal to the oriented triangle {xiXjXk) where A = (ijk) G 
T. 

Eq. ( p.22| ) defines the model of so-called crystalline surfaces [ 60] (see the lectures of 
Peliti and Nelson for details). Since T is regular we can view it as a lattice theory, either 
in terms of the variables Xj, in which case it is non-polynomial, or in terms of the normals 
riA after integrating out some of the degrees of freedom of the Xj's. If we go to the limit 
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A = it is a purely Gaussian model and it can be solved trivially. Viewed as a lattice 
model in terms of the normals it is a weird theory [ 61]. If we use a continuum notation 
and denote the lattice coordinates by ^ the normal-normal correlator is: 

for ^ 7^ 0. In this formula A is just a convenient continuum proper-time cut-off instead of 
the lattice spacingEJ. Except at coinciding points the normals are anti-correlated. This 
truly remarkable situation can only appear because the surface is totally crumpled in R^. 
It is easy to calculate the average radius of the surface in i?^ : 

(r)^^~logiVr, (5.24) 

i.e. the Hausdorff dimension dn = oo. 

With increasing A the short distance correlation between the normals becomes positive, 
but the long distance tail remains negative, which implies that dn = oo. However, for 
a finite value A = Ac the system undergoes a phase transition, the so-called crumpling 
transition, after which the normal-normal correlation becomes positive everywhere and 
the Hausdorff dimension finite. The transition seems to be second order. Large scale 
computer simulations have verified the existence of this transition and strongly suggest 
that it is second order [ 62, 63, 64]. Calculations of the /3-function /5(A) for continuum 
toy models related to the model defined by eq. (|5.22|) lead to [ 65] 

/3(A) = ciA - C2A^ Ci > 0, (5.25) 

i.e. an ultraviolet fixed point and applying the renormalization group we would have the 
flow shown in flg. 21. There should be a central charge c associated with this transition. 
Its value is not known, but it would be most interesting if c > 1. If A = we have 3 
Gaussian fields and c = 3. Since A = is an infrared fixed point (see fig. 21) one would 
naively expect from the c-theorem that Ccmmp > 3. 

The summarize: Although weird, the model fulfills the requirement for a spin model 



we can couple to quantum gravity. It is done by taking the anneal average of ( |5.22| ) over 
random triangulations (of the sphere), i.e. we arrive precisely at the random surface model 
studied earlier. Numerical simulations showed that this model had a phase transition to 
smooth surfaces. It is natural to conjecture that it is the quantum gravity version of the 
crumpling transition. If Ccmmp > 1 we will be in same situation as for the multiple Ising 
models. The numerical measurement of 7^ gives some support to this idea. It would 
be most interesting to have a better analytic understanding of the transition since it 
is relevant both from the viewpoint of string theory, as described above, and from the 



^^ To be more precise the continuum Gaussian action is used and, if A is an ultraviolet cut-ofF and L an 
infrared cut-off, the propagator is 

JL-2 



77 




A 



Fig. 21. The /3-function for the crystalline surface model (^.22) 



viewpoint of membrane physics: It is a model of a fluid membrane, not entirely physical 
since it is self-intersecting, but still of considerable interest. 



6. Euclidean quantum gravity in rf > 2 

6.1. Basic questions in Euclidean quantum gravity 

Until now we have considered various aspects of two-dimensional quantum gravity. As 

long as we restrict topology we deal with a well defined theory. There might even be some 

hope that one could define the summation over all two-dimensional topologies although 

this question is not yet settled. At the moment we address Euclidean quantum gravity in 

dimensions higher than two a frightening number of basic questions appear. Let us just 

list some of them: 

(1): How do we cure the unboundedness of the Einstein-Hilbert action in d > 27 

(2): Does the non-renormalizability of the gravitational coupling constant not spoil any 

hope of making sense of the theory? 

(3): What is the relation between Euclidean and Lorentzian signature in the absence of 

any Osterwalder-Schrader axioms to ensure that we can rotate from Euclidean space to 

Lorentzian space-time? 

(4): What is the role of topology, keeping in mind that four- dimensional topologies cannot 

be classified? 

It is possible to take these questions and our inability to answer them in a fully satis- 
factory way as an indication that there exists no theory of Euclidean quantum gravity for 
d = 4. Let us briefly discuss the problems. While (1) and (2) in principle are unrelated 
the solution of one of the problems often solves the other problem too. The unbound- 
edness of the Euclidean Einstein action is due to the conformal mode. An obvious way 
to cure this problem is to assume that the Einstein-Hilbert action is only an infrared 
approximation to the complete theory. A simple minded way to make the action bounded 
is to add higher derivative terms like R^ or R^j^R^^ . Such terms also tend to cure the 
problem of renormalizability since the propagators will contain higher inverse powers of 
the momentum. In fact adding suitable combinations of R^ terms one can get a theory 
which is renormalizable and where the Euclidean action is bounded from below [ 66]. It 



78 



is not known whether the theory quahfies as a unitary field theory if we apply the stan- 
dard rotation back to Lorentzian signature. A more general discussion along the lines 
of viewing the Einstein-Hilbert term as the first (infrared leading term) term of the full 
action goes back to Weinberg who introduced to concept of asymptotic safety [ 67] . The 
idea is that we should be able, by means of the renormalization group, to work our way 
back from the long distance region where the Einstein-Hilbert action is a good effective 
description to some non-trivial ultraviolet fixed point. In addition the associated critical 
surface is assumed to be of finite co-dimension, which means that only a finite number of 
parameters need to be fine tuned to reach the critical surface and from this point of view 
the theory will not differ in spirit from ordinary renormalizable theories. The effective 
Lagrangian description of the field theory by means of fields suitable for the infrared fixed 
point might then be an infinite series 



^ = \fg 



iOTTG 



(6.1) 



which might even be non-polynomial, but which might make sense both with Riemannian 
and Lorentzian signatures. 

A weakness in such a scenario is that the existence of the ultraviolet fixed point is 
hypothetical up to now and we have not exactly been fiooded with examples of non-trivial 
fixed points in four-dimensional field theory. However, the interesting results obtained by 
use of the expansion in 2 -|- e dimensions give some support to this idea [ 68]. 

Finally one could hope that once the correct non-perturbative formulation of the quan- 
tum theory has been found these problems will resolve themselves and it will become clear 
why one should not necessarily think in terms of perturbative expansions around fixed 
background metrics. 

Not much is known about the question of summing over topologies. As long as we think 
in terms of continuum physics and write down the path integral it offers the possibility 
of summing over different manifold structures as well as integrating over inequivalent 
Riemannian structures of a given manifold: 

Z = \^ f ^^M^ Q~S[g] (Q2) 

In two and three dimensions we do not have to worry about the meaning of Top, since there 
is equivalence between smooth manifolds and topological manifolds. In two dimensions the 
manifolds are uniquely characterized by their Euler characteristic x ^T^d the summation 
over Top is simply a summation over x or the genus g = 1 — x/2 of the surfaces. In 
spite of this simple prescription surprisingly little progress has been made in defining the 
sum in (|6.2| ) using continuum methods. Matrix models gave non-perturbative definitions, 
as described in sec.^. However, so far the results have been ambiguous. If we move to 
three dimensions we encounter a slight classification problem in the sense that no simple 
parametrization like x of the various topologies exists. However, the problem seems 
to get completely out of control when we move to four dimensions. For four- dimensional 
manifolds there is not equivalence between smooth and topological structures. Topological 
manifolds exist which do not admit smooth structures and some topological manifolds 
admit infinitely many inequivalent smooth structures. If we insist on summing over all 
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smooth structures J^Top will be rather unwieldy. To complicate the operational meaning of 
J2top it should be added that four-dimensional manifolds are not algorithmic classifiable, 
i.e. no finite algorithm in the sense of Turing exists which allows us to decide if two 
arbitrary four dimensional manifolds are equivalent. On the other hand arguments (not 
known to me) might exist which dictate a restriction of the allowed class of manifolds. 
Since there seems to be fermions in the world one could argue that the manifold should be 
a spin manifold. If one makes the additional (rather arbitrary) restriction that it should 
be simply connected it is possible to show that such a (smooth) manifold is characterized 
by its Euler number and its signature where these in four dimensions are given by: 

X(M) = j^ j^ d\^ Ra,cdRa'Wc'd'e'''^'''e'''''' (6.3) 

t{M) = ^ J^^ d'^^ R,,,,R'^\,,,e^''^''' (6.4) 

For simply-connected spin manifolds the signature r is a multiple of 16, while x is an 
integer > 2. Thus eq. (6.4) seems a minor extension compared to two dimensions where 
the summation is over x, but the restriction to simply connected spin manifolds is not 
natural at this stage of a quantum theory of space-time. 

While all these problems seem to discourage any attempt to make sense of the path 
integral of Euclidean quantum gravity, it is still our obligation to try to investigate if it 
is possible. Below I will argue that the use of dynamical triangulations, which works so 
well in two dimensions, allows us to discuss several of the above issues in some detail and 
might be a candidate for non-perturbative definition of quantum gravity even in higher 
dimensions. 

6.2. Definition of simplicial quantum gravity for d > 2 

Let us define simplicial quantum gravity as a generalization of the construction in sec. ^ 
and sec.^ for two dimensions: In d dimensions (where d = 3 ot d = 4) we construct 
all closed (abstract) simplicial manifolds from K d-dimensional simplexes. As in two 
dimensions we imagine that the length of the links in the simplexes are a (which we take 
as 1 unless explicitly stated). For such a combinatorial or, equivalently, piecewise linear 
manifold we can apply Regge calculus and in this way assign a Riemannian metric to the 
manifold. By such an assignment we see that the discretization is able, for finite K, to 



assign a meaning to the sum and integral in (|6.2|) [ 69, 70] 



The discretization has the same virtue as in two dimensions: In principle it allows a 
unified treatment of the summation over topologies and Riemannian structures. 

A few remarks should be said about the formula ( |6.5| ). First one could try as in two 
dimensions to make the gluing automatic, and in this way arrive at a generalized matrix 
model [ 70, 71, 72]. In d dimensions the task is to glue together ((i-l)-dimensional sub- 
simplexes. Let us discuss the situation in d = 3 (the generalization to higher dimensions 
is trivial). In fig. 22 we have shown a building block for three-dimensional simplicial 
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Fig. 22. Labels of a tetrahedron 

complexes: a tetrahedron (the generahzation of the equilateral triangle used up to now) 
where the edges has labels a, j3, 7, (5, e, p, which takes values 1, ..., A^, where N is an abstract 
label which we in the end want to take to infinity like in the two-dimensional case. To 
each oriented face of the tetrahedron, a,/?, 7, say, we assign a complex variable (pap-y 
whose value is invariant under even permutations of the indices but conjugated under 
odd permutations (which correspond to a reversal of orientation). To the tetrahedron we 
associate the generalization of Tr 0"^ for the triangles: 

^(0) = Y^ 4>al3^<t}-ySe4>epa4>f3pS (6.6) 

af^'ySep 

Consider the following integral: 

Z{g)= f d<Pexpl-W\<Po.pi\" + 9A{(P)\. (6.7) 




If the exponent is expanded in a power series in g, and all possible Wick contractions are 
performed on the powers of A{(j)) it corresponds to all possible gluings of the tetrahedra 
along triangles. The expansion leads to 

log^(^)~E^^''^^'^^^''^^^^ (6-8) 

where the sum is over all possible gluings T. Ni{T) denotes the number of links in T and 
Ns{T) the number of tetrahedra. Let us postpone the discussion of the terms N'^'^^'^'g'^^^'^' 
and concentrate on J2t- This kind of "blind" gluing of rf-dimensional simplexes along 
their d — 1 dimensional sub-simplexes will not in general create combinatorial manifolds, 
but only so-called combinatorial pseudo-manifolds: i.e simplicial complexes where d — 1- 
dimensional sub-simplexes which are not boundary simplexes are contained in precisely 
two (i-simplexes and where it is possible to connect any two d-simplexes by a sequence 
of d-simplexes, each intersecting along some d — 1 simplex. For such pseudo-manifolds 
the neighborhood of a vertex will not necessarily be equivalent with a rf-dimensional ball. 
Whether or not this is important for quantum gravity is not clear. In the two-dimensional 
case it did not course any problems. One could create one-loops and two-loops which can- 
not be present on a combinatorial manifold, but the pseudo-manifold still had a transpar- 
ent surface interpretation and the Euler number was perfectly well defined. In the three 
dimensional case the modification is more severe in the sense that the Euler characteristic 
can be different from zero, while x = for any odd dimensional manifold. In the following 
we will (to be conservative) restrict ourself to so-called combinatorial manifolds, i.e. com- 
binatorial pseudo manifolds where the neighborhood of any sub-simplex is homeomorphic 
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to a (i-siinplex. For closed manifolds this enforces certain relations between the number 
of i-simplexes iVj in the triangulation, known as the Dehn-Sommerville relations: 

A^. = E(-1)'^'(- + J)^^- (6-9) 

If we define A^_i = Xi the Euler characteristic of the manifold, ( p. 91 ) is valid for i = 
-l,...,d-l. 

Secondly, in light of the complicated relation between topology and diffeomorphism 
for four- dimensional manifolds one could be worried that similar problems arise when we 
compare combinatorial, i.e. piecewise linear, manifolds and smooth manifolds. However, 
for dimensions d < 7 we have equivalence between piecewise linear and smooth structures. 



Whatever subtleties might be involved in defining the sum on the Ihs of eq. (3^) it should 
be captured at the rhs of eq. ( |6.5|) . Of course eq. ( |6.5|) itself is rather formal as it stands. 
Even in the two-dimensional case the precise meaning of the Ihs was not clear from a 
mathematical point of view and in sec.^ we tried to use the rhs to define what is meant 
by the summation over all two-dimensional manifolds. The problem was non-trivial in 
the sense that the number of triangulations grows as fast as i^!, K being the number of 
triangles. No reasonable discretized action could kill this entropy factor [ 69, 70] and 

Z = J2^~^^ (6-10) 

T 

will be divergent. This problem will of course be present also when we move to higher 
dimensions. The naive expression (|6.1CI| ) makes no sense. In two dimensions the way 



out was to fix topology. This bounds the number of triangulations exponentially and eq. 
( |6.10|) will be well defined. Only afterwards the topology is allowed to fluctuate and the 
double scaling limit taken. In higher dimensions it is sensible, as a minimum, to try to 
define eq. ( |6.10| ) for a fixed topology. We need the following conjecture : 

Conjecture: The number of combinatorial equivalent (i- dimensional manifolds is an 
exponentially bounded function of the number of d dimensional simplexes. 

We call two simplicial complexes combinatorial equivalent if they have a common subdi- 
vision and when we talk about equivalence classes of piecewise linear manifolds we have 
in mind combinatorial equivalence. While the conjecture is true in two dimensions, there 
is no general proof in higher dimensionsli^l. For d > A manifolds exist which are not al- 
gorithmic recognizable, the reason being that any finitely presented group can appear as 
the fundamental group for some four-dimensional manifold. However, finitely presented 
groups cannot be algorithmically classified due to the famous "word problem". As a 
curiosity we can mention that for such manifolds their number as a function of A^4, the 
number of 4-simplexes, is not algorithmic calculable. This does not mean that the number 
cannot be exponentially bounded, only that there is no finite algorithm which allows us 
to calculate the exact number. In fact numerical simulations seem to work well for the 
simplest topology (S^) and the give support to the conjecture for S^ [ 75, 77, 78]El. 



19 Very recently there has appeared a proof [ 73] . 

^^ Recently some controversy arose [ 79], but my opinion the question is now settled. 
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After these general remarks it is natural as a first explorative step to fix the topology 
of our four-dimensional manifold to be the simplest possible, that of S'^. The Einstein- 
Hilbert action for a simplicial manifold can be calculated by Regge calculus. However, we 
do not need the full machinery in our case where the simplexes are identical and all link 
length equal. The Regge version of the Einstein action is the sum over deficit angles of 
the d — 2 dimensional sub-simplexes times their rf — 2 dimensional volume. In our case the 
deficit angle associated with a rf — 2 dimensional sub-simplex nd-2 is 27r — const. o{nd-2), 
where o{nd-2) is the order of nd~2, i-e. the number of d dimensional simplexes of which 
nd-2 is a sub-simplex. It follows that 

/ d^^v^i?) oc 5] (27r - const.o{nd-2)). (6.11) 



"d-2 



If we note that the number of li — 2-dimensional sub-simplexes in a d- dimensional simplex 

is d{d + l)/2 we can write: 



%;A,G] = |d^ev^(A-^i? 



Srikd, kd-2] = kdNd{T) - kd-2Nd-2{T), (6.12) 

where the piecewise linear manifold is given defined by the triangulation T and Nd{T) and 
Nd-2{T) denote the number of d- and d — 2-dimensional simplexes in the triangulation T. 
We can view l/kd-2 as a bare gravitational coupling constant. At first sight the action 
might seem much too simple to have anything to do with gravity. Our point of view will 
be the opposite: The fact that the action is so simple reflect the beauty and simplicity 
of quantum gravity, and hopefully this simplicity will be reflected in the solution of the 
theory. 

Our final prescription will be: 

^[^'^]-y5^voi(dii)" 



Z[kd-2.kd]= E e-^-l'^^--^'*! . (6.13) 

res'* 

We will be interested in c? = 4 and sometimes in c? = 3, in which cases we get 

Z[A;i, k^] = E e-''^mT)+k,N, (g^^4) 

Z[k2, k^] = J2 e-'=4^4(^)+^2^2 , (6,15) 

Te54 

For d = 3,4 there are only two independent coupling constants as long as we only want an 
action which depends on global quantities like Ni. It follows from the Dehn-Sommerville 
relations. If we include higher derivative terms in the action we will certainly loose the 
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Fig. 23. A hypothetical phase diagram for four-dimensional gravity. 

simplicity of eqs. (6.14) and (6.15). The higher derivative terms will contain explicit 
the order of the sub-simplexes which carry the curvature. We will not discuss the lattice 
implementation of these. 

Let us now discuss the phase diagram. Assume d = 4. Since it is easy to prove that 
N2{T) < const. A^4(T), the conjecture above implies that for each ^2, a ^4(^2) exists such 
that the Ihs of ( |6.13|) for a given k2 is well-defined for k^ > /c^l^a) and divergent for 
ki < ^4(^2)- A potential continuum limit should be taken as k^ -^ k^ from above. The 
corresponding phase diagram is shown in fig. 23. If we for a fixed /c2 approach kl we 
will be probing the infinite volume limit of the discretized system. It does not imply 
that there necessarily will be a continuum limit. Rather we should view the system as 
a lattice system where the infinite volume limit is taken. For some specific values of the 
bare couplings critical points might exist where a correlation length diverges and where a 
continuum limit exists. Such a point is tentatively indicated at the figure. Approaching 
this point in a specific way will then define the renormalized cosmological constant and 
the renormalized gravitational constant. Since we are in unchartered territory one should 
be open-minded for other possibilities, e.g. the possibility that a whole range of k2^s might 
correspond to a topological gravity where the metric, and correspondingly concepts like 
divergent distances, play no role. 

Let us for a moment return to the three-dimensional partition function given by eq. 
(6.15). Comparing with eq. ( |6.8|) derived from the tensor model one has the identification: 



A;3 = -log5f, A;i = logiV. 



(6.16) 



This relation is the same as we encountered in two dimensions and the relation will be 
valid for the matrix model in any dimension: 



kd = - logg, kd-2 = log A^. 



(6.17) 



From a formal point of view the N —>■ 00 limit corresponds to taking the "bare" gravi- 
tational coupling constant G ^ 0. In two dimensions this limit played a very important 
role: since the Einstein- Hilbert action in two dimensions is topological, an expansion in 
1/G automatically becomes an expansion in topology. However, in three and higher di- 
mensions we have not reason to expect the large A^ limit to classify topologies for us, and 
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it does not!. It has not yet been possible to construct anything hke the double scahng 
limit in higher dimensions than two. 

6.3. Observables 

It is possible to define the same critical exponents for higher dimensional simplicial quan- 
tum gravity as we have already defined for the various two-dimensional theories we have 
considered. 

First we can define an entropy or susceptibility exponent 7^. In the following it will 
always be assumed that the topology is spherical, i.e. the combinatorial manifolds are 
combinatorial equivalent to the boundary of a 5-simplex. According to the conjecture 
the number A/'(A^4) of triangulations which can be constructed from A^4 4-simplexes is 
exponentially bounded. Let us now fix ^2. According to the remarks above there will be 
a critical point ^4(^2)- This implies that 

T£S^(N4) 

where fk2{N'i) stands for subleading corrections, //the subleading correction is power-like 
we define 7s(A;2) by 

U{Ni) ~ Ar4'^('=^^-'e'^4^4(l + 0{l/N)). (6.19) 

However, it is not clear that we have such a behavior. It is possible to imagine an 
asymptotic behavior like: 

M{Ni) ~ e''lN,-c{k,)Nnk2)^-^i + 0{l/N)) (6.20) 

where < a < 1. In this case the exponential correction given by a will always dominate 
over the power-like correction determined by 7^(^:2). There are strong indications that 
there are several regions with different asymptotic behavior, depending on the value of 
k2. This will be discussed below. 

Apart from the entropy- or susceptibility exponent 7^ we can introduce the critical 
exponents v and 77 already discussed numerous times in the context of strings and random 
walks where they were determined from the properties of the two-point function in target 
space. However, here these quantities will refer to intrinsic geometry. Let us define the 
two-point function as follows: Consider the ensemble of combinatorial manifolds with 
two marked vertices ii and 12- For each manifold, build of regular 4-simplexes, we have 
by Regge's prescription a metric. This means that we can define the geodesic distance 
between the vertices ii and 12- As a rough definition one can consider paths alone the links 
in the triangulation and define the geodesic distance as the shortest link-pathlHiJ between 
ii and 12- Let us consider the sub-ensemble of manifolds where the geodesic distance 



^^ Alternatively one could talk about the distance between the 4-siniplexes, defined as by moving between 
successive 4-siniplexes via their common 3-simplex boundary. Again this is not a true geodesic distance, 
but should not differ drastically from the genuine geodesic distance. For a given combinatorial manifold 
there might be a significant deviation between link distances and 4-simplex distances, but we expect 
concepts like Hausdorff dimensions etc. to be the same. 
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between ii and 22 is fixed to be r. If we use the link-paths in the definition of geodesic 
length r will be an integer. Let us denote this ensemble of manifolds by T(2,r). We 
define the two-point function by: 

Gi^;\r)= Y. e-^^I'^*''^^! (6.21) 

Ter{2,r) 

Provided that we have the asymptotic behavior (\6. 1 i\ ) we expect a generic behavior of 
this spherical two-point function (we suppress the explicit dependence on /C2): 

Gfc,(r) ~ r-"e"'"('=4)^ r > llm{k^)- (6.22) 

G'fc,(r)~r^-^ l<r<l/m(fc4); (6.23) 



~ r~"e"''"^''''^'' 




~ r^-^ 




- \ " n, (r\ r^ 


1 


- Z^<-^fc4l^J ~ /, 


hC\'y, 



The short distance behavior of Gk^ir) arises from the generic behavior since there is an 
angular average involved in the definition. 
If m{k^ scales to zero at the critical point k^'- 

m{ki) ~ {ki - klY (6.25) 

we can define in a simple way the continuum limit: Introduce the continuum parameters 
R = r-a with the dimension of length and rriph = m{k^/a and keep them fixed as k^ — >■ k^. 
This fixes 0(^4) ~ (^4 — kf)^ . In addition one can readily prove Fischer's scaling relation: 
7, = iy{2-r]). 

As usual one can prove that the exponent p is related to the Hausdorff dimension of 
the ensemble of manifolds. However, in this case we have in mind an intrinsic Hausdorff 
dimension dh- For a manifold with N4 4-simplexes the volume ~ A^4 as long as we take 
a = 1. The average volume is therefore ~ (N^)^, where (■)^ refers to the ensemble T(2, r). 
We can now define the internal Hausdorff dimension dh by 

{N)^ ~ r'^'^ for r -^ 00, m{ki)r = const. (6.26) 

From (|6.21| ) we have 

/^^^ dhg Gk,{r) dmiki) 



dA;4 


dk^ 


and using (|6.25|) we get 




(iV)^~ri/^ i.e. 4 = 


-l/u. 



(6.28) 

As is seen from the above consideration one can apply without any problems general 
scaling considerations known from critical phenomena provided that we have the asymp- 
totic behavior (|6.19|) . It can be formulated even stronger: It is hard to imagine that any 
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continuum limit of a conventional kind can be taken unless that mass m{k4) -^ for 
^4 -^ kl{k2). 

A good test case could be to return to 2d gravity, but surprisingly it turns out that in{^) 
{n = k^ in d = 2) and the associated critical exponents have only recently been calculated 
even if they are the most fundamental scaling variables in two dimensional quantum 
gravity. Unfortunately there is no space for a detailed discussion of the calculation, which 
uses the so-called transform matrix formulation [ 80]. Details can be found in [ 81] and 
the surprising result is that the two-point function G^{r) can be calculated even at the 
discretized level: 



cosh 



:A/x)i/3r 



G^{r) = const. (A/i)^/" ^ ^ J (q_29) 



"c- 



sinh^ {Ajj)il3r 

where const, and (3 are positive constants of 0(1) and A/x = ^ — ^c 

We conclude the following: 

1 

(i) G^{r) falls of like e"^*^^^)'^** for r -^ oo, i.e. the critical exponent i^ = \ and the 
Hausdorff dimension (i/^ = 4. 



(ii) Gfj,{r) behaves like r "^ for 1 ^ r <^ A/i 4, i.e. the scaling exponent rj = 4. 

(iii) The pre-exponential factor to e"^*^^^)^'^'" is m(/i)^~^, and this factor is precisely 
the one needed if we should be able to define a continuum two-point function up to an 
allover scaling factor by the assignment (|6.26| ), which states that rm{fi) is constant when 



r — »• 00 and /i — > fie, since the short distance behavior is given by r^~ 

(iv) The exponent a giving the power correction to the exponential decay is zero. This 

can be understood intuitively from from the fact universes for r S> l/m(/i) have to be 

thin tubes, and it is not difficult to show that these have no correction to the exponential 

decay. 

(v) From Fisher's scaling relation we get 7^ = ^{2 — r]) = —1/2. This well known result 

can of course also be derived directly from 

00 
X{fj,) = Y, G^{r) = const. - c^iAfx)^ + ■■■ (6.30) 

r=l 



by use of ( |6.29| ), but it should be clear that the explicit calculation in ( |6.30[ ) is nothing 



but a specific example of the general calculation used in proving Fisher's scaling relation. 
What is somewhat unusual compared to ordinary statistical systems is that the anomalous 
scaling dimension r] > 2. r^ = is the ordinary free field result, while 77 = 2 is the infinite 
temperature limit, and for statistical systems we expect rj < 2. 
(vi) The continuum two-point function can now be defined as 

G(fl;M) = to(VSr'G,(r)~M3||^. (6,31) 

where i? is a continuum distance and M a continuum mass proportional to A^/^, the 
unusual power of the "mass" being due to dh = 4. The factor in front of G^(r) is the 
usual "wave function renormalization" present in the path integral representation of the 
propagator. 
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It is interesting to give a direct physical interpretation of the short distance behavior 
of Gk^{r) as defined by ( |6.21| )E3 (in the following we suppress the index k2). In order to 
do so let us change from the grand canonical ensemble given by (|6.21| ) to the canonical 
ensemble defined by T(2, r, A^4), the class of triangulations with N^ 4-simplexes where two 
of the 4-simplexes are marked and separated by a distance r. On this ensemble we can 
define the discretized analogue of Gk4,{r): 

G{r,N,)= Y. 1' (6-32) 

Ter(2,r,7V4) 

For r = we have the following A^4 dependence (for the r dependence see (6.37) below) 

G(0, A^4) ~ NJ^-^ e'^=^^ (6.33) 

assuming that 7^ exists. This remark is important since it is by no means obvious that that 
7s exists in four-dimensional gravity, as discussed above. The partition functionZ(A^4) ~ 
A^^"" Qf^l^-i and the one-point function for large N^ is proportional to N4^Z{Ni) since 
it counts the triangulations with one marked simplex For r = (or just small) there is 
essentially no difference between the one-point function and G{Q,N^. 
G{r,N4) is related to Gk4,{r) by a (discrete) Laplace transform: 

Gk,{r)=Y.G{r,N,)e-''^''^. (6.34) 

N4 

The long distance behavior of G{r, N4) is determined by the long distance behavior of 
Gk4,{r). Close to the scaling limit it follows by direct calculation (e.g. a saddle point 
calculationEa) that 

G'fc,(r)~e-^'('^*-'=4)^/'^^ ^ 

G(r,iV4) ~e-<^'''/^^)'^ e*^4A^^ for r^'^ > N^, (6.35) 

wherec=(4-l)/rfj^^^'^'"'^. 

On the other hand the short distance behavior of Gk4,{r) is determined by the short 
distance behavior of G{r, N4) which is simple. Eqs. ( p.26| )- (|6.28| ) define the concept 



of Hausdorff dimension in the grand canonical ensemble. A definition in the canonical 
ensemble would be: Take NJ '" 3> r and simply count the volume (here number of 4- 
simplexes) of a "spherical shell" of thickness 1 and radius r from a marked simplex, sum 
over all triangulations T]^^ with one marked 4-simplex and divide by the total number of 
such triangulations. Call this number {n{r))]y^. The Hausdorff dimension is then defined 

by 

(n(r))iV4 ~ r"^''"^ for 1 <^ r <^ nI^'^\ (6.36) 



^^ The remarks to follows are of course equally valid for G^ (r) in two-dimensional gravity. 
^^ In addition to the exponentially decaying part of G{r, N4) there is also a power correction coming from 
the quadratic integration in the saddle point approximation. We shall not consider the explicit form of 
the power correction here. 



It follows from the definitions that we can write 



G{r,Ni)r^r'^''-^N2~^e''i^^ for I <^ r ^ nI^'^\ (6.37) 



We can finally calculate the short distance behavior of G'^(r) from eq. ( |6.34|) . For A;4 — > k 
we get: 



Gfc.W -r^'^-i^iVr' e-<^'''/^0'^ -r^'^"-! (6.38) 

This shows that 

r] = 2--fsdh, i.e. 75 = 1/(2-77), (6.39) 

which is Fishers scaling relation. The above arguments highlight that the anomalous 
scaling dimension r^ is a function of the two kinds of fractal structures we can define 
on the ensemble of piecewise linear manifolds: the Hausdorff dimension and the baby 
universe proliferation probability. In addition the arguments show that the canonical and 
grand canonical definitions of Hausdorff dimension agree. 

6.4- Branched polymers 

The model of branched polymers [BP) provides us with a simple, but non-trivial example 
of the above scenario [ 61] and will play an important role in the following. In a certain 
way it can be viewed as the lowest dimensional fractal structure and it will appear as the 
limiting case of higher dimensional gravity theories. 

Let us define branched polymers as the sum over all tree graphs (no loops in the 
graphs) with certain weights given to the graphs according to the following definition of 
the partition function: 

Z(/i) = 5: -^ p(5P) e-^l^^l, (6.40) 

where \BP\ is the number of links in a BP and /x is a chemical potential for the number 
of links, while 

piBP) = n fin^), (6.41) 

if^BP 

where i denotes a vertex, rij the number of links joining at vertex i and /(nj) is non- 
negative, /(rij) can be viewed as the unnormalized branching weight for one link branching 
into Ui — l links at vertex i. Finally Cbp is a symmetry factor such that rooted branched 
polymers, i.e. polymers with the first link marked, is counted only once. 

This model can be solved [7]. It has a critical point fie (depending on /) and close to 
the critical point we have: 

Z"(/i)~(/x-/i,)-l/^ (6.42) 
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Fig. 24. The graphical representation of the two-point function for branched polymers. The dashed 
line represents the unique shortest path between the two marked vertices. The "blobs" represent the 
contribution from all rooted polymers branching out from a vertex. 

i.e. 7s = 1/2 for branched polymers. On the branched polymers we define the "geodesic 
distance" between two vertices as the shortest link path, which is unique since we consider 
tree-graphs. The graphic representation of the two-point function is shown in fig. 24. Had 
it not been for the ability to branch, the two-point function would simply be 

G^(r)=e-'^^ (6.43) 

However, the insertion of one-point functions at any vertex leads to a non-analytic coupling 
constant renormalization and the result is changed to [61] 

G^(r) = const, e'''^^^^^^ for /x - /x^ ^ 0, (6.44) 

where k is a positive constant depending on /. We can now find G{r, N) by an inverse 
Laplace transform: 

G(r, N) = const. A^-^/V e-^'"-'/^^ . (6.45) 

We confirm from this explicit expression that the (internal) Hausdorff dimension of 
branched polymers is 2 (like a smooth surface !) and that 7^ = 1/2 since the prefac- 
tor of G{r, N) for small r should be ]\["/=-'^r'^H-i^ 

6.5. Numerical simulations 

Presently it has not been possible to make much progress in analyzing (|6.13| ) by analytical 
methods except for d = 2. However, the action is well suited for the use of Monte Carlo 
simulations and the results to be discussed later for d = 3 and d = 4 have been obtained 
by such simulations. Let me shortly discuss the principles involved in such simulations 
since there are a number of interesting aspects involved compared to ordinary lattice 
simulations. 

Monte Carlo simulations usually operate by means of a stochastic process, usually taken 
to be a Markov chain. The chain is selected such that it has a stationary distribution equal 
to the desired probability distribution, in this case the probability distribution given by 
( [6.131 ). The Markov chain is a prescription for moving around in the configuration space, 
each step being independent of the former steps and chosen with a certain probability 
among a set of possible steps. Two conditions exist which, if fulfilled, are sufficient to 

90 





Fig. 25. The moves in d = 2 and d ~ 3. 

ensure that the chain converges to the desired probabihty distribution. The first condition 
is ergodicity, i.e. by applying the available steps successively one should be able to move 
between any two configurations. The other condition is detailed balance. Let the steps be 
defined by a probability P{A -^ B) for going from state A to state B and let the desired 
probability distribution be given by e"'^^"^). The equation of detailed balance state: 



-S{A) 



P{A -^ B) 



~S{B) 



P{B -^ A). 



(6.46) 



Let us first discuss the question of ergodicity. Given the class combinatorial triangula- 
tions of a manifold we want to find a sequence of so-called moves which are ergodic. Such 
moves have been known for a long time and they are known as the Alexander moves [ 82]. 
They can be simplified somewhat and in d dimensions there are d+1 oi them [ 83]. They 
can be described as follows: Given a sub-simplex i of order d+l-i the (i+l-i rf-simplexes 
which share i form a d-ball B{i). Remove i and all the higher dimensional simplexes to 
which i belongs. Instead insert an "orthogonal" d — i dimensional sub-simplex, i.e. the 
sub-simplex constructed from the d—i vertices oidB{i) which did not belong to i, together 
with all the higher dimensional sub-simplexes. In this way B{i) has been exchanged with 
a new rf-ball B'{d — i) with the same boundary: dB{i) = dB'{d — i). This is illustrated 
for d = 2 and d = 3 in fig. 25. 

This set of moves is ergodic for a given combinatorial manifold. In two dimensions it is 
seen that one of the moves preserves N2, the number of triangles in the triangulation. In 
fact it can be shown that this move alone is ergodic on the set of triangulations which a 
fixed topology and and fixed volume (i.e. A^2)- From a computational point of view it is 
very convenient to be able to keep the volume fixed. Can it be done in higher dimensions 
with a suitable set of moves of the same local nature as the ones mentioned above? The 
answer is no, at least in four- and higher dimensions, the reason being that it is known 
that in these dimensions manifolds exist which are not algorithmic recognizable. More 
precisely this means that manifolds TWq exist such that no finite algorithm in the sense of 
Turing allows us to decide if an arbitrarily chosen manifold is combinatorial equivalent to 
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the given Mq. //there were ergodic moves which kept the volume fixed one can argue that 
it takes only a finite number of operations, computable as a function of the volume, to list 
all combinatorial manifolds equivalent with A^o? and for a given combinatorial manifold 
one can now by inspection check if it is in the list. Again it can be argued that this only 
takes a finite number of operations, computable as a function of the volume [ 74] . 

It is seen that there are deep reasons which makes it impossible to have a local set 
of moves which preserves the volume and are ergodic for all manifolds. The situation is 
even more drastic: Let A^ be the volume and let us ask the question: what is the smallest 
number of moves needed in order to connect any two configuration with volume N. The 
answer is that for manifolds like Aio this number cannot be bounded by any recursively 
definable function. This result is only possible if we on the way from one configuration 
to another with volume N has to make large detours to configurations with a volume 
N' » N [ 74]. 

Let us finally discuss the implementation of detailed balance: The transition from a 
configuration A to a neighboring one B can be realized in two steps. First we pick an 
i-simplex shared hj d + 1 — i d-simplexes. The probability of doing this is l/ni{A), ni{A) 
being the number of i-simplexes in A. Afterwards we exchange i with an d — i simplex 
with some probability Pi{A -^ B). The equation for detailed balance reads: 

where i = 0,1,. ..,d. One can now choose Pj's according to some standard Metropolis 
algorithm, but it is worth to emphasize the appearance of the combinatorial prefactors 
l/rii which are usually absent in ordinary lattice simulations in the sense that they cancel 
out. They appear here because we have a dynamical lattice which change during the 
updating. 

6.6. Results 

Let me now discuss the results of the numerical simulations. They can be summarized as 
follows [ 90, 91, 94, 96, 95, 92, 93] 

(1): There seem to be two different regions, as a function of the bare inverse gravitational 
coupling constant k2'. For small or negative values of /c2 the typical quantum universe will 
be very crumpled, with almost no extension and a very large, if not infinite Hausdorff 
dimension, while the universes for large values of /c2 will be elongated with a Hausdorff 
dimension as small as two. In fig. 26 we have shown the average radius for universes of 
size 9000, 16000, 32000 and 64000 4-simplexes as a function of k2. The two phases are 
separated by a phase transition which is of order two or higher. At the transition point, 
^2, the Hausdorff dimension might be finite (the precise value is not well determined, but 
it could be close to four). 

(2): The same results are valid for three dimensional simplicial quantum gravity except 
that the phase transition seems to be of first order, rather than of higher order [ 84, 85, 
86, 87]. 

From fig. 26 it is seen that the change between the elongated region and the crumpled 
region becomes increasing visible as the size of the system increases. In addition the 
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Fig. 26. The average radius of the universes of sizes 9000 (triangles), 16000 (squares), 32000 (pentagons) 
and 64.000 (stars) versus k2- 
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Fig. 27. The pseudo critical point k2{Ni). 

critical point seems to move to higher values of k2- We observe a so called pseudo critical 
point k2{N^. The indication of convergence to a limiting value ^2 ^^ the volume iV4 ^ oo 
is shown in fig. 27 and we conclude that we have a genuine phase transition. From the 
general theory of finite size scaling the pseudo critical point A;2(^4) of 3' first order phase 
transition scales to the limiting value fcg as I/N4. Since we observe a scaling like l/^/W^ we 
conclude that the transition most likely is a higher order transition. In three dimensional 
simplicial quantum gravity it is difficult to perform this kind of measurement since we 
observe pronounce hysteresis around the transition point. This is a typical sign of a first 
order transition. 

Let us first discuss the measurements in the elongated phase, i.e. for k2 > k"^- An ideal 
quantity to measure in the computer simulations is the number of four-simplexes (ji{r))j^^ 
within a spherical shell of thickness 1, a geodesic distance r from a marked four simplex, 
the average taken over all spherical triangulations with one marked four-simplex. It will 
depend on the coupling constant k2 and it is related to the two-point function precisely 
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Fig. 28. The measured distribution (dots, the error bars smaUer than the dots), and the fit (the curve) 
using the functional form r exp{—r'^ / N4) . 



as discussed above: 

G{r,N^) ~ Ar;('=2)-2gfc|(fc2)A^4 (n(r))^„ 



(6.48) 



provided the subleading corrections are power-like. In fig. 28 we show the measured distri- 
bution ( n{r)/N4)jq^ just above the transition to the elongated phase and a corresponding 
fit to rexp(— r^/A^4). We see a very good agreement and it becomes even better if we 
move further into the elongated phase. Now recall the general scaling relations 



{n{r)) 
{n{r)) 



N4 



^^ f 



dh~l 



for 



l<r<Arl/'^'\ 



N4' 



a -c{r''h/N4)'^/('ih--i-) 



-re 



for nI^'^'' <r <N^, 



(6.49) 
(6.50) 



we conclude that dh = 2 (and a = 1). In addition we can measure the critical exponent 
7s very conveniently by baby universe counting. Again the arguments are identical to 
ones presented in the two-dimensional case. The result is shown in fig. 29 and it is natural 
from the figure to conjecture that 7 = 1/2 as long as we are in the elongated phase. 
From ( |6.48| ) and {n{r))]sr^ we know G{r,N4) and we can construct the two-point function 
Gk4{r; /C2)) by Laplace transform (as in the two-dimensional case): 



Gk4{r;k2 



00 

E 

Af4 = l 



N-^/h 



-A{k4)N4 ^g-crVA^4 _g- 



■cr\/Ak4 



(6.51) 



where A(/c4) = ^4 — kl{k2) is assumed to be small. This function is precisely the two- 
point function of the so called branched polymers, which are known to have internal 
Hausdorff dimension dh = 2. We conclude that the numerical simulations provide con- 
vincing evidence that the elongated phase of simplicial quantum gravity corresponds to 
a well known statistical theory, the one of branched polymers. The tendency to create 
baby universes is so pronounced in this phase that the geometry degenerates to the generic 
lowest dimensional fractal structure possible, i.e. that of branched polymers. 
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Fig. 29. The measure 7 in the elongated phase for various size lattices (iV4 — 9000, 16000, 32000 and 
64000). 

When we lower the value of k2 and move below the critical point ^2 the fractal structure of 
our ensemble of piecewise linear manifolds changes drastically. A glance on fig. 26 shows 
that the average radius hardly changes with the volume. This is an indication that the 
HausdorfF dimension is large or maybe infinite. If we move deep into this phase the average 
curvature is negative and in addition there are only few baby universes and they are small. 
This could lead to the idea that we entered a phase with "smooth" manifolds of negative 
curvature. For such manifolds one would expect that the volume of geodesic balls of radius 
r would grow exponentially with the radius, which is what we observe. Clearly this is a 
"fake" infinite Hausdorff dimension and indeed we should observe the dimension dh = 4: 
in the sense that {n{r))N^ ~ r^ for small geodesic distances. A closer look at "typical" 
members of the computer generated manifolds indicates that they cannot be considered 
as "smooth" . Rather they have a few vertices of very high orders which connect to almost 
any other vertex in the manifold and in such a situation it is not surprising that the linear 
extension will be small. In addition we have not been able to fit to any sensible power 
dependence r'^H-^ (qj- gniall r. A plot of log{n{r)) j^i^ shows indeed a linearly growing 
function of r up to some ro(A^4) which is not much different from the average value {r)^^. 
A reasonable fit to {r)N4 is 



{r)N, = a{k2) + b{k2)logNi 



(6.52) 



This again gives some support to the idea that the Hausdorff dimension is infinite in this 
phase since it appears as a limit of {r)N4 ~ N^^'^" for rf/^ — > 00. Finally the extrapolation 
of (6.50) to infinite Hausdorff dimension indicates that in such a case one should observe 



{n{r)) 



N4 ~ ci e" 



-mi(fc2)r 



(6.53) 



down to quite small distances r ^ A^4 , dh -^ 00, which naturally is replaced by 
r > Oi logiV4 + 02- This is indeed what we measure. 

The observation that {n{r))N^ grows exponentially from r ^ 6 out to r ^ ro and then 
falls off exponentially indicates that we deal with an infinite Hausdorff dimension at all 
distances and it is easy to get a quite good "phenomenological" fit to {n{r))N4 which 
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Fig. 30. exponential fit (curve) of the form (6.54) to the measured n(r;iV4) (dots, error bars less than 



dot-size) in the crumpled phase (iV4 = 64000, fc2 = 1.26). 
incorporates both these features by choosing e.g.: 

{n{r))N, ~ exp (-mi(fc2)r - cse-™^^'^^^) . 

It will grow like e('^2m2-mi)r ^^^ small distances and fall off like 



(6.54) 



(6.55) 



for large distances, while a N4 dependence in the coefficient C2 would explain the observed 
N4 dependence of tq. The data and a fit of the form ( |6.54D are shown in fig. 30 for /c2 = 1-26 
and A^4 = 64000, i.e. right below the transition to the crumpled phase, where the fit is 
worst. But even so close to critical point ( |6.54 ) works quite well over the whole range of 
r > 6. It should be mentioned that the coefficient in front of the second exponential in 
eq. (|6.55|) is negative. This implies that the term cannot be given the interpretation as 
an additional heavier mass excitation. However, just looking at the long distance tail the 
distribution {n{r))N4 allows us to determine mi, 7112 and C2 from ( |6.55|) . On the other hand 
we can determine 02/71.2 — '^i from the short distance exponential growth alone and find 
good agreement. This indicates that long and short distance behavior are intertwined 
in the case of infinite Hausdorff dimension, as they are in the case of finite Hausdorff 
dimension, where dn appears both in the short distance and long distance expression for 
{n{r))N, (see e.g. (6.49)-(6.50)). 

We conjecture that the internal Hausdorff dimension is infinite for ^2 < ^2- 

In accordance which this conjecture we have in the elongated phase no "mass" term 
m{ki) = (^4 — k'l{k2)Y which scales to zero as we approach ^4(^2)- The exponential 
coefficient mi{k2) is finite in the infinite volume limit ^4 -^ ^4(^2)- However, it is most 
interesting that mi{k2) scales to zero as ^2 ~^ ^2' ^^^ phase transition point between the 
crumpled and the elongated phase. This gives a strong indication that the system at the 
transition might have a finite Hausdorff' dimension, which could very well be larger that 
the generic dh = 2 found in the elongated phase. In fig. 31 we have shown the scaling of 
the mass mi{k2) as a function of /c2- 
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Fig. 31. The behavior of the masses mi{k2) and mi{k2) + '712(^2) from (6.54) in the crumpled phase as 
a function of kn — ko. 



The above mentioned numerical "experiments" suggest the following scenario: The typical 
quantum universe, determined without any Einstein action (i.e. ^2 = 0) has (almost) no 
extension. Its Hausdorff dimension might be infinite and internal distances between points 
are always "at the Planck scale" . By that we simply mean that no consistent scaling can 
be found which will be compatible with finite continuum volume and finite Hausdorff 
dimension. For a finite value of the bare gravitational coupling constant there is a phase 
transition (second or higher order) to a phase with a completely different geometry with 
pronounced fractal structures. It is tempting to view the transition between the two 
kind of geometries as a transition where excitations related to the conformal mode are 
liberated, since large /c2 is a region which formally corresponds to small values of the 
gravitational coupling constant. Right at the transition it seems as if we have the chance 
to encounter genuine extended structures with a finite Hausdorff dimension. Maybe the 
fact that the transition between the two types of geometry is of second (or higher) order 
can be used as the starting point for a non-perturbative definition of quantum gravity. 

Of course it is crucial to be able to perform high statistics simulations at the critical 
point in order to investigate this possibility in greater detail. It would be even better if we 
could perform some analytical calculations. The branched polymer picture indicates that 
mean field methods should be available in some regions of the coupling constant space. 



7. Discussion 

We have seen in some detail that it is possible to discretize reparametrization invariant 
theories and apply with success the methods known from the theory of critical phenomena. 
In this way we deal with the theories of fiuctuating geometries. However, we did not 
really answer the most interesting question: How wildly should we allow the geometries 
to fluctuate? For a fixed topology it was possible to formulate a regularized Euclidean 
quantum theory. In two dimensions one can take the scaling limit of the regularized theory 
and the corresponding quantum theory correctly describes the interaction between matter 
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and gravity. In higher dimensions it is not yet known if it leads to any interesting theory, 
but at least there is a well defined non-perturbative procedure for how to investigate 
this question. In addition the discretized action in this approach is remarkable simple 
and one could hope that it will be possible to solve the theory in the same detail as 
in two dimensions. At the moment we let loose topology we have also lost control of 
the theory. In two dimensions the double scaling limit gave us a hint that it might be 
possible to perform the summation over topologies. At least there are some prescriptions 
for how it should be done. They are not yet unambiguous and this refiects our lack 
of understanding of the physics which goes beyond a simple perturbative expansion in 
topology. Nevertheless I view it as very encouraging that one at least seems to have by 
now the tools which allow such questions to be asked. These tools have their origin in the 
discretized approach and even in three and four dimensions the discretized approach has 
the advantage that it in a natural way is able to combine the summation over Riemannian 
structures and topology with the volume of the universe acting as the cut-off. It remains 
to be seen if the approach will result in an interesting theory of gravity, but it offers 
at least a playground for a fascinating interplay between pure mathematics, theoretical 
physics and computational physics. 
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